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1.0  INTRODUCTION  AND  BACKGROUND 


1.1  Introduction 

A  parasitic  array  consists  of  a  main  element  directly  connected 
to  a  receiver  plus  additional  parasitic  elements  coupled  to  the  receiving 
element  by  the  mutual  impedance  between  the  antennas.  The  array  antenna 
pattern  is  controlled  by  adjusting  the  values  of  the  complex  impedances 
that  terminate  the  parasitic  elements.  The  value  of  each  terminating  im¬ 
pedance  determines  the  phase  and  amplitude  of  the  signals  reflected  from 
that  parasitic  element  to  the  receiving  element.  Thus,  the  main  antenna 
element  receives  the  sum  of  the  various  reflected  signals  to  form  an  array 
output.  The  use  of  such  an  array  for  directional  beam  forming  has  been 
studied  by  Harrington,  [6-10]. 

A  spatial  interference  canceller  (SIC)  uses  a  parasitic  array 
with  complex  terminations  for  directional  null  forming.  In  addition,  an 
adaptive  SIC  would  require  an  automatic  control  algorithm  for  controlling 
the  complex  terminations,  such  as:  least  mean  square  (LMS) ,  gradient  con¬ 
trol,  univariate  search,  dither,  etc. 

Zeger-Abrams  Incorporated  (ZAI)  has  performed  an  analysis  of  the 
use  of  a  parasitic  array  for  adaptive  nulling  [1,2].  Below,  the  "trans¬ 
mission-line"  model  therein  presented  is  refined,  and  an  extended  analysis 
of  the  methods  of  implementation  and  the  performance  of  using  a  parasitic 
array  for  adaptive  nulling  is  given. 

The  goal  of  the  SIC  program  is  to  develop  a  compact  adaptive 
array  (an  order  of  magnitude  or  more  smaller  than  a  conventional  array) 
using  parasitic  elements  with  variable  complex  terminations.  Such  an  ar¬ 
ray  would  provide  anti-jamming  (AJ)  capabilities  in  a  very  tiny  aperture, 
using  less  electronic  hardware  and  fewer  RF  cables. 

1.2  Background 

Many  present  and  future  receivers  (for  radio  communications, 
radar,  navigation,  IFF,  Dr,  ELINT,  etc.)  will  be  required  to  operate  In 
a  severe  electromagnetic  interference  (EMI)  environment  that  includes 
unintentional  local  and  deliberate  remote  sources  of  EMI. 

A  common  problem  in  Navy  HF  communications  is  the  presence  of 
strong  intentional  and  unintentional  sources  of  interference.  Usually 
these  sources  have  unpredictable  azimuthal  directions  of  arrival  which 
are  different  than  the  desired  signal  azimuths.  One  technique  for  reducing 


the  interference  is  to  use  an  adaptive  array  (AA)  that  is  capable  of  placing 
nulls  in  the  directions  of  interferers  and  a  main  lobe  in  the  direction  of 
the  required  signal. 

Adaptive  array  hardware  has  been  demonstrated  at  various  fre¬ 
quencies  from  VLF  to  C-band.  Previous  AA  designs  for  remotely  generated 
interferences  (e.g.,  ECM)  have  employed  an  array  of  antennas  spaced  greater 
than  a  half  wavelength  in  order  to  minimize  inter-element  coupling.  Inter¬ 
ference  cancellation  has  traditionally  taken  place  in  a  physical  device 
(directional  coupler,  transformer,  hybrid  summer)  functioning  as  a 
phasor  sub tractor. 

Prior  AA  designs  suffer  from  costly  installation  of  coaxial 
cables  to  elements,  poor  nulling  performance  due  to  inter-element  coupling 
effects,  and  great  consumption  of  available  space  on  aircraft  surfaces 
and  ship  masts  by  antennas. 

These  drawbacks  of  prior  designs  can  be  reduced  or  eliminated 
by  using  the  linear  superposition  property  of  free  space  (for  electro¬ 
magnetic  fields) -to  accomplish  spatial  interference  cancellation  (SIC). 

An  antenna  array  and  a  conventional  beamformer  are  illustrated 
in  Figure  la.  The  complex  signal  envelope  y^  incident  upon  the  nth  antenna 
is  weighted  (complex  multiplication)  by  the  complex  weight  wq.  With  an 
unweighted  main  antenna,  the  array  output  equals  yQ  +  Ew^y^,  which  is  the 
input  to  the  receiver.  Antenna  element  spacing  is  typically  X/2,  and  inter- 
element  coupling  is  small  and  is  usually  neglected. 

In  a  spacial  beamformer,  illustrated  in  Figure  lb,  a  main  antenna 
is  connected  directly  to  the  receiver,  and  a  compact  array  of  parasitic 
elements  is  terminated  in  complex  impedances.  The  value  of  the  terminating 
Impedance  z^  of  the  nth  parasitic  element  is  controlled  to  adjust  its 
complex  reflectivity  pQ.  In  this  way,  the  signal  y^  incident  upon  z,^  is 
weighted  before  being  coupled  to  the  main  antenna,  where  it  is  specially 
summed  with  contributions  from  the  other  parasitic  elements.  The  input  to 
the  receiver  is  equal  to  yn  +  Im_  p  y  ,  where  m.  is  the  dimensionless 
complex  coupling  coefficient  between  the  nth  parasitic  element  and  the  main 
antenna.  The  similarity  of  receiver  inputs  for  the  conventional  and  para¬ 
sitic  arrays  is  evident  if  one  substitutes  p  for  w  . 

On  n  n 

Array  pattern  control  in  the  conventional  array  is  achieved  by 
adjusting  the  wq.  In  the  parasitic  array  one  can  adjust  the  Pn  or  the 
If  the  complex  impedance  of  the  transmission  line  (between  the  nth  antenna 
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Figure  1  Adaptive  Beamformers 


By  combining  Equations  (2.6)  and  (2.7),  we  obtain  v_  in  terms  of  v 

“T  — oc 


v  -  Z  (Z  +  Z„)  1 v 
-^T  TVT  A  — oc 


(2.8) 


Thus  we  can  completely  characterize  the  array  during  reception 
knowing  the  open  circuit  mutual  impedance  matrices  of  the  antenna  and 
the  load  networks  and  the  vector  of  array  open  circuit  voltages,  vqc» 
due  to  an  incident  plane  wave.  Harrington  gives  computer  programs  which 

will  calculate  Z.  and  v  for  an  array  of  wire  dipole  antennas  [10,12]. 

A  — oc 

The  impedance  model  has  been  briefly  introduced  here  to  allow 
the  reader  to  become  aquainted  with  it.  The  impedance  model  is  fairly 
well  defined  in  the  literature.  In  the  next  chapter,  the  transmission 
line  model  that  Zeger-Abrams  Incorporated  has  developed  shall  be  related 
to  the  impedance  model.  This  relationship  is  based  on  the  theory  of 
scattering  parameters,  of  which  a  derivation  is  given  in  Appendix  B. 

2.3  Properties  of  Compact  Parasitic  Arrays 

The  compact  parasitic  array  has  many  differences  with  the  con¬ 
ventional  antenna  array.  These  differences  will  be  pointed  out  and  dis¬ 
cussed  in  this  section. 

The  most  important  difference  is  that  the  compact  parasitic 
array  has  strong  mutual  coupling.  In  the  conventional  array,  mutual 
coupling  is  undesired  and  efforts  are  made  to  keep  it  minimal.  The  prin¬ 
cipal  manner  in  which  this  is  done,  is  to  keep  the  antennas  far  apart. 

In  the  parasitic  arrays,  strong  mutual  coupling  is  necessary  for  oper¬ 
ation.  A  small  compact  array  (dimensions  on  the  order  of  0.1X  or  less) 
are  needed  to  ensure  strong  mutual  coupling. 

Because  the  array  is  small,  many  features  of  the  array  will 
be  broadband  and  broadbeamed.  High  directivity  will  not  be  possible  as 
the  main  lobe  will  always  be  somewhat  broad.  Similarly,  when  nulls  are 
formed,  the  null  widths  will  be  slightly  wider  than  those  obtained  in 
conventional  arrays.  For  details,  see  the  experimental  results  reported 
in  Chapters  6  and  8. 

The  frequency  characteristics  of  the  array,  particularly  with 
active  terminations,  are  difficult  to  determine.  The  reflections  off 
the  parasitic  element  terminations  and  the  mutual  coupling  create 
very  complex  feedback  that  is  difficult  to  analyze.  Even  more  funda- 


The  process  of  nulling  a  signal  is  one  of  finding  Che  PQ's  such 
Chat  y^  ■  0  in  Equation  (2.3a).  This  is  done  in  Section  4.1.  The 
resulting  Pn's  turn  out  to  be  independent  of  p^.  This  results  from  the 
fact  when  there  is  no  incident  voltage  to  reflect  at  the  termination 
(receiver)  of  the  main  antenna,  it  doesn't  matter  to  the  rest  of  the  sys¬ 
tem  what  Pq  is. 

A  final  comment,  for  the  transmission  line  model  presented  here 
it  has  been  assumed  that  there  are  transmission  lines  between  the  antennas 
and  the  terminations.  It  turns  out  that  this  is  not  necessary  (or  more 
specifically,  one  can  assume  transmission  lines  of  zero  length) .  This  is 
shown  in  Appendices  C.l  and  C.2. 

2.2  Impedance  Model 

The  impedance  model  has  been  used  by  Harrington  in  several  papers 
to  describe  a  parasitic  array  for  various  purposes  [6-10].  This  model  of 
an  array  while  receiving  is  illustrated  in  Figure  5.  is  the  open  cir¬ 
cuit  mutual  impedance  matrix  of  the  antenna  array  looking  directly  into 
the  antenna  ports.  is  the  open  circuit  mutual  impedance  matrix  looking 

at  the  terminations  through  the  transmission  lines.  In  our  array  Z T  is 
diagonal,  but  the  main  results  are  valid  for  a  general  terminating  network. 

The  open  circuit  voltage,  VQcn,  is  the  terminal  voltage  at  the  nth  port  due 
to  a  plane  wave  Incident  on  the  array  with  all  antenna  ports  open  circuited. 

As  with  of  the  transmission  line  model,  v^  will  depend  on  the  field  strength 
and  wave  vector,  k,  of  the  incident  plane  wave.  For  convenience,  the  model 
used  here  assumes  there  are  no  transmission  lines  in  between  the  terminations 
and  the  array  (or  at  least  the  effect  of  these  transmission  lines  are  included 
in  the  impedance  parameters) . 

To  compute  the  vector  of  voltages  across  the  load  terminals,  v^, 
we  must  first  find  the  current  vector,  i.. 


“DC 


=  -<ZT  +  ZJt 


therefore 


i  -  -(Z  +  Z)~ 1  v 

T  A  — oc 


But  we  know  that  i  is  related  to  vT  by  the  load  matrix  ZT> 

-Z_i 


JT— 


(2.6) 


(2.7) 
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If  the  Input  impedance  of  the  receiver  attached  to  the  main 
tches 

(2.3)  reduce  to 


antenna  matches  the  line  impedance  ■  z  »  then  prt  ■  0.  The  Equations 

o  u 


ys0  +  ^0  PZ 


+  MPv 


(2.4a) 

(2.4b) 


(with  pQ  =  0) 

Equation  (2.4b)  can  be  solved  for  the  incident  voltage  vector 
of  the  auxiliaries,  z» 

Z  =  (I  “  MP)“  ^ 


Thus  Equations  (2.4)  can  be  rewritten  as 

y0  =  ya0+*0**\ 

i  -  Q^Zg 


(2.5a) 

(2.5b) 


(with  pQ  -  0  and  Q  -  I  -  MP) 


Equations  (2.5)  could  have  been  obtained  directly  from  Equation 
(2.2)  using  matrix  identity  (A. 9)  from  Appendix  A  with  the  fact  pQ  -  0. 


Q+  (PQ  -  0) 


T 

-So  p 


I  -  MP, 


From  identity  A. 9  we  can  compute 


-I 


(P0  =  0) 


"  T  -1 
1  :»o  pq 


o  :  q 


rl 


thus 


which  is  equivalent  to  Equations  (2.5). 
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Figure  4  Illustration  ol 
of  the  Transmi! 


Equation 
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Let  us  now  examine  what  happens  during  reception.  The  vector 
of  voltages,  y_j_,  incident  on  the  terminations  will  be  the  sum  of  two  terms, 
and  M+Ppr+,  as  given  in  Equation  (2.1).  (See  Figure  4.) 

*+  “  *s+  +  M+P+?+  (2.!) 


The  term  (  is  due  to  the  incident  plane  wave.  The  term  M+P+y+  is  due 
to  the  incident  voltages  of  being  reflected  off  the  terminations  as 
and  then  coupled  via  the  mutual  coupling  matrix  M+  into  downward 
traveling  voltages  of  M+Ppr. . 

Equation  (2.1)  demonstrates  that  the  system  is  one  of  multivar¬ 
iate  linear  feedback  in  terms  of  the  Incident  voltages  £+.  The  input  is 
yfl  | .  The  "loop  gain"  is  M+P+.  Since  active  terminations  are  being  allowed 
the  reflectivity  of  the  matrix  P+  may  exhibit  gain.  Thus  it  is  reasonable 
to  ask  when  does  the  matrix  product  H.P ,  represent  positive  feedback  which 

T*  T 

renders  the  system  unstable.  This  issue  of  stability  is  addressed  in  Sec¬ 
tion  4.2. 

Given  that  the  system  is  stable,  one  may  solve  Equation  (2.1) 
for  the  voltages  incident  on  the  terminals,  y+. 


<1  -  M+P+)"lI,+ 


where  Q+  -  I  -  M+P+ 


(2.2) 


We 

for  the  main 


can  also  separate  Equation  (2.1)  into  two  equations, 
antenna  and  one  for  the  parasitic  auxiliaries. 


one 


y 
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yS0  +  P0a00y0  +-5oTp2. 

+  Wo  + 


(2.3a) 
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(2.3b) 
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Figure  3  Measurement  of  Transmission  Line  Model  Parameters 


Figure 


The  systematic  approach  considered  the  voltages  In  the  array 
after  all  the  "ping-ponging"  has  finished.  The  equations  that  described 
this  state  were  written  down  and  solved  for  the  various  traveling  voltage 
waves.  In  matrix  form,  these  equations  are  represented  by  Equation  (2.1) 
below.  In  the  initial  analysis  the  systematic  approach  yielded  little 
insight  as  to  what  values  of  the  reflectivities,  Pq,  and  mutual  coupling 
coefficients,  j •  allow  stable  operation.  The  only  restriction  that 
came  out  of  the  analysis  was 

det  |Q+|  -  det  |I  -  M+P+|  t  0  . 

This  constraint  was  not  restrictive  enough.  It  did  not  prohibit  all  sit¬ 
uations  that  were  inoperable.  In  particular  it  did  not  prohibit  the  equival¬ 
ent  of  positive  feedback. 

In  the  analysis  that  follows  the  systematic  approach  is  refined. 

With  this  new  model  several  issues  are  explored.  These  issues  include 
stability,  null  forming,  adaptive  control,  and  the  relationship  to  Harrington' 
impedance  model. 

The  refined  transmission  line  model  and  its  parameters  are  il¬ 
lustrated  in  Figure  2(a-d).  Note  it  has  been  assumed  M+  is  symmetric,  or 
in  general  m^  ■  m_.^,  This  is  true,  given  that  all  the  transmission  lines 
connecting  the  antennas  with  the  terminations  have  the  same  characteristic 
impedance.  This  is  shown  in  Appendix  C.2,  when  m^  are  calculated 
directly  from  their  definition  via  the  antenna  mutual  impedance  matrix,  Z, . 
Hereafter,  it  is  assumed  that  all  these  characteristic  impedances  are  the  same 

The  parameters  of  the  transmission  line  model  have  been  defined 
in  a  manner  that  allows  them  to  be  measured  relatively  easily  from  a 
physical  array.  Configurations  for  measuring  the  entries  of  M+  and  y^  t 
are  given  in  Figures  3a  and  3b  respectively.  Note  adjustments  in  the  cal¬ 
culations  will  have  to  be  made  due  to  the  phase  and  amplitude  effects  of 
the  directional  couplers  and  the  extra  cable  lengths. 

In  Sections  3.3  and  3.4,  formulae  for  the  computation  of  M+  and 
in  terms  of  the  open  circuit  mutual  impedance  matrix  of  the  array,  Z^, 
are  derived.  However,  calculation  or  measurement  of  Z^  is  quite  compli¬ 
cated.  Indeed,  for  measuring  Z^,  it  is  probably  easier  to  measure  M+  and 
then  calculate  Z.  from  it. 


2.0  CHARACTERIZATION  OF  PARASITIC  ANTENNA  ARRAYS 


The  two  models  of  a  parasitic  array  are  described  in  this  chapter. 
Each  model  will  be  described  in  terms  of  a  receiving  array.  Then  in  Chapter 
3,  the  relationships  between  the  models  will  be  established  using  both  recelv 
ing  and  transmitting  arrays.  The  models  have  about  equal  complexity  when 
analyzing  a  receiving  array.  The  impedance  model  is  easier  to  analyze  for 
a  transmitting  array.  The  transmission  line  model  is  easier  to  analyze  for 
control  of  the  terminations  (particularly  active  terminations)  as  an  adaptive 
nulling  array .  The  transmission  line  model  assumes  transmission  lines 
between  the  array  and  the  terminations,  although  these  lines  may  be  of  zero 
length.  The  impedance  model  may  be  used  with  transmission  lines  between  the 
terminations  and  the  array,  although  the  calculations  are  more  complex  (as 
one  sees  in  Chapter  3) . 

The  main  difference  between  a  transmitting  and  a  receiving  array 
is  in  a  transmitting  array  the  termination  is  the  source  and  the  array  is  the 
load,  while  in  a  receiving  array  the  array  is  the  source  and  the  termination 
is  the  load. 

2.1  Transmission  Line  Model 

The  transmission  line  model  of  the  array  analyzes  the  voltage 
traveling-waves  within  the  cables,  antennas  and  terminations,  the  inter¬ 
coupling  of  these  waves  through  the  antennas,  and  the  reflections  of  the 
waves  off  the  terminations. 

The  parasitic  array  with  complex  terminations  has  been  previ¬ 
ously  analyzed  by  ZAI  with  transmission  line  techniques,  [1,2].  In  this 
previous  analysis  two  analytic  approaches  were  used:  the  infinite  series 
(ping-pong)  approach  and  the  systematic  approach.  These  results  are  sum¬ 
marized  below. 

The  infinite  series  approach  considered  all  the  possible  reflec¬ 
tions  off  the  terminations  and  antennas  and  all  the  mutual  coupling  between 
antennas.  The  result  was  expressed  as  an  infinite  geometric  series  and 
summed.  However,  only  the  two  antenna  array  could  be  analyzed.  For  larger 
arrays  the  mathematics  become  to  unwieldy.  In  addition,  the  constraints 
on  the  reflectivities,  Pq,  and  the  mutual  coupling  factors,  m^,  necessary 
for  convergence  of  the  Infinite  series  are  too  restrictive.  They  disallow 
regions  where  safe  operation  is  possible.  In  particular,  negative  feed¬ 
back  of  greater  than  unity  gain.  The  general  question  of  stability  is 
very  involved  and  is  taken  up  later  in  Section  4.2. 


V 


vector  of  the  N  +  1  terminal  voltages  of  the  antenna  array 
vector  of  the  N  +  1  port  currents  of  the  antenna  array. 


i 

S  =  (N  +  1)  x  (N  +  1)  matrix  of  scattering  parameters  of  the 
N  +  1  port  antenna  network. 
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A  - 


yn  due  to  a  plane  wave  incident  on  the  array  with  all 
terminations  =  ZQ»  the  characteristic  impedance  of  the 
transmission  lines,  n  -  0,  1,...,  N.  ygn  is  dependent 
on  the  wave  vector,  of  the  incident  plane  wave. 


mutual  coupling  coefficients: 

the  ratio  of  the  downward  (toward  the  termination)  voltage 
traveling  wave  incident  on  the  jth  termination  to  the 
upward  (toward  the  antenna)  voltage  traveling  wave  trans¬ 
mitted  into  the  cable  from  the  ith  termination.  All  other 
upward  traveling  waves  being  zero.  All  terminations  “  zq. 
j  “  0.  !»•••»  N.  Note,  H  and  M+  are  symmetric  matrices. 


characteristic  impedance  of  transmission  lines  between 
antennas  and  terminations 

impedance  terminating  the  cable  from  the  nth  antenna, 
n  ®  0,  1,...,  N. 

\  0  \ 

\°  'W 

open  circuit  mutual  impedance  matrix  of  the  antenna  array 
at  the  antenna  ports  (zero  cable  lengths) . 


-:-!")■  N  l,3’0>l . N 

and  M  =  jm^j  i’  ^  =  N 


6)  The  absence  of  the  subscript  "+"  on  the  matrices  M,  P,  or  Q, 

or  the  vectors  or  will  denote  that  they  are  N  x  N  or  N  x  1 
matrices  or  vectors,  respectively.  They  will  be  composed  of 
entries  only  corresponding  to  the  auxiliary  antennas  (e.g.,  M  or 

7)  The  impedance  model  works  exclusively  with  the  total  array  of 

main  antenna  plus  auxiliary.  This  shall  be  understood  with  all 
its  associated  variables  (ZT,  Z.,  v,  S,  etc.)  without  the  sub¬ 
script  "+".  •  1  A 

8)  Subscripts  "R"  and  "I"  will  denote,  respectively,  the  real 
and  imaginary  components  of  the  designated  value. 

9)  The  superscripts  "T"  and  "*"  will  denote,  respectively,  transpose 
and  complex  conjugation. 

10)  In  the  derivations  that  are  devoted  to  the  single  element  cases, 
all  variables  become  scalars  and  will  be  denoted  by  the  appro¬ 
priate  lower  case,  non-underlined  letters.  Also,  all  indices 
will  be  dropped. 


the  wave  vector  of  a  plane  wave  "  \  ’  w^ere  ^  ^-s  ^*6 

wavelength  and  u  is  the  unit  direction  vector  of  the  plane 


reflection  coefficient  of  the  nth  termination* 
n  *  Of  1  f  •  •  •  f  N • 


Pi  0 


0  *P, 


the  complex  envelope  of  the  voltage  traveling  wave  incident 
on  the  nth  termination,  n  ■  0,  1,...,  N. 


(1.1) 


and  its  termination  z  )  has  the  value  z  ,  then: 

Tn  o 


P 

n 


(zTn  ~  zo>^ZTn  +  zo> 


The  work  described  here  introduces  the  following  changes  to  the 
previous  work  of  Harrington  [6,  7,  8,  9] : 

a)  Terminating  impedances  can  have  complex,  not  just  imaginary 
(reactive)  values. 

b)  Terminating  impedances  can  be  active  (have  gain)  so  that  the 
restriction  jp  |  <_  1  is  removed. 

c)  The  impedances  are  controlled  to  form  nulls  toward  sources 
of  interference. 

d)  A  closed  loop  adaptive  algoritnm  controls  the  z  for  pattern  control 
in  order  to  avoid  requiring  a  knowledge  of  the  array  geometry  and  to 
avoid  computing  tne  mutual  impedance  matrix  for  open  loop  pattern 
control. 

1.3  Scope 


The  SIC  program  has  been  undertaken  jointly  by  ZA1  and  Naval 
Research  Laboratory  (NRL) .  ZAI  has  completed  a  theoretical  analysis  of 
an  SIC  and  built  four  electronically  variable,  active  complex  terminations 
for  experimental  use.  NRL  has  built  an  experimental  parasitic  array.  Tests 
have  been  run  witn  this  array  using  passive  reactive  terminations  and  also 
using  the  active  complex  terminations. 

1 . 4  Glossary  of  Notation 

The  following  conventions  will  be  observed: 

1)  Upper  case  letters  will  represent  matrices  (e.g.,  M  or  P)  . 

2)  Lower  case  letters  with  an  underline  will  represent  vectors 
(e.g. ,  £  or  2^) . 

3)  Lower  case  letters  without  an  underline  will  represent  scalars 
(e.g. ,  yQ  or  pn) . 

4)  "0"  as  a  subscript  will  reference  the  main  antenna  (e.g.,  y^  or  m^) . 

5)  "+"  as  a  subscript  added  to  the  matrices  M,  P,  or  Q,  or  the 
vectors  i  or  y  will  denote  that  they  have  been  augmented  to 
include  the  corresponding  entries  of  the  ,main  antenna  in  addition 
to  those  of  the  parasitic  auxiliary  antennas.  Such  a  matrix 

or  vector  will  now  be  (N  +  1)  x  (N  +  1)  or  (N  +  1)  x  1  respectively 
(where  N  is  the  number  of  auxiliary  antennas) .  To  be  consistent 
with  Convention  4,  as  well  as  not  to  disturb  the  numbering  of  the 
auxiliary  antennas,  the  main  antenna  is  assigned  an  index  of  "0". 


mentally,  there  Is  little  information  available  on  how  the  array  para¬ 
meters,  mutual  impedances,  or  mutual  coupling  factors  of  compact  arrays 
vary  over  frequency.  Harrington  provides  computer  programs  that  can  be 
used  to  compute  these  parameters,  [12].  The  problem  is  generally  intrac¬ 
table  for  anything  but  numerical  methods. 

An  issue  that  arises  when  considering  a  parasitic  array  for 
nulling  is  whether  to  use  active  or  passive  terminations.  It  is  a  ques¬ 
tion  of  stability  versus  degrees  of  freedom  in  the  control  space. 

If  passive  terminations  are  used,  then  the  array  will  always 
be  stable.  The  mutual  coupling  between  the  antennas,  however,  is  neces¬ 
sarily  lossy.  That  is,  100%  of  the  power  is  not  coupled  to  other  antennas. 
It  is  this  loss  of  power,  through  re-radiation,  that  often  prevents  can¬ 
cellation  of  an  interference  with  only  one  auxiliary  antenna  or  N  inter¬ 
ferences  with  N  auxiliaries.  It  is  usually  necessary  to  use  two  or  more 
passively  terminated  auxiliary  antennas  to  cancel  a  single  interference. 
Thus  for  the  same  amount  of  nulling  capabilities  a  passively  terminated 
parasitic  array  will  require  more  elements  than  an  actively  terminated 
array . 

An  actively  terminated  parasitic  array  will  give  equal  or 
better  performance  with  fewer  antennas.  If,  however,  too  much  gain  is 
used  then  the  inherent  feedback  nature  of  the  parasitic  array  could 
cause  an  unstable  oscillation  to  start.  The  issue  of  stability  is  taken 
up  in  detail  in  Section  A. 2.  The  stability  analysis  determines  the  op¬ 
timal  mutual  coupling  matrix,  M+.  In  order  to  increase  the  usable  gain 
before  instability,  one  should:  1)  minimize  the  antenna  mismatch  or 
return  loss  or  equivalently  the  self  coupling  mnn  of  each  antenna,  2)  min¬ 
imize  the  mutual  coupling  between  auxiliary  antennas,  and  3)  maximize 
coupling  between  the  main  and  the  auxiliary  antennas.  The  optimal  mutual 
coupling  matrix,  M+,  is  of  the  form 


The  construction  of  such  an  ideal  mutual  coupling  matrix  can  be  approached 
by  two  methods,  antenna  element  and  array  design  and  by  connecting  various 
impedances  between  the  antennas  of  the  array. 


Another  aspect  of  parasitic  arrays  Is  the  Inherent  non-linear 
nature  of  any  control  algorithm  of  the  variable  terminations.  To  better 
understand  this,  consider  the  conventional  beamforming  array  and  the  para¬ 
sitic  array  illustrated  in  Figure  1.  Both  arrays  behave  linearly  in  tnat 
for  a  given  set  of  complex  weights  or  terminating  impedances  if  the  inci¬ 
dent  plane  wave  is  pealed  by  a  factor  c,  then  the  output  of  the  array  will 
be  scaled  by  c.  The  conventional  array  is  also  linear  in  a  control  sense, 
in  that  if  one  of  the  complex  weights  are  changed  by  an  amount  Aw,  then 
the  array  output  will  change  by  an  amount  proportional  to  Aw,  given  all 
other  weights  and  the  inputs  to  the  array  remain  constant.  It  is  in  this 
sense  that  the  parasitic  array  is  nonlinear.  If  one  of  the  reflectivities, 
p^,  were  to  change  by  an  amount  Ap,  the  array  output  in  general  will  not 
change  by  an  amount  proportional  to  Ap.  From  Equation  (2.2),  the  change 
in  the  output  (y^  of  y()  is  determined  by  the  matrix  inverse  (I  -  M+P+) 
Another  way  of  looking  at  this  is  by  the  gradient  of  the  array  output  with 
respect  to  the  control  variable,  either  the  weights  or  the  reflectivities; 
or  equivalently  be  considering  the  appropriate  individual  partial  deriv¬ 
atives.  In  the  conventional  array,  the  derivatives  of  the  array  output 
with  respect  to  the  weights,  3yg/8wn,  are  constant  and  depend  only  on  the 
incident  signals  and  array  geometry.  For  the  parasitic  array,  the  deriva¬ 
tives  of  the  output  with  respect  to  the  reflectivities  are  functions  of 
all  the  reflectivities  themselves  as  well  as  the  incident  signals  and  array 
geometry.  The  nonlinear  control  of  a  parasitic  array  is  discussed  in 
greater  detail  in  Chapter  5.  Note  a  mutual  coupling  matrix  as  in  Equation 
(2.9)  above  along  with  Pq  ■  0  would  allow  linear  control  of  the  array. 

In  this  case  (  I  -  M+P+  )  1  =  I  +  m+P+. 

Lastly,  in  this  chapter,  some  properties  of  a. transmitting  array 
are  discussed.  The  array  far  field  pattern  and  the  radiated  power  will  be 
calculated.  A  method  of  calculating  the  far  field  pattern  based  on  the  array 
impedance  parameters  is  described  by  Dinger  [4],  for  a  parasitic  array 
with  reactive  terminations.  The  results  are  easily  extended  for  general 
terminating  impedances.  Consider  the  impedance  model  of  an  array  illus¬ 
trated  in  Figure  5.  Let  the  voltage  sources,  v  ,  and  the  load  matrix, 

ZT,  be  the  N  +  1-port  Thevenin  equivalent  of  the  source  network  of  the 
array  during  transmission.  Note  that  the  polarity  of  the  voltage  sources 
is  the  opposite  of  the  usual  convention.  To  calculate  the  far  field,  one 
must  first  find  the  port  currents,  i^,  of  the  array. 
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■  -(2i  +  Vl^c 


Under  the  assumption  that  the  current  on  each  antenna  is  determined  only 
by  its  port  current,  then  one  can  use  superposition  to  add  up  the  con¬ 
tribution  of  each  antenna  in  the  far  field. 

Let  us  assume  that  all  the  antennas  are  identical  monopoles  and 
orientated  in  the  z-direction.  Then  the  far  field  pattern  as  a  function 
of  the  azimuthal  angle,  <p,  is  given  below,  [7]. 

EZ  -  -  j4-'J-kr  t  *  «Jk  C0S  *  +  Si“  W  (2.10) 

n-0 


where  r  is  the  distance  from  the  array,  n  is  the  intrinsic  impedance  of 

free  space,  and  and  y^  are  the  coordinates  of  the  nth  antenna  relative 

to  the  array  phase  center.  Note  it  is  assumed  lx  I,  |y  I  «  r. 

n  n 

A  comment  on  the  assumption  necessary  for  superposition.  It  is 
not  being  assumed  there  is  no  mutual  coupling.  The  calculation  of  the 
port  current  takes  that  into  account.  Superposition  holds  for  the  cur¬ 
rents  on  the  antennas.  But  the  calculations  of  Equation  (2.6)  give  us 
only  the  port  currents.  Thus  to  use  superposition  we  must  assume  that 
the  currents  on  an  antenna  are  due  to  its  port  current  alone  and  that  this 
antenna  current  distribution  will  be  the  same  as  that  on  an  isolated  an¬ 
tenna.  To  illustrate  further,  consider  an  array  of  a  particular  config¬ 
uration  with  a  set  of  source  voltages.  Assume  that  this  establishes, 

inclusive  of  all  mutual  coupling  effects,  a  current  i  at  the  port  of  the 

n 

nth  antenna.  Now  let  us  move  the  other  antennas  around  and  change  the 

source  voltages  but  always  under  the  condition  that  the  port  current  of 

the  nth  antenna  is  maintained  at  i  .  The  assumption  then  asserts 

n 

that  currents  on  the  nth  antenna,  which  give  rise  to  its  radiated  field, 
will  not  significantly  change.  This  assumption  is  good  for  antennas  that 
are  small  compared  to  a  wavelength  (<  X/4)  even  in  a  compact  array.  The 
port  current  alone  then  will  determine  all  antenna  currents.  If  this 
assumption  cannot  be  made  then  the  far  field  pattern  can  be  calculated 
using  Harrington’s  computer  programs  [ 12 ] . 

If  one  wishes  to  use  transmission  line  parameters.  Equation 
(2.10)  above  is  still  used  but  the  antenna  port  currents  are  calculated 
differently.  In  Appendix  C.2,  an  expression  for  the  antenna  port  cur- 
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rents  in  terms  of  the  transmission  line  parameters  is  derived.  In  the 
notation  of  Appendix  C.2,  _i(A)  is  the  vector  of  the  antenna  port  cur¬ 
rents. 


1(A)  «  \  Zo1  (D  -  D_1M+)  (Q~V  (I  -  P+)  v 


-DC 


(C32) 


where  Z  is  the  diagonal  matrix  of  the  transmission  line  characteristic 
o 

impedances,  D  is  the  diagonal  matrix  of  delays  through  these  transmission 

lines  (see  Equation  3.17)  and  v  is  the  vector  of  voltage  sources  in  the 

— oc 

terminations . 

Lastly,  we  calculate  the  power  radiated  by  a  transmitting  array. 
With  an  impedance  model  the  average  power  is  given  by  Equation  (2.11). 

P  -  1/2  Re  [i*T 
-  1/2  Re  [i*T  ZAij 


1/2  i*T  Re  [ZA]i 


(2.11) 


v-1 


where  i  *=  —  (Z  +  Zj)  v  .  The  last  step  is  true  because  Z,  is  symmetric. 
—  AT  — oc  A 

For  the  transmission  line  model  the  power  radiated  (assuming  loss¬ 
less  transmission  line)  is 

P  -  1/8  Re  [/*  (I  -  S*TSX)  v^]  (2.12) 

where  ST  is  given  in  Equation  (3.13)  or  (3.16). 
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3.0  RELATIONSHIPS  BETWEEN  TRANSMISSION  LINE  AND  IMPEDANCE  MODELS  OF  A 
PARASITIC  ARRAY 

In  Chapter  3,  we  shall  examine  the  relationship  between  the 

transmission  line  model  and  the  impedance  model  of  a  parasitic  array.  This 

relationship  will  be  illustrated  by  analyzing  various  aspects  of  receiving 

and  transmitting  arrays  with  each  model.  In  particular,  it  is  shown  how 

to  calculate  the  parameters  of  the  transmission  line  model  from  the  impedance 

parameters.  The  means  of  relating  the  two  models  will  be  scattering  parameters. 

The  derivation  of  the  matrix  of  scattering  parameters,  S,  is  given  in  Appendix 

B.  The  scattering  matrix  is  a  generalization  of  the  reflection  coefficient, 

T,  of  a  single  port  system.  For  an  N-port  network,  with  no  sources,  the 

i,jth  scattering  parameter,  s^,  is  defined  as  the  ratio  of  the  voltage 

traveling  wave  coming  out,  "reflected",  from  the  ith  port  due  to  the 

wave  incident  on  the  jth  port  with  no  other  incident  waves. 

The  matrix  of  these  parameters,  S  =  [s^],  can  be  calculated 

from  Z  and  Z  .  Here,  Z  is  the  open  circuit  mutual  impedance  matrix 
o 

of  the  N-port  network,  and  Zq  is  the  diagonal  matrix  of  the  characteristic 
impedances  of  the  media  connected  to  the  network  ports.  We  repeat  here 
equation  (B13) . 


S  *=  (Z  -Z  )  (Z  +  Z  ) 
o  o 


(B13) 


One  comment  be'ore  continuing,  in  the  previous  chapter  it  was 
assumed  there  were  no  transmission  lines  present  with  the  impedance  model. 
Hereafter,  transmission  line  will  be  present  and  their  affects  explicitly 

calculated . 

We  begin  with  a  treatment  of  the  single  element  case. 

3.1  Single  Antenna  Transmitting 


The  impedance  model  of  a  single  antenna  transmitting  is  given 
in  Figure  6a.  The  difference  between  the  transmitting  and  receiving 
versions  of  the  impedance  model  is  the  placement  of  the  voltage  source 
(compare  to  Figure  5) . 

For  a  single  antenna,  the  scattering  matrix  will  be  a  scalar, 
s,  which  is  simply  the  reflection  coefficient,  T,  at  the  same  point. 

We  shall  retain  the  notation,  s,  so  that  the  generalization  to  an  N-port 
network  will  be  seen  more  readily. 


6  a)  Impedance  Model  (with  transmission  line) 


6  b)  Voltage  Traveling  Wave  Relationships  within  the 
Transmission  Lines 


Figure  6  Impedance  Model  of  Single  Element  Array  Transmitting 


In  Figure  6b,  the  Important  relationships  among  the  various 
voltage  traveling  waves  are  illustrated,  s  •  (z  -  z  )/(z  +.  z  )  is 
the  scattering  parameter  or  reflection  coeflcient  at  the  antenna  port, 
p  m  (z^,  -  ZQ)/(Zy  +  zq)  is  the  reflection  coefficient  of  a  voltage  wave 
traveling  toward  the  source,  incident  on  the  termination.  The  trans¬ 
mission  coefficient  of  a  voltage  wave  incident  on  a  junction  is  the 
ratio  of  the  voltage  transmitted  to  the  voltage  incident.  Ignoring  any 
further  reflections,  continuity  of  the  total  voltage  means  the  trans¬ 
mission  coefficient  is  1  +  T.  Thus,  for  a  voltage  wave  toward  the  source 
and  incident  on  the  terminations,  the  transmission  coefficient  is  1  +  p. 

For  a  voltage  wave  from  the  source  and  incident  on  the  termination- 

transmission  line  interface,  the  reflection  coefficient  is  (z  -  z„,)/(z  +  zj’  * 

o  T  o  T 

ignoring  subsequent  reflections  off  the  antenna.  Thus,  for  a  similar  voltage 
wave  the  transmission  coefficient  is  I  -  p.  The  scattering  parameter,  s^,,  for 
a  voltage  wave  traveling  toward  the  antenna  and  incident  on  the  termination- 
transmission  line  interface  can  be  calculated  from  the  results  of  Appendix  C.l. 


sT  =  (zA  -  zT)  /  (zA  +  ZT)  from  Equation  (C4) 


[o  +  sA>  /  (1  -  sA)]  *0  -  *T 

7  ~  ~  T  from  Rule  2  of  (3.1) 

+  1  “  s^) J  z0  +  zj  Appendix  C.l 


where  zA  and  sa  are  respectively  the  open  circuit  impedance  and  scattering 
parameter  (reflection  coefficient)  looking  at  the  antenna  within  the 
antenna  transmission  line  just  before  the  termination.  We  can  further 
manipulate  Equation  (3.1)  to  obtain: 


3-3 


mutual  coupling  factor,  m,  and  noting  that  when  z T  **  zQ  the  termination- 

transmission  line  interface  becomes  invisible,  we  see  that  m  =  s^  (with 

z  =  z  ) .  Indeed  this  can  serve  as  an  alternate  definition  of  m. 

X  o 


m 


sT  ,  when  z T  =  zq  or  equivalently  p  *  0. 


d  K  “  0/<**  +  *  ) 


(3.3) 


(3.4) 


—  1 

where  d  *  e  ,  3  is  the  propagation  velocity  of  the  transmission,  line, 

and  i.  is  the  length  of  the  line.  Note  that  s.  is  independent  of  z_,. 

A  X 

Thus,  Equations  (3.3)  and  (3.4)  allow  us  to  calculate  the  mutual  coupling, 

m.,  from  the  impedance  parameters  z,  and  z  . 

A  o 
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We  shall  now  obtain  Equation  (3.2)  via  the  transmission  line 
model.  Consider  a  voltage  wave  traveling  toward  the  antenna  and  incident 
on  the  termination-transmission  line  interface.  Let  us  call  such  a  vol¬ 
tage  wave  (see  Figure  7).  The  voltage  that  gets  transmitted  into  the 


transmission  line,  ignoring  subsequent  reflections,  is  then  (1  -  p)  v 


(+)* 


This  voltage  wave  toward  the  antennas  becomes  a  voltage  wave  toward  the 

source  Incident  on  the  termination  of  m  (1  -  p)  v^  ,  again  ignoring 

further  reflections.  The  resultant  (after  all  reflections)  voltage  wave 

Incident  on  the  termination  from  the  antenna,  y,  can  be  computed  from 

Equation  (2.2)  with  y  =  m  (1  -  p)  v,,.  . 

s  (+; 


y  =  (1  -  mp) 


-1 


-1 


=  (1  -  mp)  m  (1  -p)  v 


(+) 


(3.5) 


The  total  downward  voltage  incident  on  the  termination  transmits 
a  voltage  (1  +  p)y  into  the  termination,  traveling  toward  the  source. 

There  is  another  component  of  this  voltage,  -pv^  ,  due  to  the  initial 


reflection  of  the  original  incident  voltage  v 


(+) 


Thus,  we  can  write: 


V(+) 


(1  +  p)y  -  pv 


(+) 


-1. 


[(1  +  p)(l  -  mp)  (1  -  p)  -  p ] v 


(+) 


ana 


sT  =  (1  +  p)  ( 1  -  mp)  ^(1  -  p)  -  p 


(3.6) 


(1  -  P  )m 


1  -  mp 


m  -  p 


-  P 


1  -  mp 


(3.7) 


3-5 


----- 


(+) 


(l-P)v 


(+) 


"  ^ 
Antennna 
Impedance 


Figure  7  Transmission  Line  Model  of  Single 
Element  Array  Transmitting 
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But  we  know  that  m  —  s.»  therefore  Equation  (3-7)  is  identical  to 

A 

Equation  (3.2). 

In  this  Section,  we  have  demonstrated  that  m  *  ■  B^Cp) |p«g 

and  have  calculated  s^  through  both  the  impedance  and  transmission  line 
models . 


3*2  Single  Antenna  Receiving 

The  impedance  model  of  a  single  antenna  during  reception  is 

given  in  Figure  8a.  Figure  8b  illustrates  the  accompanying  traveling 

wave  parameters.  Let  us  find  the  voltage  across  the  termination,  v^. 

Let  v  be  the  voltage  across  the  antenna  port.  Then  we  can  find  v 
A  A 

from  the  current  out  of  the  port,  i  . 

A 


V 


A 


A 

-2 


T^ 


(3.8) 


But,  iA  is  given  by 


■  -<ST  +  *A>~\C 


(3.9) 


therefore,  by  combining  Equations  (3.8)  and  (3.9),  we  obtain 


v  =  z  (z  +  z  )  *v 
A  T  T  A  oc 


(3.10) 


From  Rule  2  of  Appendix  C.l,  we  know  that 


2T  =  (1  +  p)(l  -  p)  ^ 

1  o 


where  p  is  the  reflection  coefficient  seen  looking  into  the  transmission 
line  at  the  anter.na.  Substituting  this  into  Equation  (3.10),  we  obtain 


Transmission  Line 
Impedance  =  z 

o 


O 


z 

A 


8a)  Impedance  Model 


8b)  Voltage  Traveling  Wave  Relationships  witnin  the 
Transmission  Lines 


Figure  8  Impedance  Model  of  Single  Element  Array  Receiving 


vA  =  (1  +  d)(l  -  p)"1zo  [(1  +  £)(1  -  P)_1zo  +  Za]_1voc 

=  (i  +  [d  +  2)z  +  (i  -  *US\C 

o  o  A 

■  +  «zo[(zA+Zo)  -  »(zA  -  Zo)flvoc 

-  (i  +  5)zo(Za  +  ./‘u  -  «*A  -  *0uzA  +  *„>'1riv 


But  s  ■=  (z  -  2  )(2,  +  z  )  1  and  z  (z,  +  z  )  1  -  1/2(1  -  si,  therefore, 
A  A  o  A  o  oAo  A 

vA  =  1/2(1  +  p)(l  -  sA)(l  -  PsA)_1voc  (3.11) 

To  relate  to  vT>  we  must  consider  the  traveling  waves  within 
the  transmission  line.  Let  y  be  the  voltage  traveling  toward  the  term¬ 
ination,  measured  at  the  antenna  end  of  the  transmission  line.  Since 
there  is  no  source  in  the  termination,  the  total  voltage  at  the  antenna 
end  of  the  transmission  line  is  y  plus  the  corresponding  reflected 
voltage  py.  Therefore, 


v  “  (1  +  P)y  and 

A 

y  =  (1  +  P)_1vA 

If  y  is  the  voltage  traveling  toward  the  termination,  incident  on  the 
termination,  then  y  =  dy,  where  d  =  e  is  the  delay  associated  with 

the  transmission  line  of  length  1.  Finally,  the  total  voltage  at  the 
termination,  v^,  is  the  sum  of  incident  and  reflected  voltages  at  the 
termination. 
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vT  =  (1  +  P)y 

*  (1  +  P)dy 

=  (1  +  p)d(l  +  ^)_1vA 


=  1/2(1  +  p)d(l  -  s  )  (1  -  psA)_1v 

A  A  oc 


*  1/2(1  +  p)(l  -  s.p)  ^d ( 1  -  s . ) v 

A  A  oc 


(3.12) 


where  p  =  d  p  and  s. 


If  we  compare  v  to  the  definition  of  y  ,  we  see  that 
i  s 

ys  =  VT  With  ZT  =  Z0  °r  etJuivalently  p  =  0.  Similarly,  as  in  the  case 
of  m,  this  could  serve  as  an  alternate  definition  of  y  . 


ys  =  VT 


p=0 


-  1/2d<‘  -  *A>  Voc 


(3.13) 


Let  us  examine  Equation  (3.13).  From  Equation  (C6)  of  Appendix 

C.l  we  know  that  the  traveling  voltage  from  the  antenna,  incident  on  the 

antenna-transmission  line  interface,  is  simply  (1/2)vqc.  The  reflection 

coefficient  seen  from  the  antenna  is  -s.  and  the  transmission  coefficient 

A 

is  then  1  -  s^.  Thus,  the  traveling  voltage  wave  that  gets  transmitted 
through  this  interface  is  (1/2) (1  -  s.)v  .  y  is  simply  this  voltage 

A  OC  S 

multiplied  by  d  to  account  for  the  delay  in  the  length  of  the  transmission 
line. 

Equation  (3.12)  can  also  be  quickly  derived  through  the  trans¬ 
mission  line  model.  Let  the  argument  in  the  preceding  paragraph  establish 
Equation  (3.13).  Then  the  total  voltage  in  the  transmission  line  inci¬ 
dent  on  the  termination,  y,  is  given  by  Equation  (2.2)  as  applied  to  a 
single  antenna  case. 

y  =  1/2(1  -  mp)_1d(l  -  s  ) v 

A  oc 


3-10 


The  vector  vi  is  physically  dimensionless.  While  its  physical  meaning  is 
not  obvious,  u  does  contain  all  the  information  of  the  incident  signals 
necessary  to  form  the  nulls.  From  Equation  (4.7)  we  can  write 

P(Mu  -  n>0  )  =  u  ,  '4.8 

Since  P  is  a  diagonal  matrix.  Equation  (4.8)  represents  N  equations  of 
the  form 


n  ”  [Mu]  -  m  (4.9) 

—  n  on 

where  [*]  represents  the  nth  entry  of  the  vector 
n 

argument. 

There  are  possible  sets  of  initial  signal  vectors,  y  ,  which 

‘■sn 

will  result  in  a  u  that  causes  the  denominator  of  Equation  (4.9)  to  be 
zero.  This  means  the  array,  for  that  combination  of  the  N  incoming 
signals,  loses  a  degree  of  freedom  in  the  control  space  and  will  only 
be  able  to  null  N-l  signals.  This,  however,  is  not  a  particularly 
scandalous  result.  Similar  behavior  is  observed  in  a  conventional  adap¬ 
tive  array  with  an  unweighted  main  antenna.  Consider  the  3  element  array 
in  Figure  13.  Since  the  arrival  directions  of  the  signals  are  symmetric 
about  an  axis  through  the  two  auxiliary  antennas,  the  relative  phase 
between  the  two  signals  will  be  the  same  at  both  antennas.  If  the  two 
antennas  have  identical  azimuthal  gain  patterns,  then  the  signal  at  one 
antenna  will  be  linearly  dependent  on  the  other,  i.e. ,  the  signals  will 
be  equal  except  for  a  complex  scalar.  In  this  case  only  one  of  the  sig¬ 
nals  can  be  nulled. 

The  concept  of  degrees  of  freedom  is  a  useful  concept  in 
adaptive  arrays.  In  a  conventional  array  with  an  unweighted  main  antenna, 
the  number  of  degrees  of  freedom  is  generally  the  number  of  weighted 
auxiliaries.  Similarly  in  a  parasitic  array,  with  a  constant  termination 
on  the  main  antenna,  the  number  of  degrees  of  freedom  is  generally  equal 
to  the  number  of  active  complex  terminations .  The  number  of  signals  that 
can  be  canceled  is  equal  to  the  number  of  degrees  of  freedom. 

A  degree  of  freedom  can  be  lost  in  several  ways.  In  the 
example  above  in  Figure  13,  a  degree  of  freedom  was  lost  because  the  sig¬ 
nals  were  linearly  dependent  at  the  weighted  antennas  but  the  signal  at 
the  main  antenna  was  independent  of  the  weighted  antenna  signal.  Also 
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Now  the  equations  of  (4.1)  can  be  rewritten  as 

T  T  T  T 

^0  =  ^so  +  moopoio  +  2.0  PY 


Y  =  YS  +  EoPoio  +  MPy 


(4.2) 


If  the  reflectivities  needed  for  a  null  are  used  then  =  0,  and 
Equations  (4.2)  become 


-  T  T 

°  =is  0  +H0  Py 


(4.3a) 


Y  =  Y  +  MPY 
S 


(4.3b) 


Let  us  solve  Equation  (4.3b)  for  Y  (compare  to  Equation  2.5b) 


Y  =  (I  -  MP)  Y 


(4.4) 


Substitute  Equation  (4.4)  into  Equation  (4.3a). 


0  =  is0T  +  HL0Tp  (I  “  MP)~lys 


(4.5) 


If  Y  is  not  invertible,  that  means  the  set  of  vectors  of 
s 

initial  auxiliary  voltage  waves  is  linearly  dependent.  The  system  is 
not  fully  constrained.  Two  signals  may  be  arriving  form  the  same  direc¬ 
tion.  In  the  process  of  nulling  N-l  signals  the  Nth  will  be  nulled  for 


"free". 


Assuming  that  all  the  signal  vectors  are  independent  let  us  post 


multiply  both  sides  of  Equation  (4.5)  by  Y 

s 

T  -1  T  -1 

0  =  y  Y  +  m  -  P(I-MP) 

—  sO  s  —0 


Now  let  us  post  multiply  both  sides  of  Equation  (4.6)  by  (I-MP) . 


T  -1  T 

0  =  y„Y  (I  -  MP)  +  m  P 
—  sO  s  —  0 


t  —  1  T  T  -1 

=  y  „  Y  +  (in  -  y  Y  M)P 
—  sO  s  —0  —  sO  s 


T  ,  T  T  „ 
=  u  +  (m  -  u  M)  P 


(4.7) 


T  T  -1  -1 

where  u  =  y  „  Y  oru=(Y  )  y  „ 

-  —  sO  s  -  s  —  sO 
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4.0  NULL  FORMING  CAPABILITY  OF  AN  ACTIVELY  TERMINATED  PARASITIC  ARRAY 


4.1  Forming  N  Nulls  on  N  Simultaneous  Signals 

In  this  section  on  forming  multiple  nulls,  it  has  been  assumed 
that  all  the  signals  are  at  the  same  frequency.  This  was  required  so 
that  the  array  coupling  matrix,  M+,  would  be  the  same  for  all  the  signals 
desired  to  be  canceled.  In  Chapter  5  when  adaptive  control  is  discussed 
it  will  also  be  assumed  that  these  signals  are  uncorrelated.  This  con¬ 
straint  is  applied  to  keep  the  analysis  tractable.  Cancellation  of 
signals  of  different  frequencies  is  possible  and  is  discussed  briefly 
later,  but  a  detailed  analysis  was  not  possible  under  this  contract. 

Given  a  parasitic  array  of  a  main  antenna  and  N  auxiliaries, 
let  us  find  the  reflectivities,  p^,  necessary  to  null  N  signals.  Let  each 
incident  signal  be  expressed  by  the  voltage  waves  that  the  signal  pro¬ 
duces  in  each  transmission  line  incident  on  the  termination,  in  the  absence 
of  the  other  signals  and  with  all  terminations  «  zQ  (the  transmission  line 
characteristic  impedance) .  Let  the  vector  of  the  voltage  waves  at  the  N 
auxiliary  antennas  from  the  nth  signal  be  v  q  and  the  main  antenna  vol¬ 
tage  y  .  Similarly,  let  be  the  vector  of  final  voltage  waves  inci¬ 
dent  on  the  auxiliary  termination  and  yQn  the  final  voltage  wave  incident 
at  the  main  antenna  termination  for  the  current  set  of  reflectivities 
due  to  the  ntn  signal.  Tnese  quantities  can  be  related  by  the  following 
equations  (see  Equation  2.3). 


(4.1) 

signal  source  respectively, 
following: 

Y  =  £z  2  I  •  •  •  I  Z  N  3  -  : 

and  YS  =  [ysl|...|ysN]  . 


yoi  =  ysoi  +  moopoyoi  +  2oTp*i 
*1  =*sl  +^0P0y01  +MP*1 

T 

y0N  ‘  ysON  +  "OOVOH  *  =  0  P*» 

The  i,j  indices  of  each  y. .  denote  the  antenna  and 
To  express  these  more  compactly,  let  us  define  the 


*'•  -  ••  •.*  •_  '.*  % 
•  .  ■  »*  »  ’  •  ■  -  '  •  „•  -  .•  »“ 

,  w.v.  * ■■  in, 


iVi.VArt.Vi 
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the  antenna  sources  v  ,  incident  on  the  transmission  line  antenna  port 

—  oc 

interfaces,  will  be  h  v  .  The  voltages  transmitted  into  the  trans- 

—  oc 

mission  lines  will  be  ^  (I  -  S„)v  .  y  is  simply  the  vector  of 

A  —  oc  s+ 

transmitted  waves  as  they  appear  at  the  other  end  of  the  transmission 
lines,  i.e.,  phase-delayed  by  D. 
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Now  knowing  SA  (the  scattering  matrix  of  the  antenna)  =  (z^  -  zqI) ( z A  +  zqI) 
and  using  identities  A. 3  and  A. 6  reduces  Equation  (3.26)  further. 

-  MU-  Pt)(I  -  Sjf.)"1!!  -  SA)voc  (3.27) 

The  voltages  of  interest  however  are  v^. 


Zr  -  (i  +  p+>2. 

=  (i  +  p+)d£ 


(3.28) 


2.  is  the  vector  of  voltage  traveling  waves  entering  the  transmission  lines 

at  the  antenna  ends  and  flowing  toward  the  terminations;  therefore, 

v  =  (I  +  P.) v,  and  y  =  (I  +  P  )  ^v  .  We  now  substitute  this  expression 
— A  +  +  — ft 

for  ^  into  Equation  (3-28) . 

vT  -  (I  +  P+)  D  (I  -  P+)_1vA 

=  Ml  +  P  )  D  (I  -  S  P  ) -1  (1  -S  ) v 

+  A  +  A  — oc 

=  *a(I  +  P  )  D  (I  -  S  DP  D)_1(I  -  S  )v 

+  A  +  A  — oc 

«  Ml  +  P.)  (I  -  DS  DP  )  _1D(I  -  S  )  v 

+  A  +  A  — OC 

-  Ml  +  P.)  (I  -  S  P  }  _1D(I  S  )  v 

+  A  +  A  —  oc 


Ml  +  P+)  q"  D(I  -  SA)V 


A  —  OC 


(3.29) 


Equation  (3.29)  shows  the  relationship  of  the  voltage  across  the  termina¬ 
tions  v,j,  to  the  open  circuit  voltages  using  transmission  line  parameters. 
By  comparing  Equations  (3.28)  and  (3.29),  one  can  find  an  expression  for 
=  Q+X  ter®s  °f  v  oc  . 


*s+ 


*jD  (I  -  S. )  v 

A  —  oc 


(3.30) 


Let  us  analyze  this  equation  a  little.  In  a  manner  similar  to  the 
derivation  of  Equation  (C18)  it  can  be  shown  that  the  voltage  waves  from 
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Figure  12  Impedance  Model  of  Multiple  Element  Array  Receiving 


The  impedance  model  of  a  multiple  antenna  array  receiving  is 

given  in  Figure  12.  As  in  Section  3.3  the  indices  for  N  +  1  antenna 

elements  shall  be  0/  . ..,  N,  where  0  is  the  main  antenna.  Figure  12 

illustrates  the  various  transmission  line  parameters  for  a  receiving  array. 

The  usual  quantities  of  interest  for  a  receiving  array  are  the  voltages 

across  the  terminations  v„. 

— T 

Within  the  transmission  line  model,  v^  is  the  vector  of  voltages 
transmitted  to  the  termination;  therefore, 

Zj.  -  <i  +  VV  (3*22) 

Equation  (2.2)  tells  us  ^  =  Q+  ^~y^ |  .  This  can  be  substituted  into 
Equation  (3.22)  to  obtain  v^  in  terms  of  the  transmission  line  parameters. 

v„.  -  (I  +  P+) Q+’1XS+  •  (3.23) 


For  the  impedance  model  the  calculation  is  more  involved,  unless 
the  transmission  lines  are  of  zero  length,  in  which  case  Equation  (2.8) 
may  be  employed.  With  transmission  lines  one  must  use  the  equations  of 
Appendix  C.2.  First  one  must  find  v^,  the  vector  of  port  voltage  of  the 
array . 


v  =  2  (2 
-A  T  T 


+  ZJ  1v 
A  — OC 


(3.24) 


where  ZT  is  the  diagonal  matrix  of  the  terminating  impedances  as  seen 

from  the  antenna  end  of  the  transmission  lines.  From  Rule  2  of  Appendix  C.2, 

we  find 

ZT  =  (I  +  P+)(I  -  P+) (3.25) 

A 

where  P+  is  the  diagonal  matrix  of  reflection  coefficients  of  the  termi¬ 
nations  translated  down  the  transmission  lines  to  the  antennas  and  z 

o 

(a  scalar)  is  the  characteristic  impedance  of  the  transmission  lines. 
Substituting  Equation  (3.25)  into  (3.24)  and  manipulating  the  equations 
somewhat  yields 


v*  =  (I  +  P+)  tl  + 


<V  +  - 


Z- ) P  ]z  (z  I  +  Z  ) 
A  +  o  o  A 


V 

— OC 


(3.26) 
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At  this  point,  let  me  restate  the  conventions  on  the  use  of  "+" 

as  a  subscript,  as  put  forward  in  the  Glossary  of  Section  1.4.  With  the 

variables:  y ( ,  M+,  P+,  or  Qf,  the  subscript  signifies  that  the 

variable  describes  the  entire  array  of  main  and  auxiliary  antennas.  When 

"+"  is  used  with  parentheses  with  the  variables:  v^  or  v ^  ,  it  signifies 

that  the  variable  represents  a  traveling  wave  or  waves,  and  in  particular 

voltage  waves  traveling  toward  tne  antennas  (see  Appendix  C) . 

The  individual  total  voltage  wave,  y  ,  incident  on  each  termination 

n 

interface  transmits  a  voltage  wave  of  (1  +  p^Jy^  through  the  interface  and 
into  the  termination.  This  is  expressed  for  all  the  voltage  waves  in  vector 
form  by  (I  +  P| )y( . 

The  total  voltage  wave  within  the  nth  termination  traveling  away 

from  the  termination- transmission  line  interface  is  the  sum  of  (1  +  p  )y 

n  n 

and  an  initial  reflection  from  v  ,  .  incident  upon  the  termination-transmission 

n(+) 


line  interface.  This  initial  reflection  is  -p  v  .  , 

n  n(+) 


sum 


is  (I  +  P+)^  - 


P  v 


(+) 


In  vector  form,  this 
However,  we  also  know  that  this  total  voltage 


wave  is  simply  (by  definition)  S^v  . 


£A+) 


-  11  +  V*.  '  p^c+, 


(3.20) 


■[ 


-i. 


Thus  we  have  found  an  expression  for  S^,. 

ST  =  (I  +  P+)(l  -  M+P+)_1M+(I  -  P+)  -  P+  (3.21) 

=  (i  +  p+)Q+"1m+u  ■  V  "  p+ 

It  an  be  shown  that  this  form  for  ST  is  equivalent  to  Equation  (3.16) , 
noting  that  M+  =  §A-  This  is  done  in  Appendix  D. 

3.4  Multiple  Element  Array  Receiving 

The  material  presented  in  Section  3.2  is  generalized  in  this 
section  for  the  multiple  antenna  array.  As  in  Section  3.3  voltage  and 
currents  will  now  become  vectors  of  voltages  and  currents.  The  reflection 
and  transmission  coefficients  will  become  matrices. 


Propagation  Velocity 


Transmission  Line  Model  of  Multiple  Element  Array  Transmitting 


(Jj 

V 


Note  that  §A  is  independent  of  the  terminating  empedances,  z^'s.  Thus 
M+  13  and  be  used  to  calculate  the  mutual  coupling  matrix  from  the 
antenna  impedance  matrix,  Z  ,  the  characteristic  impedance  z  ,  and  the 

A  O 

transmission  line  length  £  's. 

n 

In  addition,  one  can  solve  Equation  (3.17)  for  in  terms  of 

A 

S.  *=  M, .  If  desired  one  can  calculate  Z„  from  M. . 

ZA  **  ZQ  [I-D~1M+D_1]  [I  +  (3.18) 

We  shall  now  obtain  Equation  (3.16)  from  the  transmission  line 

model .  Consider  a  voltage  wave ,  vn  ^  +j  ^  from  each  source ,  traveling  toward 

the  antenna  and  incident  on  the  termination-transmission  line  interfaces. 

Let  all  these  voltage  waves  be  combined  into  the  vector  v ^ .  At  the  nth 

interface,  the  transmitted  wave  (ignoring  further  reflections)  is  then 

(1  -  Pr)  vn(+)*  Considering  all  N  +  1  interfaces,  in  vector  notation  this 

becomes  (I  -  P  )  v  where  P  is  the  diagonal  matrix  of  the  N  +  1  reflec- 
+  \ + j  s  + 

tion  coefficients  (see  glossary  in  Section  1.4) .  These  voltage  waves  are 
illustrated  in  Figure  11.  These  voltage  waves  toward  the  antennas  are 
coupled  (or  scattered  or  reflected)  through  the  antennas  to  produce  reverse 
traveling  waves  toward  the  terminations.  The  reverse  voltage  waves  inci¬ 
dent  upon  the  terminations  from  a  single  round  trip,  in  vector  notation 
are  M+  (I  -  P+)  v  ^ .  Let  us  consider  these  initial  reverse  traveling 
waves  as  equivalent  initial  waves  produced  by  an  imaginary  incident  plane 
wave.  Let  the  vector  of  these  voltage  waves  be  ^+.  Fro®  Equation  (2.2), 
the  final  (after  all  reflections)  voltage  waves  traveling  toward  the 
terminations  and  incident  on  the  termination- transmission  line  interface 
will  be  where 

-  <1  -W"1  X,.  <2-2> 

=  (I  -  M+P+)_1  M+  (I  -  P+)  v(+)  (3.19) 
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the  cables, 
the  scattering 

(3.14) 

Thus  Equation  (3.13)  can  be  rewritten  as 

ST  '  [[I  -  SAr‘  11  +  K]  \  -  h]  [!I  '  §Arl  [I  +  §A]  +  Zl]'‘  <3-15) 

Equation  (3.15)  is  the  generalized  form  of  Equation  (3.1).  Let  us  manipulate 
further  this  form  for  S^.  From  Identity  A. 3  of  Appendix  A,  we  see  that 

[i  -  sA]_1  [i  +  sA]  -  [i  +  saHi  -  sAr1. 

ST  *  t11  +  SAJ[I  -  §Arl  *0  -  *l]  [[I  +  §A1IJ  -  §Ar‘  ’o  +  Zl]  _l 

-  [[I  +  SA]  zo  -  ZT  [I  -  Sj)  [I  -  SArl  [I  -  SA]  [[I  +  SA]zo  +  ZT[I  -  SA 

‘  K  +  *oll§A  -  tZT  '  ‘a1’]  [[ZT  +  “o11  *  lZT  *  V*A]-1 

■  [ZT  +  “o11  (SA  '  [ZT  +  Zoirl  [ZI  -  ‘o11] 

.  [i  -  [ZT  +  *0I]_1EZt  -  r0USA]  ’Szj  +  x0U_1 

By  applying  Identity  A. 3  again,  [ Z T  +  zqI]  1 [ZT  -  zoI]  «*  [ZT  -  zoII [ZT  +  zqI]  1 
■  P+,  and 

ST  -  (ZT  +  zQIJ (SA  -  P+] [I  -  P+SA]_1[ZT  +  zoI]_1  (3.16) 


Z.»  the  impedance  matrix  seen  at  the  termination  end  of 

^  Ak  A 

can  be  calculated  from  S.  and  Rule  3  of  Appendix  C.2,  where  S.  is 

A  .  A 

matrix  seen  from  the  termination  ends  just  inside  the  cables. 


ti  -  sAi 


-1 


U  +  SA]  zo 


Here  Equation  (3.16)  is  the  generalized  form  of  Equation  (3.2).  As  in  the 
single  antenna  case  let  us  now  set  all  the  terminating  impedances  to  zfl, 
i.e.,  ZT  ■  zqI,  or  equivalently  P+  **  0. 

=  2ZoI[gA  ■  °HI  -  °-Sa]-1[2  zoI]-1 


But  ST  when  Z T  ■  zqI  is  the  definition  of  M+,  the  matrix  of  mutual  coupling 
factors.  Therefore  it  is  always  the  case  that  M+  “  §A> 

Ve  can  also  calculate  SA  from  Rule  3  of  Appendix  C.2. 


P,=0 


S  -  D  [Z.  -  z  I] [Z.  +  z  I]-1D 


(3.17) 


Figure  9  Impedance  Model  of  Multiple  Element  Array  Transmitting 


The  voltage  across  the  termination,  v^,  by  continuity  is  the  total 
voltage  at  the  end  of  the  transmission  line.  The  total  voltage  is 
the  sum  of  the  incident  and  reflected  traveling  waves.  Thus,  we  can 
write 

vT  -  (1  +  p)y 

-  (1  +  p) (1  -  mp)-1yg 

“  4  (1  +  p)(l  -  mp)  1d(l  -  s  )v 
l  A  OC 

If  we  remember  that  m  *  s^,  then  we  see  that  the  above  equation  is  indeed 
Equation  (3.12). 

In  this  section  we  have  derived  the  relationships  that  exist 

between  the  transmission  line  model  and  the  impedance  model  of  a  single 

antenna  while  transmitting.  The  voltage  across  the  termination  was 

derived  through  both  models,  and  we  obtained  alternate  definitions  for 

y  ,  the  incident  voltage  on  a  matched  termination  due  to  a  far  field 
8 

source,  and  for  m,  the  mutual  coupling  coefficient. 

3.3  Multiple  Element  Array  Transmitting 

The  material  presented  in  Section  3.1,  is  generalized  in  this 
section  for  the  multiple  antenna  array.  Voltages  and  currents  will  now 
become  vectors  of  voltages  and  currents.  The  reflection  and  transmission 
coefficients  will  become  matrices. 

The  impedance  model  of  a  multiple  antenna  array  transmitting  is 
given  in  Figure  9.  To  remain  consistant  with  the  rest  of  the  text,  we  shall 
assume  an  array  of  N  +  1  elements  and  use  the  indices  0,  1,  ...,  N.  Figure 
10  illustrates  the  various  transmission  line  parameters  while  the  array  is 
transmitting. 

The  transmitting  array  is  characterized  by  ST»  the  scattering 
matrix  seen  by  the  terminations  looking  into  the  transmission  lines. 

An  expression  for  the  matrix  can  be  calculated  from  Equation  (Cll)  of 
Appendix  C.2. 

ST  "  '  ZT][2A  +  ZT]_1  (3,13) 
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Figure  13  Scenario  Which  Removes  a  Degree  of 
Freedom  in  Conventional  Array 


above,  it  was  shown  for  a  parasitic  array  certain  combinations  of  inci¬ 
dent  signals  could  cause  the  loss  of  a  degree  of  freedom.  As  will  be 
discussed  in  the  following  section,  the  need  to  suppress  an  instability 
will  also  remove  a  degree  of  freedom. 

The  number  of  degrees  of  freedom  is  the  dimensionality  of 
the  usable  control  space.  Initially,  the  dimensionality  is  equal  to 
the  number  of  control  variables,  i.e.,  the  number  of  complex  weight  or 
variable  reflectivities.  In  the  derivation  above  it  was  assumed  the 
array  was  maximally  constrained,  that  the  N  control  variables  were  trying 
to  null  N  signals.  If  there  were  fewer  than  N  signals  to  null,  then  there 
would  not  be  particular  values  for  the  reflectivities  as  given  by  Equation 
(4.9),  but  a  range  of  suitable  values. 

The  reflectivity  values  necessary  to  form  nulls  can  be  cal¬ 
culated  from  Equation  (4.9) .  Before  cancellation  can  be  achieved,  how¬ 
ever,  one  must  also  determine  that  the  array  is  stable  for  those  reflec¬ 
tivity  values.  Stability  is  a  separate  issue  and  is  addressed  in  the 
following  section. 

The  issue  of  nulling  signals  at  different  frequencies  will  be 
briefly  discussed.  First  note  that,  at  the  different  frequencies  the 
mutual  coupling  coefficients  will  in  general  be  different.  Let  us  define 
T 

m  Q  ^  and  to  represent  the  mutual  coupling  factors  at  the  frequency  of 
the  nth  signal.  Rewriting  Equation  (4.1)  and  applying  the  constraint 
that  the  main  antenna  voltage  is  0  at  a  null,  we  obtain  the  Equations 


Note  these  equations  cannot  be  combined  as  was  done  before.  However, 
the  can  be  solved  for  and  substituted. 


0  '  ysoi  +=oiTp  11  ' 


sl 


0  ■  yS0K  +  20NTp  11  - 


(4.11) 


This  is  a  set  of  N  nonlinear  equations  in  N  unknowns,  p  j, . . . ,  The 

solutions  of  these  equations  are  not  easily  obtained. 

4.2  Stability  of  a  Parasitic  Array  With  Active  Complex  Terminations 

There  are  two  types  of  stability  that  can  be  analyzed  in  ref¬ 
erence  to  a  parasitic  adaptive  array,  RF  stability  and  control  stability. 
The  issue  of  HF  stability  is  independent  of  whatever  control  algorithm 
is  being  used.  RF  stability  determines  whether  the  system  is  stable  for 
a  given  set  of  reflectivities.  This  section  deals  with  RF  stability. 
Control  stability  is  the  issue  of  stability  and  convergence  of  the  algo¬ 
rithm  controlling  the  reflectivities.  Control  stability  is  discussed  in 
Chapter  5. 

It  was  pointed  out  in  Section  2.3,  that  given  a  set  of  reflec¬ 
tivity  values,  a  parasitic  array  is  linear  in  terms  of  plane  wave  inputs 
relative  to  an  output  to  the  receiver.  The  domain  of  RF  stability  is 
simply  those  sets  of  reflectivity  values  for  which  this  linear  system  is 
stable.  The  determination  of  that  domain  is  extremely  complex  and  shall 
be  accomplished  in  two  steps.  First  a  necessary  condition  for  stability 
is  derived  and  discussed.  This  condition  is  shown  to  be  sufficient  in 
certain  cases.  Several  examples  are  analyzed  to  give  insights  into  the 
behavior  of  an  actively  terminated  parasitic  array.  Then  a  general  suf¬ 
ficient  condition  is  derived. 

The  necessary  condition  is  derived  for  a  specific  frequency  of 
operation.  In  general  the  mutual  coupling  factors,  hl  . ,  will  be  functions 
of  frequency.  The  reflectivities  will  also  not  be  constant  over  frequency 
due  to  characteristics  of  the  electronics.  The  necessary  condition  re¬ 
quires  that  at  the  frequency  of  operation,  the  array  does  not  exhibit 
positive  feedback. 
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Let  us  first  study  a  single  antenna.  Figure  14a  illustrates  the 
signals  in  a  single  antenna.  Figure  14b  illustrates  these  signals  in  a 
system  block  diagram.  From  the  equation  for  simple  feedback  the  output 
is  computed  to  be  y  =  (1-mp)  ^  yg,  given  the  system  is  stable.  It  is  a 
well  known  result  that  a  necessary  condition  for  stability  is  Re[mp]  <  1. 
Stating  this  condition  another  way,  Re[mp]  >  1  is  a  sufficient  condition 
for  instability. 

Looking  at  the  complete  array  we  obtain  Equation  (2.1). 

L  e  ^s+  +  “M  <2-« 

As  long  as  M+P+  does  not  represent  positive  feedback  in  some  manner,  the 
array  will  be  stable.  To  study  the  multidimensional  feedback  represented 
by  M  P  it  will  be  convenient  to  introduce  the  concept  of  eigenvalues  and 
eigenvectors . 

For  a  matrix,  A,  there  are  certain  vectors  called  eigenvectors, 
that  when  operated  on  by  the  matrix  A  return  a  scaled  version  of  them¬ 
selves.  The  scale  factor,  X,  for  an  eigenvector,  x,  is  called  an  eigen¬ 
value.  Equation  (4.12)  illustrates  this  relationship  for  an  eigenvalue, 

X,  and  an  eigenvector,  x. 

Ax  =  Ax  (4.12) 

For  an  NxN  matrix  there  are  N  eigenvalues  not  necessarily  dis¬ 
tinct.  These  eigenvalues  are  found  as  the  roots  of  the  polynomial 
Det  (Xl-A) .  For  each  eigenvalue  there  is  an  eigenvector,  x*  determined 
to  within  a  scalar  factor.  If  the  eigenvalues  are  distinct  then  the  set 
of  eigenvectors  is  linearly  independent  and  in  fact  mutually  orthogonal. 
The  matrix  of  these  eigenvectors,  T  =  [xQ  |...|  xN]  ,  will  diagonalize  A. 

TAT_1  = 


If  the  eigenvalues  are  not  distinct,  then  the  set  of  eigenvectors  may 
or  may  not  be  linearly  independent.  If  they  are  not,  a  set  of  generalized 
eigenvectors  can  be  found  which  is  linearly  independent,  [16,  Chapter  2], 
The  matrix  of  these  generalized  eigenvectors,  in  a  manner  identical  to 
that  above,  will  transform  A  into  what  is  called  Jordan  form.  Jordan  form 


has  the  eigenvalues  on  the  main  diagonal,  l's  in  particular  spots  on 
the  diagonal  directly  above  it  and  0's  every  where  else.  In  the  dis¬ 
cussion  that  follows,  for  computational  ease,  it  will  be  assumed  that 
M+P+  has  distinct  eigenvalues.  The  results  derived,  however,  can  be 
shown  valid  for  general  matrices  by  considering  generalized  eigenvectors 
and  their  accompanying  Jordan  forms  [16] . 

Let  Ag  ,  ...,1^^  be  the  N  distinct  eigenvalues  of  M+P+  and 
the  corresponding  eigenvectors.  Finally,  let  T  =  [xQ  |  —  |x^  ] 
Since  the  eigenvectors  are  linearly  independent,  the  vector  of 
Equation  (2.1)  can  be  expressed  as  a  weighted  sum  of  them. 


N 


r+  ‘  E  ai-i 

i=0 


=  Ta  where  a  = 


(4.13) 


From  Equation  (4.13) ,  the  weighting  value  for  x.  is  seen  to  be  the  ith 

-1  1 

component  of  the  vector  T  Similarly,  can  also  be  written  as  a 

weighted  sum  of  the  eigenvectors. 


(4.14) 


With  Equations  (4.13)  and  (4.14)  one  can  rewrite  Eaiation  (2.1)  as  below. 
N  N  N 


a.x.  =  b.x,  +  MP  a.x. 


i=0 


x—  1 


.  _  i—i  +  +  .  _  i—  i 

1=0  1=0 


E 

i=0 


(b.  +  A. a . )  x  . 
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(4.15) 


We  see  that  Equation  (4.15)  really  represents  N  +  1  independent  equations 
of  the  form, 


Here  the  N  +  1  dependent  variables  y0,...,yN  have  been  trans- 
formed  into  the  independent  variables  a„,...,a„.  Thus  the  feedback  nature 
of  the  array  can  be  expressed  by  N  +  1,  independent  subsystems  of  the 
form  illustrated  in  Figure  14b,  where  the  loop  gain  is  for  the  ith  sub¬ 
system.  We  can  apply  the  previously  derived  stability  condition  to  each 
of  these  subsystems  and  thereby  generalize  the  result  for  an  N  +  1  element 
array.  Given  the  array  parameters,  m^  and  p^,  as  measured  for  a  par¬ 
ticular  frequency,  one  can  find  N  +  1  eigenvalues  of 

M+P  .  Then  a  necessary  condition  for  stability  is  Re[X^]  <  1  for  all 
i  *=  0, . • . »N.  As  previously  stated,  this  result  can  be  shown  valid  for 
general  M+P+. 

With  a  parasitic  array,  a  useful  and  practical  constraint  is 
often  that  the  reflectivity  of  the  main  antenna  be  zero.  This  means 
the  receiver  is  well  matched  to  the  transmission  line.  Under  this 
condition  one  eigenvalue  of  M+P+  is  necessarily  0,  and  the  rest  are  the 
eigenvalues  of  MP.  Thus  if  p  ■  0,  one  need  only  consider  the  array  of 
auxiliary  elements  in  a  stability  analysis. 

It  is  shown  in  Appendix  E  if  the  reflectivities  contain  a  single 
stage,  narrow  bandpass  filter,  such  that  the  mutual  coupling  coefficients 
appear  constant  over  the  passband,  then  this  stability  condition  is  both 
necessary  and  sufficient.  Such  a  filter  was  designed  into  the  active 
terminations  built  by  ZAI. 

Under  the  assumption  that  a  single  stage  narrowband  filter  is 
in  use,  one  can  perform  some  insightful  stability  analysis.  In  Appendix  F 
a  stability  analysis  is  made  with  regard  to  a  two  element  array . 

Consider  an  example  of  a  two  element  array  of  auxiliary  antennas. 

It  is  assumed  that  the  first  antenna  is  well  matched  to  its  transmission 
line  so  that  its  self  coupling,  or  return  loss,  is  zero.  It  is  also 
assumed  that  the  reflectivity  and  self  coupling  of  the  second  antenna  is 
such  that  by  itself  it  would  be  unstable,  i.e.,  ^[“*22^2^  1  X.  It  is 
shown  in  Appendix  F  under  certain  conditions,  Re[m22P2]  <  2,  that  an  appro¬ 
priate  choice  of  p^  will  stabilize  the  array. 

The  most  general  stability  criterion  is  obtained  from  [16,  p.  376]. 
Assume  that  the  mutual  coupling  matrix  M+(s)  is  known  as  a  function  of 
complex  frequency,  s.  In  particular,  that  each  entry  of  M+(s)  is  a  rational 
polynomial  in  s.  Similarly,  assume  that  the  diagonal  entries  of  P+(s)  are 


.■  * ■  j 


rational  polynomials  in  s.  Let  the  polynomial  A(s)  be  the  least  common 
denominator  of  all  the  minors  of  all  orders1  of  M+(s)P+(s).  Then  the  arrav 
is  stable  if  and  only  if  det[l  -  M+(*)P+(“) ]  /  0  and  the  polynomial 
A(s)det[l  -  M  (s)P  (s)]  has  no  right  half  plane  zeros. 

T  T 

In  Appendix  E,  it  is  shown  from  this  general  criterion,  that  if 
f(s)  is  a  single  stage  narrowband  filter,  then  the  necessary  stability 
condition,  RetX^]  <  1  ,  is  both  necessary  and  sufficient. 


A  minor  of  the  kth  order  is  formed  by  taking  the  determinant  of 
all  the  elements  common  between  k  distinct  rows  and  k  distinct  columns. 
A  minor  of  order  1  is  an  element.  Minors  of  order  n-1  are  the  minors 
people  are  most  acquainted  with.  These  are  formed  by  "blocking  out" 


a  row  and  a  column  and  then  taking  the  determinant.  There  are 

number  of  minors  of  the  kth  order,  where  ^  y  . 

Note,  A(s)  is  not  necessarily  the  denominator  of  the  determinant, 
the  multivariate  feedback  represented  by 


Consider 


G  (s) 


1  -s 

(s+1) (s-1)  (1+s) 

1 

(s-1)  -1 


A(s)  =  (s+1) (s-1)  . 


Then  detG(s) 


1 

s+1 


but 
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5.0  ADAPTIVE  CONTROL  OF  PARASITIC  ARRAY  NULLS 

For  a  parasitic  array  to  be  useful  in  practice,  adaptive  control 
of  the  terminations  or  reflectivities  must  be  employed.  There  are  several 
adaptive  control  algorithms  and  several  ways  of  implementing  each  algorithm. 
He  shall  discuss  the  gradient  algorithm  and  various  search  and  dither 
algorithms  and  also  possible  methods  of  implementing  each. 

The  gradient  algorithm  would  be  the  fastest  and  probably  achieve 
deeper  nulls.  The  gradient  tells  not  only  when  to  change,  but  also  in 
which  direction.  It  will,  however,  use  more  hardware  and  cabling  and  may 
require  a  training  pilot  signal. 

There  are  several  possible  dither  algorithms,  some  of  them  closely 
related  to  the  gradient  algorithm  and  some  more  akin  to  a  search  algorithm. 
Performance  and  hardware  needs  will  vary  greatly  depending  on  the  algorithm 
find  its  implementation.  The  best  compromise  of  performance  over  hardware 
might  be  some  kind  of  dither  algorithm. 

Search  algorithms  generally  will  be  slower.  They  may  achieve  as 
deep  of  nulls.  Search  algorithms,  however,  will  require  the  least  amount 
of  hardware  and  cabling.  The  only  connections  to  the  terminations  needed 
are  DC  control  and  power. 

5.1  Least  Mean  Square  Control  of  Active,  Complex  Terminations 


In  the  search  for  an  adaptive  control  algorithm  of  a  compact 
array,  one  might  first  ask  why  not  take  the  outputs  of  the  antennas  and 
put  them  into  a  conventional  adaptive  array  combiner  as  in  Figure  1.  One 
could  use  such  a  system,  but  there  are  a  few  drawbacks. 

First,  the  practical  variable  RF  weight  has  a  high  reflection 
coefficient  for  much  of  its  operating  range.  Such  weights,  unbuffered 
with  the  strong  mutual  coupling  of  a  compact  array  would  have  an  unpredictable 
and  possibly  uncontrol lable  effect  on  the  array.  To  buffer  the  antennas  from 
the  weights,  one  would  have  to  take  a  high  degree  of  insertion  loss  or 
include  active  elements  like  an  amplifier. 

Second  and  more  importantly,  the  mutual  coupling  can  cause  blind 
spots  in  the  array  pattern  [19].  The  effect  of  these  blind  spots  might  be 
the  array  loses  some  degrees  of  freedom  for  jammers  arriving  from  particular 
directions,  or  there  might  be  some  direction  which  is  always  nulled  and  the 
desired  signal  cannot  be  received  from  that  direction.  Note  the  effects  of 
these  blind  spots  would  be  independent  of  the  weight  values,  assuming  the 
weights  have  been  buffered. 
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In  a  parasitic  array,  terminated  in  variable  impedances,  blind 
spots  due  to  the  mutual  coupling  cure  not  as  great  a  problem.  One  has 
control  over  them  through  the  variable  terminations.  It  is  possible,  as 
was  described  in  Section  4.1,  that  certain  combinations  of  jammers  can 
cause  the  loss  of  a  degree  of  freedom  of  operation,  but  for  an  array  of 
several  elements  this  is  highly  unlikely. 

5.2  Gradient  Control 

Adaptive  control  of  a  conventional  array  is  a  linear  problem. 
The  change  in  the  output  of  the  array  is  proportional  to  the  change  in 
any  of  the  complex  weights.  Adaptive  control  of  a  parasitic  array,  how¬ 
ever,  is  a  nonlinear  problem.  How  the  array  output  changes  with  a  change 
in  a  reflectivity  will  depend  not  only  on  the  current  value  of  that 
reflectivity,  but  also  on  the  values  of  all  the  other  reflectivities. 

A  generalization  of  the  LMS  algorithm  for  nonlinear  problems 
is  the  gradient  control  algorithm.  The  gradient  control  law  for  complex 
variables  is  derived  in  Appendix  G.2.  The  control  law  is  stated  in 
Equation  (G17) .  This  control  law  is  applied  to  adaptive  control  of 
parasitic  termination  in  Appendix  G.3.  The  result  is  given  in  Equation 
(G33)  . 


y0  lc,n  V* 


where  a. 


T  -1 
e  Q  in 
—  n  —  0 


T  -1 

1  -  p„  (ni  +  m  _  PQ  m  _) 
K0  00  —  0  —  0 
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and  e 
—  n 


nth 

place 


The  two  expressions  for  can  be  shown  to  be  equivalent.  The  first  is 

more  convenient  for  calculations.  A  physical  understanding  can  be  gleaned 

more  easily  from  the  second  form.  Consider  the  factor  ay  =  3y_/3p  •  It  is 

n  n  0  n 

the  proportionality  constant  between  a  small  change  in  p^  and  the  resultant 
change  in  y^.  Let  us  assume  a  small  change  Ap^,  in,  p^  and  see  how  it 
changes  yQ.  For  the  moment,  also  assume  that  *  0  so  there  are  no 


5-2 


One  can  consider  the  initial  change  of  y  Ap  in  the  upward  travei- 

n  n 


reflections  off  the  main  receiver.  The  change  Ap  in  p  changes  the  upward 

n  n 

traveling  wave  (toward  the  antenna)  reflected  of.f  the  nth  termination  by 

y  Ap  .  This  signal  is  coupled  to  the  main  antenna  by  the  factor  m.  .  Thus 
n  n  „  On 

T 

y.  changes  by  the  amount  m.  y  Ap  .  Let  us  write  this  as  m„  e  y  Ap  .  But 
this  is  not  the  end  of  the  story.  The  signal  y^Apn  also  couples  to  the 
other  antennas,  changes  their  voltages,  and  this  change  couples  to  the  main 
antenna 

ing  wave  at  the  nth  antenna  as  an  initial  change  in  the  vector  of  upward 

traveling  voltages  of  e  y  Ap  .  This  signal  couples  to  the  other  auxiliary 

— n  n  n 

antenna  by  the  coupling  matrix  H.  The  initial  vector  of  downward  voltages 

due  to  a  change  Ap  in  p  is  then  Me  y  Ap  .  This  is  like  a  change  in  y  . 

n  n  ~ n  n  n  ^  “-s 

The  final  change  in  the  downward  voltages,  y,  is  then  Q  Me  y  Ap  .  This 

— n  n  n 

signal  gets  reflected  by  P  and  coupled  to  the  main  antenna  by  m^.  Thus 

the  change  in  y  due  to  the  change  in  the  auxiliary  voltages  caused  by 
0  T  -1 

the  change  in  p  is  PQ  Me  y  Ap  .  If  we  sum  these  two  terms  and  divide 
by  Ap^,  we  obtain  the  numerator  of  the  second  expression  for  a^. 

If  p0  t*  0,  this  creates  a  feedback  loop  to  the  main  antenna 
through  the  self  coupling  “>00»  and  the  coupling  to  the  auxiliaries  n^. 

The  loop  gain  of  this  feedback  is  the  sum  of  two  terms.  The  first  term. 

The  second  term  takes  the  gain  of  the 

As  above,  this 


from  the  self  coupling,  is  p^. 

reflectivities  pQ  and  couples  it  to  the  auxiliaries  by  m^. 
signal  can  be  considered  an  initial  signal  and  is  turned  into  a  final  sig¬ 
nal  by  Q  1.  This  final  signal  gets  reflected  by  P  and  coupled  back  to 

the  main  antenna  again  by  m.  completing  the  loop.  Thus  the  gain  of  the 

T  -1  ~ 0  T  -1 

second  loop  is  m^  PQ  “IqPq*  The  total  loop  gain  is  P^®^  +  2,,  PQ  n^)  • 


0’  00 


We  see  that  the  denominator  of  a  represents  positive  feedback  with  this 


loop  gain. 


One  might  ask  if  the  feedback  loop  is  stable.  It  can  be  shown 
with  the  formula  for  the  determinant  of  a  partitioned  matrix  (Identity  A. 10) 
that  this  denominator  is  an  eigenvalue  of  Q+  or  equivalently  the  loop  gain 
is  an  eigenvalue  of  M+P+.  The  stability  conditions  for  this  loop  are  there¬ 
fore  a  special  case  of  the  stability  conditions  described  in  Section  4.2. 

Gradient  control  can  be  implemented  with  basically  the  same  equip¬ 
ment  as  LMS  control,  except  a  method  of  obtaining  or  calculating  the  factor 
is  needed.  A  basic  implementation  is  illustrated  in  Figure  15.  Note  the 
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Figure  15  Basic  Gradient  Control  Implementation 


main  antenna  signal  must  be  sent  to  the  control  circuitry  of  each 

auxiliary.  Also  note  that  an  extra  complex  weight  will  probably  be 

needed  to  implement  the  multiplicative  factor  a  ,  or  a  depending  on 

n  n 

where  it  is  placed. 

In  implementing  gradient  control  in  this  manner,  it  is  only 

necessary  to  determine  the  phase  of  a  •  The  amplitude  of  a  appears 

n  n 

simply  as  a  modification  to  the  control  loop  gain  k.  Furthermore,  small 
errors  in  the  estimation  of  the  phase  of  are  not  significant.  Theo¬ 
retically,  the  control  loop  will  still  work  with  up  to  89°  of  phase  error. 
Phase  errors  manifest  themselves  by  reducing  the  loop  gain  by  a  factor 

of  the  cosine  of  the  error.  In  practice,  if  one  is  not  too  concerned 

o 

with  settling  time,  phase  errors  of  up  to  40  can  be  tolerated.  All 

this  shows  that  whatever  process  is  used  to  estimate  the  factor  a  ,  it 

n 

need  not  be  extremely  accurate. 

Three  methods  of  obtaining  the  factor  fire  presented.  The 
first  method  uses  a  pilot  signal,  close  in  frequency,  sent  out  of  the 
main  antenna.  In  appendix  H,  it  is  shown  that  such  a  pilot  signal  (with 
the  condition  PQ  =  0)  will  be  received  at  the  nth  auxiliary  antenna  mul¬ 
tiplied  by  the  factor  a  .  This  is  measured  and  then  used  to  complete  the 

n 

gradient  algorithm.  The  second  method  is  to  determine  all  the  mutual 
coupling  factors  of  the  array  over  all  frequencies  as  well  as  how  the  fil¬ 
ter  characteristics  of  the  reflectivities  vary  over  frequency.  This  in¬ 
formation  would  be  stored  and  used  to  digitally  calculate  the  necessary 

values  of  the  a  .  These  values  would  then  be  implemented  with  a  digitally 
II 

controlled  weight.  The  third  method  uses  a  dither  to  approximate  which 
is  then  used  in  a  gradient  control  system.  This  method  also  requires 
pg  =  0.  The  description  of  this  algorithm  is  included  in  this  section  be¬ 
cause  it  implements  gradient  control  as  expressed  in  Equation  (5.1) . 

Figure  16  illustrates  an  analog  implementation  of  the  pilot  al¬ 
gorithm  with  the  factor  determined  and  applied  by  an  adaptive  weight. 
There  are  several  points  to  make  about  this  design.  First,  we  give  a  de¬ 
scription  of  the  general  operation.  The  reflectivity  at  the  auxiliary 
antenna  is  implemented  in  polar  form,  as  is  currently  being  done  on  the 
NHL  test  array.  The  phase  control  is  obtained  by  reflecting  the  signal 
off  a  variable  reactance.  The  amplitude  control  is  obtained  with  a  bi- 


Figure  16  Design  for  SIC  with  Pilot  and  Polar  Reflectivities 


of  the  reactances  while  observing  the  receiver  output  responding  to  one  or 
more  signals  incident  on  the  array.  That  is,  the  pattern  was  formed  adaptively. 
Some  preliminary  attempts  were  made  to  use  the  array  in  a  deterministic  manner 
by  using  the  theory  developed  by  Harrington  [6]  to  calculate  the  set  of 
reactance  values  that  comes  closest  to  synthesizing  a  given  pattern.  Using 
a  commercial  RF  bridge  to  measure  the  reactance  value,  the  reactive  termin¬ 
ations  were  then  adjusted  to  have  the  corresponding  calculated  reactance. 
However,  the  approximations  in  the  theory  and  the  inability  to  measure  the 
reactance  values  sufficiently  accurately  usually  resulted  in  a  measured 
pattern  that  differed  substantially  from  the  desired  pattern.  Therefore, 
we  ended  all  efforts  to  use  the  array  in  a  deterministic  manner. 


6.2  Experimental  Measurements 

The  manual  adaptation  investigations  were  performed  using  two  tech¬ 
niques  to  adjust  the  reactance  values.  In  the  first  technique  a  single  CW 
source  was  installed  approximately  70  m  from  the  array.  The  center  element 
was  connected  to  a  HF  receiver,  tuned  to  the  source,  whose  output  was  observed 
on  an  oscilloscope.  We  then  adjusted  the  six  terminating  reactances  sequen¬ 
tially  and  iteratively  to  produce  a  minimum  in  the  receiver  output.  The 
azimuthal  antenna  pattern  was  then  measured.  In  Figure  22  we  show  a  typical 
pattern  formed  in  this  manner.  There  are  two  noteworthy  features.  First,  a 
sharp  null  has  indeed  been  formed  in  the  direction  of  the  source,  as  desired, 
that  has  a  depth  approximately  30  dB  below  the  pattern  main  lobe,  and  that 
has  a  width  of  about  10  degrees.  The  remainder  of  the  pattern  is  relatively 
featureless;  except  for  the  small  dip  near  240  degrees,  the  pattern  away  from 
the  null  varies  no  more  than  3  dB. 

The  second  noteworthy  feature  is  that  the  gain  of  the  array  in 
azimuthal  directions  away  from  the  null  is,  on  the  average,  only  4.7  dB  less 
than  the  gain  for  a  monopole  resonant  at  20  MHz.  In  the  sector  from  30 
degrees  to  180  degrees,  the  difference  decreases  to  2.4  dB.  Hence,  the 
sensitivity  of  the  array  is  nearly  that  of  a  monopole  and  does  not  suffer  the 
large  degradation  that  one  might  expect  from  such  a  small  array.  Patterns 
similar  to  Figure  10  (that  is,  an  essentially  omnidirectional  pattern  with  a 
notch  in  the  direction  of  an  incident  signal)  could  be  formed  for  an 
incident  signal  of  any  bearing. 
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Figure  21  Antenna  Array  Test  setup 
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6.0  EXPERIMENTAL  RESULTS  WITH  A  PASSIVE  REACTIVELY  TERMINATED  PARASITIC  ARRAY 

NRL  has  built  an  electrically  small  parasitic  array  for  experimental 
testing.  At  first  the.  elements  of  the  array  were  terminated  with  only  variable 
reactances.  Test  results  with  these  terminations  are  detailed  in  [2-5]  and 
are  summarized  here.  Tests  with  active  complex  terminations  are  discussed 
in  Chapter  8. 

6.1  Experimental  Parasitic  Array 

A  seven  element  array  was  fabricated  and  tested  by  NRL  to  evaluate 
the  pattern  control  achievable  with  reactively  terminated  parasitic  arrays. 
Figure  19  is  a  photograph  of  the  array.  This  array  has  one  central  element 
connected  to  a  receiver  and  six  parasitic  elements.  All  elements  are  uniformly- 
wound  helices  with  a  5  cm  diameter  and  an  overall  height  of  one  meter.  The 
antennas  are  constructed  to  have  a  resonant  frequency  of  about  20  MHz, 
and  the  spacing  between  the  central  element  and  the  surrounding  parasitic 
elements  is  40  cm  (0.027  X  at  20  MHz) .  The  top  plate  visible  in  Figure  19 
is  non-metallic  and  is  solely  for  rigid  support  of  the  elements. 

The  electronically-controllable  variable  reactances  are  mounted  in 
the  boxes  at  the  base  of  each  parasitic  element.  A  schematic  diagram  of 
the  reactance  is  given  in  Figure  20.  Continuous  control  of  the  reactance 
termination  is  achieved  by  adjusting  the  bias  voltage  on  a  varactor  diode; 
the  continuous  control  range  mid-point  is  determined  by  the  combination  of  the 
series  inductors  switched  into  the  reactance  circuit.  At  the  20  MHz  resonant 
frequency  of  the  elements,  the  terminating  reactance  can  be  varied  over  a 
range  from  -450  ohms  (capacitive)  to  +450  ohms  (inductive) . 

The  terminations  shown  in  Figure  20  provide  control  only  of  the 
phase  of  the  incident  waves  reflected  from  a  parasitic  element.  Amplitude 
control,  using  variable  attenuators  alone  or  in  conjunction  with  amplifiers, 
had  not  yet  been  implemented  in  the  array.  There  is  a  similarity  between  the 
parasitic  array  described  here  and  various  "phase-only"  adaptive  arrays  de¬ 
scribed  in  the  literature  [ 22] . 

Measurements  of  the  array  antenna  pattern  were  made  at  NRL's 
Brandywine  Antenna  Range  (see  Figure  21) .  The  antenna  range  includes  a  large 
ground  plane  and  a  rotating  platter  for  measuring  the  azimuthal  radiation 
pattern.  All  of  the  measurements  reported  here  were  made  by  manual  adjustment 
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accelerated  rAndom  search  (GAPS)  algorithm  is  discussed  in  the 
references  [22].  A  guided  random  search  always  has  the  possibility  of 
choosing  a  new  point  far  away  from  the  current  operating  region.  At 
first  glance  this  aspect  may  seem  a  hinderance  to  convergence,  but 
because  of  this  property  a  guided  random  search  can  only  get 
temporarily  hung  up  at  a  local  extrema. 


it  may  then  go  into  a  tracking  mode. 

We  shall  discuss  two  types  of  random  search  algorithms,  the 
random  search  and  the  guided  random  search.  These  algorithms  are 
somewhat  akin  to  the  exhaustive  and  guided  univariate  search  algorithms 
discussed  above.  We  say  somewhat  because  the  guided  random  search  is 
not  univariate. 

The  random  search  assigns  a  probability  distribution  to  the 
whole  control  space.  It  may  be  uniform  or  it  may  be  shaped  by  certain 
apriori  information  of  the  control  system  and/or  the  optimal  point. 

The  algorithm  then  selects  an  operating  point  at  random  based  on  this 
distribution  and  evaluates  its  PA.  If  the  PA  sastifies  the  threshold, 
then  it  stops.  Otherwise  it  obtains  another  randomly  generated 
operating  point.  Random  search  algorithms  can  be  run  with  or  without 
replacement.  In  an  algorithm  with  replacement  the  probability 
distribution  is  unchanged  from  one  trial  to  the  next.  In  an  algorithm 
without  replacement,  if  an  operating  point  selected  in  a  trial  does 
not  pass  the  threshold  test,  then  that  point  is  removed  from  the  control 
space  and  the  distribution  is  renormalized.  If  it  is  a  continuous 
distribution,  then  a  region  around  the  point  is  removed.  Generally, 
one  would  want  to  use  an  algorithm  without  replacement  unless  the 
system  is  not  time-invariant  over  the  time  it  takes  to  converge.  The 
random  search  without  replacement  is  very  much  like  the  exhaustive 
search,  except  the  random  search  stops  when  it  finds  a  satisfactory 
operating  point.  It  is  conceivable,  however,  that  the  random  search 
would  search  through  all  the  points  before  finding  the  right  one. 

Just  as  the  univariate  and  guided  univariate  searches  were 
attempts  to  expedite  the  exhaustive  search,  the  guided  random  search 
is  an  attempt  to  expedite  the  random  search.  In  the  guided  random  search 
the  probability  distribution  described  above  is  used  for  the  first  trial. 
To  choose  a  point  for  a  second  trial  a  random  direction  from  the  first 
point  is  selected.  Then  a  random  distance  out  in  this  direction  is 
selected  from  a  distribution  whose  mean  might  be  the  interval  used  in 
the  guided  univariate  search  above.  After  the  second  point,  one  now 
squints  the  direction  probabilities  distribution  either  in  the  direction 
of  the  second  point  or  in  the  opposite  direction,  depending  on  whether 
the  PA  was  better  or  worse.  One  can  also  accelerate  a  guided  random 
search  just  as  the  guided  univariate  search  was  accelerated.  A  guided 


completed  with  it,  etc.  This  algorithm  has  advantages  over  the  ex¬ 
haustive  search.  First,  as  the  algorithm  procedes  the  output  becomes 
more  and  more  minimized.  He  reap  the  benifits  of  previous  computations 
as  we  go  along.  Second,  the  algorithm  converg  s  faster,  particularly 
if  the  array  is  underconstrained.  In  practice,  usually  after  several 
iterations  negligible  further  increase  in  performance  is  observed. 

This  is  the  method  NRL  is  currently  using  in  their  computer  run  adaptive 
algorithm. 

A  third  search  algorithm  is  the  guided  univariate.  This  al¬ 
gorithm  is  the  same  as  the  univariate  except  an  attempt  is  made  to  mini¬ 
mize  the  number  of  operating  points  analyzed  within  each  iteration. 

When  the  algorithm  begins  a  new  iteration  of  a  control  variable,  a  point, 
a  certain  interval  away  from  current  operating  value  is  chosen.  If  the 
output  at  this  point  is  less  than  at  the  initial  point,  the  algorithm 
continues  in  this  direction  until  the  output  stops  decreasing.  The 
algorithm  then  zeros  in  on  the  local  optimal  value  of  this  variable. 

If  the  output  of  the  first  guess  was  higher  than  at  the  initial  point, 
the  algorithm  continues  in  the  opposite  direction  until  the  output  stops 
decreasing  and  then  zeros  in  on  the  local  optimum.  A  way  to  reduce  the 
number  of  operating  points  evaluated  is  to  make  the  step  sizes  between 
the  evaluated  points  variable.  If  a  big  change  in  the  output  was  noticed 
between  the  last  two  evaluations,  then  make  a  big  step  in  the  control 
variable.  If  only  a  small  change  was  noticed,  then  make  a  small  step. 
This  is  called  accelerating  the  algorithm.  The  guided  univariate 
search  is  a  more  complex  algorithm  than  the  univariate  search,  but  it 
will  converge  faster  because  fewer  operating  points  are  evaluated. 

One  problem  that  both  the  univariate  search  and  the  guided 
univariate  search  have  is  that  they  can  get  hung  up  on  a  local  minimum 
of  the  output,  or  a  local  maximum  of  S/I.  Random  search  algorithms 
may  temporarily  get  hung  up  on  a  local  output  minimum,  but  will 
eventually  find  the  global  optimal  operating  point  or  at  least  a 
satisfactory  operating  point. 

In  all  random  search  algorithms,  one  computes  a  performance 
assessment  (PA)  statistic  for  each  new  operating  point.  In  our  case 
this  might  be  the  array  output  power  (to  be  minimized)  or  the  S/I  (to 
be  maximixed) .  Whenever  the  algorithm  finds  a  point  for  which  the  PA 
is  less  than  (if  minimizing)  or  greater  than  (if  maximizing)  a 
particular  threshold,  then  the  algorithm  stops  searching.  If  need  be, 


5.4  Search  Algorithms 


Search  algorithms  attempt  to  search  the  control  space  for  the 
optimal  operating  point.  Given  the  value  of  the  function  to  be  minimized 
at  the  current  operating  point  and  the  past  history  of  the  function  at 
previous  operating  points,  a  search  algorithm  makes  a  guess  as  to  where 
the  optimal  operating  point  might  be.  All  search  algorithms  are  sequential, 
i.e.,  they  choose  an  operating  point,  evaluate  it  and  then  choose  another, 
etc.  Unless  the  evaluation  and  selection  of  operating  points  is  done 
manually,  search  algorithms  must  be  implemented  digitally. 

An  advantage  of  search  algorithms  is  they  require  mi  nim»i 
amount  of  hardware  of  all  the  algorithms,  neglecting  the  computer .  other 
than  the  computer,  they  require  only  D.C.  control  lines  to  the  termina¬ 
tions  and  a  signal  to  interference  ratio  (S/I)  detector.  There  is  no 
need  for  correlators,  adaptive  control  loops,  pilots,  or  extra  complex 
weights.  Furthermore,  search  algorithms,  because  of  their  digital  im¬ 
plementation,  are  more  flexible.  If  one  desires  to  modify  the  algorithm, 
he  need  only  change  the  computer  program.  The  main  disadvantage  of 
search  algorithms  is  that  they  tend  to  be  slow. 

There  are  many  types  of  search  algorithms.  We  shall  describe 
several  here.  The  simplest  in  theory  is  the  exhaustive  search  algorithm. 

In  this  algorithm  every  operating  point  is  analyzed.  The  one  which  gives 
optimal  performance  is  noted  and  then  applied.  In  most  systems  the  num¬ 
ber  of  possible  operating  points  is  too  great  for  this  algorithm  to  be 
practical.  This  is  the  case  here.  Consider  an  array  two  complex  ter¬ 
minations  controlled  by  two  8  bit  signals  each.  Thus  the  number  of  oper- 
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ating  points  to  evaluate  is  2  raised  to  the  fourth  power.  This  is  2"**  or 
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approximately  4  x  10  .  If  we  could  evaluate  an  operating  point  every 
micro-second,  it  would  take  over  an  hour  to  evaluate  them  all.  The  one 
advantage  of  the  exhaustive  search  is  that  it  will  always  select  the 
optimal  operating  point. 

A  second  type  of  search  algorithm  is  the  univariate  search 
algorithm.  In  the  univariate  search  one  of  the  control  variables  is  selected. 
An  exhaustive  search  is  conducted  upon  this  variable  along  with  all  the 
other  variables  held  constant.  The  value  of  this  variable  which  gave 
the  minimum  output  (or  maximum  S/I)  over  its  range  is  then  applied  to 
it.  A  second  variable  is  then  selected  and  a  similar  procedure  is 


Most  dither  algorithms  attempt  to  measure  the  partial  derivative 
of  the  function  to  be  minimized  with  respect  to  some  other  variable  of  the 
system.  The  measurement  procedure  of  the  partial  derivative  is  taken 
straight  from  the  definition  of  the  partial  derivatives  the  ratio  of  a 
small  change  in  output  to  the  small  change  in  the  input  that  caused  it. 

A  particular  variable  of  the  system  is  "dithered"  about  its  present  value 
with  a  small  amplitude  A.C.  signal,  e.g.  a  sine  wave,  square  wave,  or  PN 
coded  signal.  The  output  signal,  or  the  signal  to  minimized,  is  correlated 
with  the  A.C.  signal.  The  correlator  output  is  the  partial  derivative. 

The  distinction  being  made  here  between  dither  algorithms  and  gradient 
algorithm  is  that  a  gradient  algorithm  uses  a  coherent  (i.e.  inphase  and 
quadrature)  correlator,  while  a  dither  algorithm  simplex  detect  the 
change  in  the  array  output  power. 

The  dither  algorithm  can  be  implemented  in  an  analog  or  digital 
manner.  In  an  analog  implementation  the  introduction  of  the  A.C.  dither 
signal  and  correlation  with  the  output  would  be  accomplished  with  analog 
components .  The  subsequent  adjustment  of  the  undithered  value  of  the  control 
variable  could  also  be  accomplished  with  analog  circuitry  driven  by  the 
correlator  outputs.  In  a  digital  implementation  the  dither  of  the  control 
variable  and  measurement  of  the  change  in  the  output  could  all  be  accomplished 
and  controlled  digitally.  For  example  if  the  control  variable  is  a  digitally 
controlled  weight,  the  dither  could  be  implemented  by  alternating  between 
to  adjacent  settings. 

Consider  the  case  where  the  various  control  variables  are  dithered 
individually  and  then  all  variables  are  changed  simultaneously  to  new  values. 
The  amount  that  each  variable  is  changed  is  proportional  to  the  amount  the 
output  changed  during  the  dither  for  that  variable,  i.e.  the  correlator 
output.  Such  an  algorithm  implements  sequential  gradient  control. 

The  advantage  of  dither  algorithms  over  gradient  algorithms  is 
the  reduction  in  the  amount  and  complexity  of  the  analog  hardware.  In 
fact,  besides  the  variable  weights,  a  dither  algorithm  could  be  implemented 
with  a  computer,  a  power  detector,  an  A/D  convertor,  and  2xN  D/A  convertors. 
The  cost  for  this  reduction  in  hardware  is  an  increase  in  the  adaption  time 
for  the  array  and  under  some  conditions  less  mull  depth.  In  a  changing 
RF  environment  the  array  may  never  completely  adapt. 


Figure  17  Design  for  SIC  with  Dither  Driven  Gradient 


measurement  errors  get  compounded  into  significant  final  errors  by  the  many 
mathematical  operations  of  the  computation.  Here  again  the  limitation  is 
more  severe  in  a  larger  array  where  there  will  be  more  operations  in  the 
calculation . 

A  third  implementation  of  the  gradient  algorithm  is  a  cross  be¬ 
tween  a  gradient  algorithm  and  a  dither  algorithm.  The  algorithm  is  based 
on  the  following  derivative 
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for  n  =  1,  ...,  N. 

These  derivatives  are  derived  in  a  manner  similar  to  that  used  to  derive 
3yQ/3p^.  The  derivatives  are  approximated  by  dithering  pQ  about  pQ  *  0 
and  monitoring  the  changes  in  y^.  A  design  for  an  implementation  of  this 
method  is  given  in  Figure  17.  The  reflectivity  pQ  is  dithered  by  alternately 
terminating  the  appropriate  port  on  the  directional  coupler  with  a  short  and 
open  circuit.  The  remaining  circuitry  is  similar  to  that  used  in  the  pilot 
implementation  to  estimate  a .  .  The  dither  basically  replaces  the  pilot 
signal.  The  advantage  of  this  method  is  that  one  dither  signal  allows 
the  calculation  of  all  the  gradient  law  constants,  a^.  There  is  no  need  for 
multiple,  mutually  orthogonal  dither  signals. 

One  possible  disadvantage  to  this  method  is  that  dithered  components 
of  the  desired  signal  may  enter  the  receiver.  For  some  combinations  of  desired 
signal  modulation  and  dither  signal  this  may  corrupt  reception.  If  the 
frequency  of  the  fundamental  component  of  the  dither  signal  is  greater  than 
the  information  bandwidth  of  the  signal  than  this  distortion  will  be  minimal. 

If  the  amount  of  distortion  is  too  great  an  extra  adaptive  loop  at  the  main 
antenna  can  be  used  as  in  Figure  18.  This  loop  will  cancel  the  dither 
component  that  might  enter  the  receiver.  The  adaption  of  this  loop  cam 
proceed  concurrently  with  that  of  the  array,  but  will  fully  adapt  only 
after  the  array  has  acquired.  The  adaption  of  the  array,  however,  can 
proceed  completely  independently  of  this  main  antenna  loop  (assuming  p  “  0)  . 

The  gradient  can  also  be  obtained  by  measuring  the  individual 
partial  derivatives  directly.  Most  dither  algorithms  take  something  of  this 
approach.  These  algorithms  are  discussed  in  the  next  section. 


polar  attenuator  in  the  return  path  of  this  reflected  signal.  The  factor 
necessary  for  the  gradient  algorithm  is  implemented  by  multiplying  the 
main  antenna  signal,  y^,  by  before  the  correlators  that  drive  the  com¬ 
plex  reflectivity.  This  is  exactly  what  is  illustrated  in  Figure  15, 
except  now  the  reflectivity  is  implemented  in  polar  form.  Note  for  a  polar 
weight  the  reference  signal  for  the  correlator  must  come  from  a  point  after 
the  variable  phase  control. 

In  the  design  of  Figure  16  a  second  adaptive  complex  weight  is  used 
to  both  determine  (actually  a*)  and  multiply  y^  by  it.  This  is  done 
simultaneously  because  the  input  to  the  weight  is  a  sum  of  the  two  signals,  y^ 

and  y  .  The  weight  and  its  control  are  structured  so  that  the  weight  value 

«  th 

needed  for  y^  to  cancel  the  signal  y^  from  the  n  auxiliary  is  the 

weight  value  desired  applied  to  yQ  before  the  correlators  of  the  complex 

reflectivity.  This  weight  value  is  a*.  Note  the  combiner  after  the 

weights,  w  and  w  ,  invert  the  signal  from  w  before  adding.  This  means  the 

signal  out  of  the  combiner  is  w  (yQ  +  y^)  where  w  =  w^  +  jw^.  In  order 

for  the  signal  ay  to  be  cancelled  in  the  combiner,  w*  must  equal  a  . 

n  P  *  n 

The  adaptive  control  therefore  drives  w  =  a^.  Two  types  of  tunable  filters 
are  used  in  the  design.  These  filters  are  needed  to  keep  undesired  signals 
out  of  the  correlators.  TF1  rejects  y  ,  and  TF2  rejects  yR  or  y^.  Note 
quadrature  hybrids  are  generally  considerably  more  expensive  than  power 
splitters  and  an  effort  has  been  made  to  use  as  few  as  possible. 

The  second  method  of  implementing  gradient  control  uses  a 
computer.  It  is  assumed  that  all  the  coupling  parameters  of  the  array  have 
been  measured  over  all  the  desired  frequencies.  From  these  parameters  and 
the  known  or  measurable  reflectivity  one  can  digitally  confute  the  vector 
Q  ^m^.  This  is  the  vector  of  the  needed  for  gradient  control  (with  p  **  0)  . 
One  can  then  use  a  programmable  complex  weight  to  either  multiply  y^  by 
an  or  to  multiply  yQ  by  each  a^*  as  in  Figure  15.  This  method  would  work 
particularly  well  when  the  reflectivities  are  digitally  controlled  as  well. 

The  are  functions  of  the  reflectivity  values  and  when  these  values 
change  the  must  be  recomputed.  Thus  this  method  is  limited  by  the  speed 
at  which  the  can  be  computed.  As  the  size  of  the  array  increases  the 
computation  time  increases  rapidly  and  this  limitation  becomes  more 
significant.  A  second  limitation  is  the  accuracy  of  the  calculation.  Small 
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Figure  22  Antenna 
Pattern  After  Nulling 
an  Interference  Source 


Figure  23  Antenna 
Pattern  After  Manual 
Adaptation  to  a  Desired 
Signal  S  and  Two  Inter- 
ferers  I.  and  I_ 


a 


90° 


The  second  technique  by  which  the  reactance  values  were  adjusted 
used  a  desired  signal,  in  whose  azimuthal  direction  a  main  lobe  was  required, 
and  one  or  more  interference  signals,  in  whose  direction  a  pattern  null  was 
required.  The  desired  signal  and  the  interference  signals  were  slightly 
offset  in  frequency  (M.  kHz)  ,  so  that  they  all  fell  within  the  16  kHz  band¬ 
width  of  the  receiver,  but  were  separated  sufficiently  that  a  subsequent 
narrow  band  filter  (10  Hz  bandwidth)  could  be  tuned  only  to  the  desired  signal. 
By  using  a  hard  limiter  between  the  receiver  output  and  the  narrow  band 
filter,  the  output  of  the  narrow  band  filter  was  a  voltage  proportional  to 
the  output  signal-to-interference  ratio  (SIR) ,  when  SIR  «.  1.  This  voltage 
was  monitored,  and  the  reactances  were  manually  varied  to  maximize  it. 

Two  examples  of  patterns  produced  in  this  manner  are  Bhown  in  Figures 

23  and  24.  In  Figure  23  a  desired  signal  was  incident  on  the  array  at  0  degrees 
and  two  strong  interferers  were  incident  at  180  degrees  and  252  degrees. 

After  the  adjustment  of  the  reactances,  1^  and  I  were  decreased  28  dB  and 
26  dB,  respectively,  relative  to  the  desired  signal.  The  pattern  in  Figure 

24  was  generated  during  a  series  of  measurements  to  investigate  the  minimum 
angular  separation  necessary  to  resolve  the  desired  signal  from  the  interference. 
The  interference  was  separated  20  degrees  from  the  desired  signal  and  was 
suppressed  by  24  dB  relative  to  the  signal. 

The  array  and  the  variables  terminations  were  also  interfaced  with 
a  PDP-11  computer  to  allow  automatic  control  by  the  computer.  The  control 
algorithm  was  a  univariate  search  (see  Section  5.4)  .  The  automatically 
controlled  array  formed  nulls  like  those  with  manual  control.  A  tracking 
program  was  also  written  and  implemented  with  the  computer.  This  allowed 
the  array  to  maintain  cancellation  of  one  or  more  jammes  as  the  platform  was 
slowly  rotated,  representing  a  changing  signal  environment.  Results  of  the 
computer  controlled  array  are  described  in  [4,5]. 


6-7 


7.0  DEVELOPMENT  BY  ZA  OF  ACTIVE  CONTROLLABLE  COMPLEX  TERMINATIONS 
7.1  Design 

In  designing  the  active  variable  complex  termination,  two 
ground  rules  were  established.  These  ground  rules  were  established  to 
make  interfacing  with  the  current  array  easier,  make  any  necessary  debugging 
and  testing  of  the  terminations  or  array  easier,  and  to  increase  the  flexi¬ 
bility  in  array  configurations.  These  ground  rules  were:  1)  that  the 
terminations  should  use  the  existing  variable  reactances  as  variable  phase 
shifters  and  2)  that  the  "loop"  of  the  reflectivity  be  easily  opened.  A 
block  diagram  in  Figure  25  illustrates  the  design  used.  The  termination  was 
assembled  in  three  separate  small  aluminum  boxes  in  order  to  meet  the  ground 
rules  established. 

A  signal  is  captured  by  the  antenna  and  is  transferred  to  the  load 
by  the  matching  network.  The  matching  network  need  not  be  ideal  over  all 
frequencies  but  it  should  be  good  enough  to  prevent  the  reverse  signal  of 
the  mismatch  from  dominating  the  reflectivity.  After  the  matching  network 
is  a  two-way  splitter-combiner.  This  device  is  being  used  as  a  three-port 
circulator  (even  though  it  has  3  dB  loss  each  way)  .  The  forward  flowing 
signal  enters  both  the  upper  and  lower  paths  (in  Figure  25) ,  but  is  blocked 
by  the  isolation  amplifier  in  the  lower  path.  In  the  upper  path,  the  signal 
receives  some  amplification  and  then  enters  another  two-way  splitter-combiner 
again  acting  as  a  3-port  circulator.  Coming  out  of  the  splitter-combiner  the 
signal  is  reflected  off  the  variable  reactance  to  vary  the  phase  for  phase 
control.  The  reverse  flowing  signal  reenters  the  splitter-combiner  and 
enters  both  the  upper  and  lower  paths.  In  the  upper  path  it  is  stopped 
by  the  isolation  amplifier.  In  the  lower  path  it  passes  through  a  bipolar 
amplitude  modulator  for  amplitude  control  and  then  is  amplified.  The  signal 
is  then  narrowband  filtered  and  returned  to  the  first  splitter-combiner, 
which  sends  it  out  the  antenna  again  through  the  matching  network. 

Ther  termination  comprises  three  boxes:  the  antenna  interface  box, 
the  amplitude  control  box  and  the  phase  control  box.  A  complete  schematic 
of  the  antenna  interface  box  and  the  amplitude  control  box  is  given  in 
Figure  26.  The  antenna  interface  box  is  simply  a  four  to  one  broadband 
transformer  (16  to  1  impedance  convertor)  and  a  two-way  splitter-combiner. 
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Block  Diagram  of  Active  Termination 


The  transformer  attempts  to  match  the  real  part  of  the  antenna  impedance  to 
50ft.  The  splitter-combiner  is  included  in  the  antenna  interface  box  so  that 
the  "loop"  of  reflectivity  may  be  opened. 'for  analysis  and  testing.  A  2  d£ 
attenuator  is  included  between  the  transformer  and  the  splitter-combiner  to 
improve  the  impedance  match  the  splitter- combiner  sees.  The  power  loss  is 
not  a  problem  because  of  the  gain  within  the  termination. 

Upon  exiting  the  splitter-combiner  the  signal  also  exits  the  antenna 
interface  box.  Before  entering  the  amplitude  control  box  the  signal  is 
passed  through  a  manually  variable  attenuator  that  varies  from  0  to  25  dB 
and  is  mounted  on  the  side  of  the  amplitude  control  box.  The  attenuator  is 
bidirectional  so  either  port  can  be  used  for  input.  The  attenuator  will 
account  for  variations  in  the  strength  of  the  antenna  mutual  coupling  in 
different  array  configurations.  For  a  given  array  configuration  the 
attenuator  is  set  so  the  array  is  usually  stable.  (Remember  in  Section 
4.2  it  was  shown  that  an  unstable  array  can  often  be  stabilized  if  an  extra 
degree  of  freedom  is  available  in  the  control  space.)  For  a  given  con¬ 
figuration  this  attenuation  need  not  be  changed.  In  practice  the  correct 
value  of  attenuation  would  be  built  into  the  termination. 

The  signal  now  enters  the  amplitude  control  box.  It  first  passes 
through  a  32  MHz  lowpass  filter  to  reject  out-of-band  signals.  Then  the 
signal  passes  through  an  MWA-130  amplifier  for  both  gain  and  reverse  signal 
isolation.  After  the  amplifier  the  signal  is  sent  to  the  phase  control  box 
via  a  splitter-combiner  and  a  3  dB  pad  to  improve  the  impedance  match 
seen  by  the  splitter-combiner. 

The  signal  is  reflected  off  the  variable  reactance  and  it  reenters 
the  amplitude  control  box,  passing  again  through  the  3  dB  pad  and  into 
the  splitter-combiner.  The  signal  is  then  amplified  by  another  MWA-130 
amplifier.  After  the  amplifier  the  signal  passes  through  a  2  dB  pad  and 
then  into  a  Mini-Circuits  PAS-3  bipolar  variable  attenuator. 

The  D.C.  control  for  the  variable  attenuator  comes  from  the 
computer  and  should  be  sent  to  the  amplitude  control  box  through  a  twisted 
pair  of  wires  (or  like  means)  to  reduce  the  effect  of  RF  pickup.  The  control 
signal  is  received  differentially  inside  the  amplitude  control  box  and  then 
sent  through  a  linearization  circuit  to  account  for  the  nonlinear  diode 
effects  of  the  PAS-3.  This  circuit  improves  the  linearity  between  the 
RF  amplitude  out  of  the  PAS-3  to  the  applied  D.C.  control  voltage. 


The  output  of  the  linearization  circuit  is  sent  to  the  PAS-3  through 
a  relay.  The  relay  is  controlled  externally  to  enable  or  disable  the 
termination.  If  +15v  is  connected  to  the  relay  control  line,  the  relay 
will  energize  and  allow  the  control  of  the  variable  attenuator  to  pass 
through.  If  the  relay  control  line  is  open  circuited  or  grounded  then  the 
output  of  the  linearization  circuit  is  open  circuited  and  the  control  line 
to  the  attenuator  is  grounded.  This  causes  the  attenuator  to  apply 
maximum  attenuation  so  virtually  no  signal  passes  through.  This  disables 
both  amplitude  and  phase  control  becaues  the  phase  control  preceeds  the 
attenuator . 

Following  the  attenuator  is  a  final  MWA-130  amplifier,  which  is 
followed  by  a  tunable  narrowband  filter.  The  filter  is  composed  of  a 
series  inductor  and  variable  capacitor.  With  the  current  component  values 
the  center  frequency  of  the  narrowband  filter  can  be  varied  from  12  to 
23  MHz  with  a  constant  3  dB  bandwidth  of  about  900  KHz. 

After  the  narrowband  filter  the  signal  is  returned  to  the  antenna 
through  the  first  splitter-combiner  and  the  matching  network. 

Let  us  now  conduct  a  gain  analysis  of  the  reflectivity  loop.  The 
termination  was  designed  so  that,  with  minimum  attenuation  in  both  the 
manually  variable  attenuator  and  the  electronically  variable  bipolar 
attenuator,  the  gain  of  the  signal  returned  to  the  antenna  would  be  10  dB. 
Figure  27  details  the  gains  and  losses  of  the  termination  in  this  situation. 
Note  to  get  the  desired  amount  of  gain  two  MWA-130  amplifiers  were  not 
sufficient  so  three  were  used.  This  necessitated  the  distribution  of  extra 
attenuation  throughout  the  termination.  This  distributed  attenuation  was 
also  desired  to  improve  the  impedance  match  between  the  various  components. 

Finally  one  must  consider  the  problem  of  added  noise.  The 
added  noise  comes  mainly  from  the  MWA-130  amplifiers.  These  amplifiers 
have  a  noise  figure  of  7  dB.  Consider  first  the  configuration  discussed 
above,  i.e.  that  of  minimum  attenuation.  The  termination  looks  like  a 
a  gain  of  10  amplifier.  Using  the  7  dB  noise  figure  for  the  amplifiers,  one 
finds  the  termination  (with  10  dB  gain)  has  a  15  dB  noise  figure  overall. 

This  noise  figure  for  the  termination  is  not  constant  over  all  control 
settings.  It  is  probably  more  instructive  to  look  at  the  noise  power  inde¬ 
pendent  of  the  signal.  Note  the  output  noise  power  will  be  independent  of 
the  manually  variable  attenuator.  The  output  noise  power  will  vary  with 
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the  electronic  bipolar  attenuator.  The  maximum  output  noise  will  be  when 
the  electronic  attenuator  has  minimum  attenuation,  for  which  the  output 
noise  is  25  dB  above  thermal  noise.  The  minimum  output  noise  will  be  when 
the  electronic  attenuator  has  maximum  attenuation.  Then  the  output  noise 
comes  totally  from  the  last  amplifier  and  appears  at  the  output  13  dB 
above  thermal. 

7.2  Laboratory  Experimental  Data 

In  the  laboratory  at  Z-A  three  experiments  were  conducted  on  the 
amplitude  control  box  of  the  active  termination.  These  experiments 
determined  dynamic  range,  verified  simple  phase  control,  and  determined  the 
bandwidth  of  the  amplitude  control. 

The  setup  for  the  dynamic  range  experiment  is  given  in  Figure  28. 

The  control  voltage  was  varied  for  minimum  and  maximum  attenuation  of  the 
RF  signal.  Minimum  attenuation  occurs  when  vc  =  Ov  or  +10v  (in  the  picture 
vc  ■»  +10v) .  Maximum  attenuation  occurs  when  v^  -  +5v.  The  pictures  of 
Figure  29a  and  b  show  a  dynamic  range  or  83  dB. 

The  phase  control  experiment  is  diagramed  in  Figure  30.  In  the 
experiment  the  reactance  port  was  alternately  open  then  short  circuited. 

The  photos  in  Figure  31a  and  b  show  the  results.  We  see  that  the  output 
does  indeed  shift  by  180°  as  expected. 

The  experiment  to  determine  the  bandwidth  of  the  amplitude  control 
is  diagramed  in  Figure  32.  Four  tests  were  run.  Control  signals  of  sine 
waves  and  square  waves  were  used,  each  at  1  KHz  and  40  KHz.  The  results  of 
these  tests  are  given  in  the  photos  of  Figures  33  and  34.  In  Figure  33a 
a  control  signal  of  1  KHz  sine  wave  is  used  and  in  Figure  33b  a  control  signal 
of  1  KHz  square  wave  is  used.  In  each  of  the  photos  the  low  frequency  control 
signal  is  superimposed  on  the  resultant  RF  output  signal. 

The  feature  to  note  in  Figure  33a  is  that  the  zero  crossings  of  the 
control  signal  and  of  the  RF  modulation  occur  together.  Thus,  at  a  control 
frequency  of  1  KHz,  group  delay  is  not  significant.  In  Figure  33b,  the 
important  feature  is  the  nonlinear  relationship  shown  between  the  control 
voltage  and  the  resultant  RF  amplitude  modulation,  in  particular  the  soft 
limiting  nature  of  the  relationship. 


a)  Minimum  Attenuation 


b)  Maximum  Attenuation 
(note  scale  change) 


Figure  29  Dynamic  Range  Test  Results 


a)  Reactance  Port  Open  Circuit 


b)  Reactance  Port  Short  Circuit 

Top  Trace  =  Output 
Bottom  Trace  =  Input 


Figure  31  Phase  Control  Test  Results 
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Figure  32  Amplitude  Control  Bandwidth  Test  Setup 


Sine  Wave 


b)  Square  Wave 


Figure  33  Amplitude  Control  Bandwidth  Test  Results 
wivl  I  KHz  Control  Signal 
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a)  Sine  Wave 


b}  Square  Wave 


Figure  34  Amplitude  Control  Bandwidth  Test  Results  with  40  KHz 

Control  Signal  (note  shift  in  time  scale  from  Figure  33) 
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In  Figure  34a  one  notices  that  the  sine  wave  control  signal  and 
the  resultant  RF  amplitude  modulation  are  45*  apart.  In  particular,  the 
modulation  lags  the  control  by  45°.  Thus  40  KHz  must  be  the  3  dB  band¬ 
width  of  the  amplitude  control  system. 

7-3  Experimental  Results  of  Parasitic  Array  with  Active  Complex  Terminations 

Zeger-Abrams  built  for  NRL  four  active  terminations  as  described 
in  the  previous  two  sections.  The  first  termination  was  delivered  to  NRL 
in  June  1982.  At  this  time  it  was  demonstrated  that  in  a  small ,. closely 
coupled,  two  element  array,  a  single  active  termination  on  one  of  the  elements 
can  effectively  form  a  null  in  a  desired  direction.  Figure  35  shows  the 
adapted  pattern  in  one  such  case. 

The  remaining  three  active  terminations  were  delivered  to  NRL 
later  in  1982.  Work  at  NRL  has  been  ongoing  to  fully  determine  the 
capabilities  of  active  terminations  in  arrays  with  more  than  2  elements. 


8.0  DEVELOPMENT  OF  A  CONTROL  ALGORITHM  FOR  ACTIVE  TERMINATIONS 

An  algorithm  is  proposed  here  for  control  of  an  actively  terminated 
parasitic  array  to  find  a  stable  null.  This  algorithm  is  a  modification  to 
the  univariate  search  algorithm  currently  used  by  NRL.  The  algorithm 
presumes  the  use  of  oscillation,  or  instability,  detectors  as  described 
in  Section  8.1.  These  detectors  are  minor  modifications  to  the  active 
terminations design  given  in  Chapter  7  and  would  be  incorporated  into  each 
active  termination.  The  algorithm  also  presumes  the  use  of  a  performance 
assessment  indicator  (PAI) ,  possibly  a  signal  to  interference  ratio  detector. 
An  improved  version  of  the  current  detector  used  by  NRL  is  given  in  the 
following  chapter. 

The  algorithm  searches  for  a  stable  null  by  restricting  its 
search  to  that  region  of  the  control  space  for  which  the  array  is  stable. 

The  stability  criterion  of  Chapter  4  (Re  [all  eigenvalues  of  M+P+]  <  1)  is 
difficult  to  apply  in  practice,  because  M+P+  is  hard  to  measure  and  finding 
the  eigenvalues  is  an  involved  computation.  This  is  where  the  oscillation 
detector  comes  of  use.  It  allows  us  to.  find  the  stable  region  of  the 
control  space  without  knowing  M+P+  or  its  eigenvalues. 

8.1  Oscillation  Detector  (Unstable  States) 

The  oscillation  detector  is  based  on  the  following  principle: 
the  final  amplifier  in  at  least  one  of  the  active  terminations  will  saturate 
when  the  array  goes  unstable  (starts  to  oscillate  at  some  frequency)  . 

This  fact,  along  with  the  knowledge  that  the  signals  encountered  in  normal 
operation  (interference  or  desired)  should  not  saturate  this  amplifier, 
allows  us  to  detect  when  an  array  is  unstable.  We  can  detect  an  unstable 
array  by  noting  when  any  of  these  final  amplifiers  saturate. 

Figure  36  shows  a  modification  to  the  design  of  the  active 
termination  of  Figure  26.  This  modification  detects  when  the  peak  amplitude 
of  the  MWA-130  output  RF  crosses  a  certain  threshold.  The  op  amp  is 
configured  as  a  comparator.  The  1.4  Mfi  feedback  resistor  provides  positive 
feedback  for  hysteresis  to  prevent  chatter.  The  op  amp  is  one  of  two 
unused  op  amps  currently  within  the  termination.  Laboratory  experiments 
have  determined  that  when  the  op  amp  threshold  is  set  for  1.4v,  the 


V.* 


■  .0 


8-1 


output  RF  is  +15  dBm.  We  selected  to  trip  the  threshold  with  +15  dBm  output 
of  the  MHA  amplifier,  because  this  is  higher  than  any  expected  operational 
signal  (interference  or  desired)  and  yet  slightly  less  than  the  saturation 
point  of  the  amplifier  (its  1  dB  compression  point  is  +18  dBm) .  With 
the  threshold  set  in  this  manner,  we  can  detect  when  the  array  is  near 
oscillation  but  the  MWA-130  is  still  linear.  The  output  of  the  op  amp  is 
+15v  if  no  oscillation  is  present  and  -15v  if  an  oscillation  is  present. 

If  an  inverted  output  is  desired,  the  final  unused  op  amp  can  be  employed 
to  invert  this  signal.  In  practice,  eight  of  these  signals  would  probably 
be  combined  to  represent  one  I/O  byte  to  the  computer. 

Finally,  if  it  is  found  that  the  detector  will  trigger  on  oper¬ 
ational  signals  (in  particular  strong  interferences)  and  not  just  on 
oscillations,  one  can  obtain  some  extra  operating  signal  range  with  some 
modifications.  First,  one  may  raise  the  threshold  of  the  comparator 
slightly.  Secondly,  the  2  dB  pad  and  the  splitter-combiner  within  the 
antenna  interface  box  can  be  replaced  with  a  10  dB  directional  coupler.  In 
this  case  the  ports  on  the  antenna  interface  box  will  no  longer  be  symmetric. 
The  antenna  is  connected  to  the  IH  port  of  the  directional  coupler 
(through  the  matching  transformer) .  The  COUPLED  port  of  the  coupler  is 
output  to  the  termination.  The  OUT  port  of  the  coupler  is  the  return 
signal  from  the  termination.  In  this  configuration,  one  gains  about  4  dB 
more  range  for  operating  signals.  The  sacrifice  is  a  corresponding  increase 
in  added  noise. 

8.2  Modified  Univariate  Search  Algorithm 

The  simplest  way  to  control  the  active  terminations  will  be  to 
modify  the  current  algorithm,  which  is  a  univariate  search.  The  simplest 
modification  of  all  is  to  conduct  the  univariate  search  as  before;  however, 
consider  only  those  operating  points  for  which  no  termination  is  oscillating. 
More  specifically,  while  ramping  one  of  the  control  variables  (either 
amplitude  or  phase)  and  just  before  evaluating  the  performance  assessment 
indicator  (PAI)  for  each  new  operating  point,  one  checks  to  see  if  any 
of  the  terminations  are  oscillating.  If  am  oscillation  is  present  in  any 
of  the  terminations,  then  this  operating  point  is  rejected  (without  bothering 
to  check  the  PAI)  just  as  if  its  PAI  had  been  insufficient.  Upon  reaching 
a  null,  the  same  process  of  checking  for  oscillations  is  continued  even 
in  tracking  mode. 


The  main  advantage  of  this  algorithm  is  its  ease  of  implemen¬ 
tation.  The  main  disadvantage  is  inherent  to  the  univariate  search 
algorithm,  namely,  it  is  slow.  However,  the  current  algorithm  use  two 
different  step  sizes  for  when  it  varies  the  control  signal.  A  large 
step  size  when  cancellation  is  poor  and  a  small  step  size  for  fine  tuning 
a  null.  This  arrangement  greatly  reduces  the  number  of  operating  points 
evaluated  while  cancellation  is  still  poor  and  greatly  speeds  up  the 
algorithm.  Thus  it  is  felt  that  more  speed  is  not  essential  to  the  current 
experiments  with  parasitic  arrays. 

A  second  possible  problem  with  this  algorithm,  which  also  is  a 
possible  problem  with  any  non-random  algorithm,  is  that  it  may  get  hung 
up  at  a  local  minimum  which  is  not  an  adequate  null.  This  is  only  called  a 
possible  problem  because  it  becomes  unlikely  with  excess  degrees  of  freedom  in 
the  array,  i.e.,  when  there  are  more  actively  terminated  elements  than  there  are 
interferences  to  be  nulled.  It  should  be  noted  that  it  may  be  desirable 
to  have  more  than  one  excess  degree  of  freedom.  Excess  degrees  of  freedoms 
may  also  be  used  to  stabilize  an  otherwise  unstable  array.  Finally,  if 
the  array  still  forms  an  inadequate  null,  even  with  the  extra  degrees  of 
freedom,  the  current  algorithm  will  recognize  this  and  start  over  from  a 
new  location  which  hopefully  will  lead  to  a  suitable  null. 
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9.0  DEVELOPMENT  OF  AM  IMPROVED  PERFORMANCE*  ASSESSMENT  INDICATOR  (PAI) 

Every  control  algorithm  needs  some  method  to  evaluate  the  current 
utility  of  the  system  being  controlled.  The  device  which  does  the  evaluating 
is  the  performance  assessment  indicator  (PAI) .  For  a  conventional  adaptive 
array  with  LMS  control  the  PAI  is  simply  the  array  error  signal  which  one 
then  tries  to  minimize.  For  an  SIC  the  usual  PAI  is  the  signal  to  inter¬ 
ference  ratio  at  the  array  output  (to  be  maximized) . 

In  designing  a  PAI  for  the  current  experimental  array,  working 
at  an  IF  frequency  would  greatly  ease  the  complexity  of  the  circuit.  Since 
the  current  PAI  for  the  array  uses  a  Watkins-Johnson  WJ8718  receiver,  it 
was  felt  that  the  same  receiver  should  be  incorporated  into  the  new  design. 
Figure  37  hhows  a  basic  design  of  a  PAI  using  the  WJ8718  receiver.  The 
main  antenna  is  connected  to  receiver  through  an  antenna  matching  network 
(possibly  as  simple  as  a  transformer) .  The  receiver  is  set  up  to  convert 
the  signal  to  a  455  KHz  IF  with  16  KHz  bandwidth.  At  an  IF,  the  operation 
of  the  signal-interference  discriminator  is  greatly  simplified. 

One  may  ask:  since  the  discriminator  separates  the  desired 
signal  from  the  interference,  why  have  an  adaptive  array  at  all?  The  problem 
is  that  the  discriminator  is  not  a  perfect  separator.  The  discriminator 
simply  uses  some  feature  of  the  desired  signal  which  distinguishes  it  from 
tne  interference  and  then  produces  two  outputs  such  that  one  of  the  outputs 
looks  more  like  the  signal  and  less  like  the  interference  than  the  other. 

The  amount  of  discrimination  is  far  from  enough  to  allow  one  to  attach  a 
receiver  to  the  desired  signal  port  but  is  enough  to  control  an  adaptive 
array. 

The  two  signals  out  of  the  discriminator  are  each  sent  through 
detectors,  either  power  or  peak  detectors.  The  outputs  of  these  detectors 
are  D.C.  signals  which  are  then  each  input  to  a  log  ratio  amp.  A  log  ratio 
amp  outputs  a  voltage  proportional  to  the  logarithm  of  the  ratio  of  the  two 
inputs.  In  this  case  the  output  is  an  attempt  at  measuring  SIR  in  d£. 

The  key  element  of  this  design  is  the  discriminator.  Three  types 
of  signal  discriminators  will  be  introduced.  The  recommended  type  will 
be  described  in  detail.  Tne  three  different  signal-interference  discriminants 
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are:  1)  the  signal  and  jammer  are  narrowband  (or  even  CW)  and  slightly 
separated  in  frequency,  2)  the  signal  is  modulated  with  a  psuedo-noise  (PN) 
sequence  and  the  jammer  can  be  arbitrary,  and  3)  the  signal  is  narrowband 
(doublesided  bandwidth  6  KHz)  and  the  jammer  is  wideband. 

The  jammer-signal  environment  required  for  the  type  1)  discrim¬ 
inator  is  really  a  fairly  benign  environment.  It  is  not  likely  to  be  the 
jammer-signal  environment  encountered  in  practice.  This  discriminator 
would  be  used  only  with  experimental  arrays  to  demonstrate  nulling  capabilities. 
This  discriminator  would  typically  be  composed  of  two  narrowband  filters; 
however  with  a  single  notch  filter  (at  the  desired  signal  frequency)  one 
could  use  the  same  technique  as  described  below  for  the  type  3)  discriminator. 

Type  2)  discriminator  has  the  advantage  of  working  against 
any  jammer  that  is  uncorrelated  with  the  PN  sequence.  The  total  signal 
going  into  the  discriminator  is  split  into  two  paths.  Both  paths  are 
multiplied  by  the  PN  sequence.  One  path  (the  signal  path)  then  passes 
through  a  narrowband  filter  at  the  IF  frequency.  The  other  path  (the 
interference  path)  is  passed  through  a  notch  filter  at  the  IF  frequency. 

e  outputs  of  these  filters  are  then  sent  to  detectors  as  in  Figure  37. 

The  disadvantages  with  the  PN  discriminator  are  the  increased  hardware 
complexity  and  the  need  to  synchronize  to  the  PN  sequence  of  the  desired 
signal  (which  will  have  come  from  a  distant  transmitter  in  practice) . 

The  type  3)  discriminator  is  the  type  recommended.  It  requires 
little  more  hardware  than  the  type  1)  and  may  be  employed  against  a  far 
greater  variety  of  jammers.  In  fact,  it  can  be  used  anytime  type  1)  can 
be  used,  as  long  as  the  frequency  separation  is  wider  than  3  KHz.  The 
type  3)  discriminator  assumes  that  the  jammer  has  wider  bandwidth  than  the 
signal,  or  at  least  has  a  good  fraction  of  its  power  out  of  this  band.  The 
input  signal  (S  +  I)  is  split  and  one  part  is  sent  through  a  notch  filter 
at  the  IF  frequency.  The  power  remaining  in  this  signal  will  be  proportional 
to  the  interference  power.  We  now  have  two  signals:  one  is  S  +  I  and 
the  other  representative  of  I.  We  would' like  to  have  two  signals  such  that 
one  is  representative  of  I  and  the  other  of  S  alone.  If  we  use  the  signals 
I  and  S  +  I  in  the  detectors  of  Figure  37,  then  the  output  of  the  log  ratio 
amplifier  will  be  log[(S  +  I)/l].  This  output  cannot  be  used  for  a  PAI 
because  if  I  »  S  the  output  is  =  0.  The  PAI  will  change  little  even  for  large 
changes  in  I . 


A. 14  If  y  is  a  scalar  function  of  n  variables,  x4 ,  arranged  into  a  column 

X 

(row)  vector  x  (x  ),  the  expression 

3y/3x  (3y/9xT) 

denotes  a  column  (row)  vector  with  elements  3y/3x^. 

A. 15  If  £  is  a  column  vector  with  m  elements,  each  a  function  of  n 

T  T 

variables,  x^,  arranged  into  a  row  vector  x  ,  the  expression  3y/3x 
denotes  a  matrix  with  m  rows  and  n  columns,  with  elements  dy^/dxy 


A. 16  3UV/3x  -  (3U/3x)V  +  U(3V/3x) 


A.  17 


3Y/3y^j  -  E±J  “  e^e^  where  E±j  is  a  matrix 
for  the  i,  j-th  element  which  =  1;  and  e^  is 
except  for  the  i-th  element  which  ■  1. 


of  all  zeroes,  except 
a  vector  of  all  zeroes 


A.  18  3aTx/3x  -  a  and  d^x/dx?  «  aT  (if  the  elements  of  a  are  not 
functions  of  any  x^) . 

A. 19  3Ax/3xT  “  A 

A.  20  Chain  Rule:  3y_/3xT  =  ( 3^/ 3^T)  (3^/3xT) 

A.  21  3Y_1/3x  -  -y'^SY/SxIY-1 


-1 


(B  -  CE_1D)_1 

•  -1  -1  -1 
-B  C(E  -  DB  C) 

-E_1D(B  -  CE'V1 

(E  -  DB^C)”1 

.  A 


This  can  be  verified  by  multiplying  A  and  A  ^  and  showing  that  the 
result  reduces  to  the  identity. 


A.  10  Determinant  of  a  partitioned  Matrix  given: 


and 


then  | A 1  =  | E | • | (B  —  CE  1D) j  if  E  1  exists, 

and  =  | B | • | (E  -  DB  ^C) |  if  B  1  exists. 


MATRIX  DIFFERENTIATION  FORMULAE 


A. 11 


If  the  elements  of  a  matrix  Y(m  x  n)  are  functions  of  a  scalar,  x, 
the  expression 

3Y/3x 


denotes  a  matrix  of  order  (m  x  n)  with  elements  3y^/3x. 

T 

A.  12  If  the  elements  of  a  column  (row)  vector  £  (x  )  are  functions  of  a 
scalar,  x,  the  expression 

32./3x  Ol?/3x) 

denotes  a  column  (row)  vector  with  elements  3y^/3x. 

A.  13  If  y  is  a  scalar  function  of  m  x  n  variables,  x^,  arranged  into  a 
matrix,  X,  the  expression 

3y/3X 


denotes  a  matrix  with  elements 


dy/Sx^ . 
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A. 4  If  (A  +  B)  *  exists,  then  ABT  is  symmetric  (i.e.,  AB^  -  (ABT)T) 


if  and  only  if  (A  +  B) 
proof  of  "if": 


-1 


(A  -  B)  is  symmetric. 


(A  +  B)  - 1  (A  -  B) 
(A  -  B) (AT  +  BT) 

aat  -  bat  +  abt  -  bbt 
abt 


[<A  +  B)_1(A  -  B)]T  ■ 

T  T 

(A  +  B)(AA  -  B1) 

T  T  T  T 

AA  +  BA  -  AB  -  BB 


-  BA 


T  T 
-  (AB1)1 


(A1  -  bW  +  B1)' 


proof  of  "only  if":  Reverse  the  above  argument. 

-1  T 

A. 5  If  (A  +  B)  exists,  then  A  B  is  symmetric  if  and  only  if 
(A  -  B)(A  +  B)  *  is  symmetric. 

This  can  be  proved  with  minor  variations  to  the  proof  of  (A.  4). 


A. 6  If  (A  +  B)  *  exists,  then: 

(A  -  B) (A  +  B)-1 

2B(A  +  B)_1 


(A  +  B  +  2B) (A  +  B)  1  -  I  -  2B(A  +  B) 

I  -  (A  -  B)(A  +  B)-1 


A. 7  If  (A  +  B) 


-1 


exists,  then: 

-1 


(A  -  B)  (A  +  B) 


or 


2A(A  +  B) 


-1 


[2A  -  (A  +  B)  ] (A  +  B)*"1  -  2A(A  +  B)"1  -  I 

I  +  (A  -  B)  (A  +  B)-1 


The  following  two  identities  and  the  matrix  differentiation 
formulae  are  taken  from  CRC,  Standard  Mathematical  Tables. 

A. 8  Given  A  ■  B  +  UV,  where  A  and  B  are  nxn,  Uisnxp,  Vispxn, 
and  B  exist,  then: 

A-1  -  B-1  -  B_1U(I  +  VB~1U)“1VB_1 

This  can  be  verified  by  multiplying  A  and  A  !  and  showing  that 
the  result  reduces  to  the  identity. 


A. 9  Formula  for  inverting  a  partitioned  matrix.  Given: 

/  B  ;  G  \ 

A  *  I . I  ,  where  A  is  (p  +  q)  x  (p  +  q) ,  B  is  p  x  p, 

\":  V 

'  -1  _i 

,  D  is  q  x  p,  Eisqxq,  and  B  and  E  exist,  then: 


C  is  p  x  q 


APPENDIX  A 


MATRIX  IDENTITIES  AND  DIFFERENTIATION  FORMULAE 

Identities  (A.l)  through  (A.  7)  are  specialized  identities 
found  to  be  useful  in  the  text.  They  are  given  with  proof.  Identities 
(A.  8)  and  (A. 9)  and  the  differentiation  formulae  are  taken  from  CRC. 
Standard  Mathematical  Tables.  21st  Edition  [12],  pp.  125-136. 

MATRIX  IDENTITIES 


In  Identities  (A.l)  through  (A. 8),  A  and  B  are  square 


matrices . 


If  A  *  exists. 

then  A  *B  ■ 

ba" 

*  if  and  only  if  AB 

proof  of  "if": 

AB 

BA 

A_1 ABA_1 

- 

A_1 BAA_1 

BA-1 

- 

A_1B 

proof  of  "only 

if": 

A_1B 

BA_1 

AA_1BA 

m 

ABA”1 A 

BA 

m 

AB 

BA. 


A. 2 


(A  -  B)  (A  +  B)  -  (A  +  B)  (A 
proof  of  "only  if" : 

(A  -  B)  (A  +  B) 

2  2 
A  -  BA  +  AB  -  B 

-BA  +  AB 

AB 


-  B)  if  and  only  if  AB 

-  (A  +  B) (A  -  B) 

-  A2  +  BA  -  AB  -  B2 

-  BA  -  AB 

-  BA 


BA. 


proof  of  "if":  Reverse  the  above  argument. 

A. 3  If  (A  +  B)  1  exists,  then  (A  +  B)  *(A  -  B)  -  (A  -  B) (A  +  B) 
if  and  only  if  AB  ■  BA. 
proof:  Apply  (A.l)  to  (A. 2). 
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Holt,  Rinehart,  and 


jammer  is  of  a  wider  bandwidth  than  the  signal,  or  at  least  contains 
significant  energy  outside  the  signal  bandwidth  of  ±3  KHz  of  the  center 
frequency.  This  also  works  as  the  same  discriminant  that  NHL  has 
currently  been  using  (namely  frequency  separation)  although  at  a  somewhat 
greater  frequency  difference  (3  KHz  instead  of  1  KHz) .  The  final  statistic 
that  drives  the  array  is  basically  log[(S  +  I)/I^]  where  y  >  1.  The 
optimum  value  of  Y  is  probably  between  1.5  and  2. 

Zeger-Abrams  developed  a  new  model  for  analysis  of  a  parasitic 
array,  particularly  useful  for  analysis  of  active  terminations  and  adaptive 
control.  Zeger-Abrams  developed  active  terminations  to  be  used  with  a 
compact  parasitic  array  at  HF  and  developed  circuitry  and  an  algorithm  for 
their  control.  Finally,  Zeger-Abrams  developed  am  improved  PAI  circuit 
to  drive  an  adaptive  parasitic  array.  Possible  areas  for  future  development 
and  improvement  include:  1)  improvements  in  array  stability  through 
improvements  in  the  design  of  the  elements,  the  terminations,  and  the 
algorithm,  2)  analytic  studies  of  the  array  performance  that  include 
noise  and  studies  of  combining  the  element  outputs  in  the  manner  of  a 
conventional  array,  3)  further  improvements  in  design  of  a  PAI,  and 
4)  the  study  and  design  of  improved  control  algorithms. 


10.0  SUMMARY  AND  CONCLUSIONS 


The  parasitic  array  is  a  promising  method  of  building  an  adaptive 
array  in  a  small  package.  This  is  a  particularly  appealing  prospect  at 
HF,  where  a  conventional  array  would  be  measured  in  100's  of  feet  or 
even  100's  yards. 

The  theory  of  an  adaptive  parasitic  array  was  presented  with 
emphasis  on  active  terminations.  Two  theoretical  models  of  a  parasitic 
array  were  discussed,  the  Impedance  model  and  the  transmission  line 
model.  The  array  was  described  in  terms  of  both  models  and  then  the 
relationship  between  the  models  was  established,  should  one  wish  to 
transform  from  one  model  to  the  other.  With  this  relationship  established, 
we  used  the  transmission  line  model  throughout  the  remainder  of  the  report. 

An  explicit  solution  for  the  reflectivity  values  (of  the 
transmission  line  model)  necessary  to  form  nulls  on  jammers  was  derived. 

The  issue  of  stability  for  an  actively  terminated  array  was  discussed. 

The  concept  of  degrees  of  freedom  for  the  array  was  introduced.  The 
possibility  of  using  excess  degrees  of  freedom  to  stabilize  the  array  was 
demonstrated.  Several  ideas  on  how  to  make  the  array  less  likely  to  go 
unstable  where  suggested. 

Experimental  results  obtained  in  adaptive  nulling  with  a  passive 
reactive  parasitic  array  was  reviewed.  These  results  were  obtained  at 
HF  on  the  NKL  Brandywine  Antenna  Range  with  an  array  about  one  meter  high 
and  a  half  meter  in  diameter. 

Zeger-Abrams  designed  and  built  four  active  terminations  for 
the  NRL  array.  Preliminary  results  of  using  these  terminations  in  the 
parasitic  array  were  described. 

As  noted  above,  an  actively  terminated  array  may  have  stability 
problems.  Zeger-Abrams  developed  an  instability  detector  (actually  a 
near- instability  detector)  to  be  used  in  conjunction  with  a  control 
algorithm.  A  simple  modification  to  use  this  detector  with  the  current 
NRL  computer  control  program  was  described. 

Finally,  in  order  to  improve  the  performance  of  the  NRL  array, 
Zeger-Abrams  designed  an  improved  performance  assessment  indicator  (PAI) 
that  drives  the  adaptive  control  algorithm.  This  PAI  assumes  that  the 
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the  output  of  the  difference  amp  is  the  PAI  statistic  negated.  If  V 

O  T  i 

is  the  output  voltage  of  the  upper  log  amp  and  Vx  is  the  output  voltage 
of  the  lower  log  amp,  then  the  true  output  of  the  op  amp  difference  ampli¬ 
fier  is  2(YV1  -  Vg  +  j) .  The  gain  of  the  lower  (non-inverting)  input 
is  determined  by  the  10K  potentiometer.  If  the  pot  is  set  so  there  is  a 
short  circuit  between  the  log  amp  output  and  the  non-inverting  terminal, 
then  the  gain  for  the  non-inverting  signal  is  twice  that  of  the  inverting 
signal,  y  *  If  the  pot  is  set  so  the  non- inverting  terminal  is  shorted 
to  ground,  then  (quite  obviously)  the  gain  for  the  non-inverting  signal 
is  zero,  y  ■  0.  Finally,  if  it  is  desired  that  the  PAI  be  positive  going 
for  increasing  S/I  then  one  of  the  excess  op  amps  can  be  used  to  invert  it. 

Several  possible  performance  assessment  indicators  have  been 
discussed.  The  recommended  PAI  has  been  described  in  detail.  A  schematic 
of  the  recommended  design  was  presented  in  Figure  39. 


Table  1 

Notch  Filter  Inductors  Specifications 


Inductor 


Approximate  Number 
of  Windings 


Wire  Gauge 


6.41  mH 

160 

32 

5.71  mH 

150 

32 

3.0  mH 

110 

32 

1.85  mH 

88 

32 

1.76  mH 

75 

32 

45  yH 

14 

26 

43.6  yH 

I 

13 

26 

All  cores  ai 

re  Ferroxcube  pot  cores 

1408C-A250-3B7, 

which  requii 
set  each. 

re  the  1408F1D  bobbin  a i 

id  1408H  hardware 

■  _•  1  !■  l"  ■  -■  .■  I"  IP" IF  ' 


r  -  i  ■ — .  * — . 


more.  This  is  clearly  impractical .  To  get  the  selectivity  desired,  it  is 
recommended  that  the  signal  at  455  KHz  be  converted  to  a  lower  IF  —  in  parti¬ 
cular  50  KHz.  This  could  be  accomplished  by  mixing  with  an  L0  at  505  KHz  to 
produce  outputs  at  50  KHz  and  1.06  MHz  (see  Figure  39)  .  The  1.06  MHz  component 
is  then  easily  filtered  out.  The  filtering  is  accomplished  with  a  third  order 
Butterworth  lowpass  filter  with  a  cutoff  at  70  KHz.  This  filter  has  less 
than  1  dB  loss  at  50  KHz  and  ideally  70  dB  loss  at  1.06  MHz.  The  final 
concern  is  the  necessary  stability  of  the  505  KHz  LO.  If  the  LO  is  in¬ 
accurate  then  the  desired  signal  is  not  centered  in  the  notch  and  some  power 
may  leak  through  the  filter.  This  leak  through  power  may  corrupt  the  PAI 
measurement  and  degrade  the  final  achievable  S/I.  100  Hz  accuracy  in  the 
ocsillator  is  felt  sufficient  to  make  this  degradation  negligible.  100  Hz 
accuracy  in  a  505  KHz  ocsillator  corresponds  to  a  frequency  stability  of 
200  PPM  or  0.02%. 

The  notch  filter  is  a  design  previously  used  by  Zeger-Abrams .  It 
had  been  built  and  tested  and  found  to  work  very  well.  Table  1  gives 
specifications  for  the  inductors.  It  is  suggested  that  several  extra 
windings  be  placed  on  each  inductor  and  then  resonate  the  inductor  with  a 
known  capacitance  to  "tweek"  the  inductance  value. 

The  signal  into  the  notch  filter  first  goes  through  a  voltage 
divider.  This  lowers  the  input  impedance  seen  by  the  filter  and  makes  the 
necessary  inductors  smaller.  There  also  is  a  signal  attenuation  associated 
with  this  divider,  but  the  signal- to-noise  ratio  of  the  signal  is  assumed 
to  have  been  well  established  by  the  WJ8718  receiver  and  this  is  not  deemed 
a  problem.  The  output  op  amp  provides  an  infinite  output  impedance  to  the 
filter  and  restores  the  gain  lost  in  the  input  resistor  divider. 

The  output  of  the  notch  filter  as  well  as  a  sample  of  the  input 
signal  are  then  sent  to  separate  peak  detectors.  Since  the  signal  levels 
may  be  smaller  than  a  diode  voltage  drop,  it  was  necessary  to  use  ideal 
peak  detectors.  The  KC  network  at  the  output  of  the  op  amp  sets  the  fall 
time  constant  of  the  peak  detector.  This  time  constant  is  currently  set 
for  1  ms. 

The  outputs  of  the  peak  detectors  go  into  log  amplifiers.  The 
outputs  of  the  amplifiers  are  then  sent  into  a  difference  amplifier  to 
form  the  PAI  statistic  log[(S  +  I)/I^]  ■  log  (S  +  I)  -  ylog(I) .  Note 
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Consider,  however,  the  circuit  diagramed  in  Fiugre  38.  The  output  of 
this  circuit  is  log[(S  +  I)/I^],  or  equivalently  log[(S  +  I)/l]  -  (y  -  l)log(I). 
For  large  I  this  becomes  -(y  -  l)log(I).  Since  y  >  1,  this  term  becomes  increas¬ 
ingly  more  negative  as  I  increases.  Thus  this  circuit  will  function  as  a 
useful  PAI.  There  is  a  compromise  to  be  made  when  using  this  circuit.  When 
1  >  S,  we  would  want  y  =  2  so  that  the  control  is  driven  by  =  -log (I) .  If 
y  were  less,  then  we  would  be  reducing  the  loop  gain  and  increasing  the  time 
it  would  take  to  reduce  a  large  interference.  However  when  S  >  I,  we  would 
want  y  -  1  so  that  the  control  is  driven  by  -  log(S/I) .  For  larger  y,  the 
array  seeks  to  reduce  I  with  a  single-mindedness  that  is  counterproductive. 

It  will  strive  to  reduce  I  even  to  the  detriment  of  S/I.  In  particular,  if 
Y  *  2,  the  array  will  reduce  S/I  by  as  much  as  1  dB,  in  order  to  reduce  I 
by  1  dB.  Note  this  is  better  than  no  desired  signal  maintainance  at  all  as 
one  would  get  by  using  -log(I)  alone  as  a  PAI.  It  is  thought  that  the 
optimum  value  for  y  is  between  1,5  and  2.  Finally,  it  may  be  possible  to 
make  y  itself  adaptive,  although  this  is  not  recommended  in  the  initial 
implementation.  One  such  variable  y  might  be  y  *=  2  -  f£log[(S  +  I)/l]j 
where  f  [x]  =  x  if  x  <  1  and  1  if  x  >  1.  Here  log[(S  +  I) /I]  would  be 
obtained  by  finding  the  difference  of  the  log  amp  outputs  before  multiplying 
by  y. 

The  block  diagram  in  Figure  38  is  basically  the  recommended  PAI 
design.  The  critical  part  of  the  design  is  the  notch  filter.  Analysis 
shows  that  the  maximum  achievable  signal- to-interference  ratio  out  of  the 
array  is  the  ratio  of  the  input  to  output  signal-to-interference  ratios  of 
the  notch  filter,  i.e. 


max  (-)  =  (3/1)  filter  input 

\I/  array  (S/I)  filter  output. 


Assuming  one  wants  a  maximum  S/I  output  of  the  array  to  be  at  least  30  dB, 
then  one  would  want  the  desired  signal  rejection  ratio  of  the  notch  filter 
to  be  at  least  30  dB  without  significantly  reducing  the  nearby  interference. 

A  notch  filter  with  a  6  KHz  3  dB  notch  bandwidth  at  center  frequency 
of  455  KHz  and  sufficiently  steep  skirts  to  satisfy  the  above  requirements 
was  considered.  The  design,  however,  required  inductors  with  Q's  of  800  and 


APPENDIX  B 


RELATING  S-PARAMETERS  TO  Z- PARAMETERS 

For  any  N-port  network  there  are  two  equivalent  representations 
S  (scattering) -parameters  and  Z  (impedance) -parameters.  Scattering  para¬ 
meters  are  defined  in  Equation  (B5)  and  analyze  the  system  through  the 
forward  (incident)  and  reverse  (reflected)  voltage  traveling  waves. 
Impedance  parameters  are  defined  in  Equation  (Bl)  and  analyze  the  system 
through  the  port  currents  and  voltages.  This  appendix  derives  the 
relationship  between  these  two  represations. 

Consider  the  N-port  network  in  Figure  B.l,  terminated  with 
transmission  lines.  Let  Z  =  [z^]  be  the  open-circuit  impedance  matrix 
for  the  network.  The  entries  of  this  matrix  are  defined  by 

i  =  0,  for  all  k  ^  j  .  (Bl) 

We  define  the  vectors  i  and  v 


From  the  definition  of  z  and  superposition  we  see  that 

v  =  z_i  .  (B2) 

The  transmission  lines  connected  to  the  nth  port  has  a  characteristic 

impedances  of  z 

on 

z  = 

o 


B-l 


T 


The  voltage-current  telegraph  equations  for  each  of  the  trans¬ 
mission  lines  can  be  solved  (14,  17,  18]  to  show  that  each  voltage  and 
current  is  the  sum  of  two  components: 


v  —  V  +  V 

Vn  V (+) n  v(-)n 


i  —  i .  +  i  .  . 

n  (+)n  (-)n 


and  in  vector  form 


(B3) 


v 


i 


=  V  +■  V 

-(+)  -{-) 


=  i(+)  +  V) 


(B4) 


Here  v^+jn  ^(+)n  are  forward  voltage  and  current  "waves"  traveling 

toward  the  N-port  network  and  v ,  .  and  i  .  ,  are  the  reverse  waves 

(-)n  (-)n 

traveling  away  from  the  N-port.  The  scattering  parameter  matrix  S  =  [s^] 
is  defined  by 


s.  .  =  (-)  i 

^  - — : 

^ 3  V (+) k  =  0,  for  all  k  ?  j. 
From  this  definition  and  superposition,  we  see  that 


(B5) 


^(-)  =  S^<+) 


(B6) 


It  can  also  be  shown  from  the  telegraph  equations  that  these  traveling 
waves  are  related  by  the  following  equations: 


V  V 

(+)n  (-) n 

i  i  "  zon 

(+) n  (-) n  (B7) 


B-3 


Substituting  Equations  B7  into  B3  yields: 


Vn  =  V(+) n  +  v(-)n 
i  -  V (+) n  _  V(~)n 

Q  2  Z 

on  on  (B8) 


Upoh  solving  Equation  (B8)  for  v(+)n  an&  v(_)n  we 


v .  .  =  >5  (v  +  z  i  ) 

(+)n  n  on  n 


v.  .  =  (v  -  z  i  ) 

(-)n  n  on  n 


and  in  vector  form 

v  =  *3  (v  +  Z  i) 
-(+)  -  o- 

v,  ,  =  h  (v  -  Z  i) 
— (-)  —  o— 


(B9) 


(BIO) 


Equation  (B2)  can  be  solved  for  £  and  substituted  into  Equation  (BIO) . 


v(t)  -■.<!♦  V"1^ 

V>  '  '  (I  -  zo2'l,i 

where  I  is  the  matrix  identity 


(Bll.a) 

(Bll.b) 


Equation  (Bll.a)  can  be  solved  for  v  and  substituted  into  Equation  (Bll.b) 


V)  ■  (I  "  v"1'  u  + 


(B12) 


By  comparing  Equation  (B12)  to  (B6) ,  we  see  that 


s  =  (i  -  z  z_1)  (i  +  z  z'V1 

o  o 

=  (z  -  z  )z-1  z(z  +  z  )-1 

o  o 

=  (Z  -  Z  )  (Z  +  Z J'1  (B13) 


Using  the  matrix  identies  A. 6  and  A. 7  of  Appendix  A,  we  find  the  following 
equivalent  forms  for  S : 

S  =  I  -  2Z  (Z  +  Z  )_1  (B14) 

o  o 

=  2Z  (Z  +  Z  )_1  -  I  (B15) 

o 

Note  if  all  the  z  's  are  the  same,  then  with  the  fact  that  Z 
on 

is  symmetric,  S  must  also  be  symmetric. 

Finally,  Equation  (B13)  can  be  solved  for  Z  in  terms  of  S  and  Z^. 

S  =  (Z  -  ZQ)  (Z  +  ZQ)  1 

S(Z  +  Z  )  =  (Z  -  Z  ) 

O  o 

(I  +  S)Z  =  (i  -  S)z 

O 

z  =  (I  -  S) _1 (I  +  S)Zo  (Bl6) 


APPENDIX  C 


SCATTERING  PARAMETERS  IN  TERMINATED  TRANSMISSION  LINES 

4 

C. 1  Single  Transmission  Line 

Derivations  of  the  theory  used  in  this  appendix  are  given  in 
[14,  18].  To  solve  transmission  line  problems  one  applies  three  rules: 


1)  v (x)  =  v(+)(x)  +  v^_j(x),  i (x)  =  i(+j(x)  +  i  (x)  and 


x(x) 


W° 


‘hw 


are  everywhere  continuous ,  where  x  is 


the  distance  along  the  transmission  line. 


1  +  8  ( X.) 

2)  z(x)  =  z  (x)  *5 - —  and  s(x) 

o  1  -  s(x)  ' 


z(x)  -  z  (x) 
_ o 

z  (x)  +  z  (x) 
o 


where  zq(x)  is  the  characteristic  impedance  at  the  point  x. 


ZQ(x)  - 


X  >  l 
0  <  x  <  l 
x  <  0 


z(x)  is  simply  the  impedance  the  transmission  line  could  be 
terminated  with  at  the  point  x,  and  not  change  the  load  as  seen 
by  the  source. 


3) 


8(X1) 


ej2&(x1  xQ)  8ix^j  where  £  «  u/v^  and  is  the  propa_ 


gation  velocity  of  the  transmission  line.  This  is  just  the 
change  in  phase  due  to  the  change  in  position  along  the  trans¬ 
mission  line.  Note  the  change  in  phase  is  calculated  from 
twice  the  change  in  distance  because  it  is  a  round  trip. 


To  apply  these  rules  to  the  transmission  line  network  in  Figure 
C.l,  we  start  at  the  load,  z  .  Here  z  ( l +)  *  z  and  therefore  by 

A  A 

Rules  1  and  2 

(Cl) 


By  Rule  3  we  can  "move"  this  reflection  coefficient  down  the  transmission 
line. 

s(0+)  -  e"j2^  s(i,-)  .  (C2) 


8  <£-) 


za -)  "  zq 
2  a-)  +  zo 


z%  -  z 
A  o 

z  +  z 
A  o 


We  can  now  compute  z(0+)  from  Rule  2. 


z(0+) 


z 

o 


1  +  s  a-) 

1  -  e-jB£  s(Jl-) 


(C3) 


With  the  continuity  of  Rule  1,  and  the  second  equation  of  Rule  2  we  obtain 
s(O-)  from  Equation  (C3) . 


s(0+) 


z (0+)  -  zT 
z (0+)  +  zT 


(C4) 


We  can  relate  these  variables  to  those  used  in  Section  3.1  of  the  text. 


z  =  z (0)  =  z (0+)  =  z(0-) 
A 

*A  =  s(H 

=  s ( 0+ )  =  m 

A 

s  =  s(0-) 

T 


In  this  section  of  the  appendix,  we  would  also  like  to  show 

v,  . (0-)  =  v  .  From  the  definition  of  traveling  waves  and  characteristic 

(+)  g 

impedance  we  have: 


v(0)  =  v  (0-)  +  v  (0- ) 
i  (0)  =  i  (0-)  +  i  (0-) 


(C5) 


C-  2 


But  from  "Figure  C.l,  we  see  that 


v 

g 


z  i  (0)  +  v  (0) 

T 

•  V1-II0‘)  +  v(+)  (°->  + 


(C6) 


Finally,  we  would  like  to  show  that  the  model  makes  sense  as 
the  length  of  the  transmission  line  approaches  zero,  i.e.,  in  the  absence 
of  any  line.  In  the  absence  of  any  transmission  line,  one  would  confute 
the  reflection  coefficient  of  the  antenna  as  seen  by  the  source  as 
T  *  (zA  -  zT)  /  (zA  +  z^,)  .  What  needs  to  be  shovn  is  that  s(O-)  approaches 
T,  as  £  approaches  0.  It  is  sufficient  to  show  that  as  £  approaches  0,  then 
z(0)  approaches  for  then  by  Rule  2  and  the  continuity  of  z(x) ,  s(O-)  will 
be  r  “  <ZA  “  zt^^za  +  ZTJ  *  As  *  approaches  0,  then  e~^2^  approaches  1  and 
from  Equation  (C3)  we  can  write 


lim 

£-K) 


z(0) 


z 

o 


1  +  s (£-) 
1  -  s  (£-) 


(C7) 


Note  from  Equation  (Cl)  that  s(£-)  does  not  depend  on  £.  By  combining 
Equations  (Cl)  and  (C7)  ,  we  find  the  limit  for  z(0)  as  £  approaches  0. 


C.2  Multiple  Transmission  Lines 


In  this  section  the  results  of  the  previous  section  are  extended 
for  the  case  of  N  transmission  lines  terminated  by  an  N+l  port  network 
(see  Figure  C.2).  First,  let  us  define  the  following  quantities  to  sim¬ 
plify  notation.  In  keeping  with  the  text,  we  shall  use  the  indices 
0,  1,  . . . ,  N. 


S (x)  matrix  of  s-parameters  relating  v,  .  (x)  and  v,  .  (x) 
—  — (+)  —  — (-) 


The  three  rules  stated  in  Section  C.l  are  restated  below  for 
multiple  transmission  lines.  (In  evaluating  these  or  any  following  expression 
the  components  of  x  must  be  all  in  the  same  region,  i.e.,  all  in  the 
transmission  lines,  all  in  the  load  ,  or  all  in  the  termination.) 

1)  v(x) ,  i(x)  and  Z(x)  are  all  continuous  in  x,  where 

v  (x)  =  v  (x)  +  v(_}  (x)  ,  i(x)  =  i  (x)  +  i.  (x)  and 

Z  (x)  is  defined  by  (x)  =  Z(x)^+}  (x)and 


Transmission  Lines 


Multiple  Terminated  Transmission  Lines 


>  •  V  V*  *  •  *  -  '7  "  •  ^ 


2)  Z(x)  =  [I  -  S(x)]”1  [I  +  S(x)]Zo(x) 


and  S(x) 


[Z  (x)  -  Z  (x)  ]  [Z  (x)  +  Z  (x)] 
—  o  ~  —  o  — 


where 


Z  (x) 
o  — 


Z  x.  >  £.  for  all  i 

Ail 

Z  0  <  x.  <  £ ,  for  all  i 
o  11 

Z_  x.<  0  for  all  i 

T  l 


3) 


S  (x^)  =  Dt^  -  2Sq)  s  (Xq)  -  x^) 


where  D(x) 


0 


±*li' 


x. .  <  H.  for  all 
Oi  —  l 


and  6 .  = 


1  V  . 

pi 


for  all  i 


As  in  Section  C.l,  we  start  at  the  load  Z 


Z(£+)  =  Zft 


(C9) 


Then  from  Rules  1  and  2  we  can  compute  S(£~) 


S  ( £-)  =  \Z(l~)  -  Z  U)][Z(£-)  +  Z  (£)J 

—  o—o 

*  l2A  -  Zo,|ZA  *  Zorl 


(CIO) 


By  Rule  3,  we  can  move  this  scattering  matrix  down  the  transmission  lines 
to  x  =  0+ . 


S(0+)  =  D(£)  S  (£,-)  D  (£) 


(Cll) 


From  S(0+)  and  Rule  2,  we  can  compute  Z(0+). 

Z(0+)  =  [I  -  S(0+)]-1(l  +  S(0+)]z 

o 


(C12) 


C-7 


With  Rule  2  one  can  obtain  S(O-)  since  7,(0-)  =  Z(0+)  from  Rule  1. 


S(O-)  =  [ Z (0+)  -  ZT]  [Z(0+)  +  z  ]‘ 


(C13) 


As  in  Appendix  C.l,  we  can  relate  these  variables  to  those  used  in  the 
text  in  Section  3.3. 


=  Z(0)  =  Z(0+)  =  Z(O-) 

=  S  (_£-) 

=  S (0+)  =  M+ 

=  S(O-) 


(C14) 

(C15) 

(C16) 

(C17) 


In  the  remaining  portions  of  this  appendix,  several  results 
useful  in  the  text  will  be  derived.  First,  Equation  (C6)  of  Appendix  C.l 
will  be  generalized  for  the  N-port  case.  We  have  by  definition 


V(0)  =  v  (0-)  +  v  (0-) 

i<o)  =  i(+)  <°->  +  i(_} (o-) 


By  Rule  1,  we  also  have 


i(0)  =  zT_1v(+) (0-)  +  zT_1v(_}  to-) 


But  from  Figure  C.2,  where  v  is  defined,  we  see  that 

-g 


v  =  Z  i(0)  +  v (0) 

-g  t-  - 

=  v(  +  )(0-)  -  v(_) (0-)  +  — (  +  )  (0-)  +  v(_)(0-) 


Therefore,  it  can  be  stated  in  the  general  N-port  case  that 


v  (0-)  =  Jjv 
-(+)  -g 


(C18) 


We  would  also  like  to  generalize  the  result  shown  in  Appendix 
C.l,  that  in  the  limit  as  l  (or  in  the  general  case  |&J  )  approaches  0, 
the  model  behaves  as  if  no  transmission  lines  were  present.  Mathematically 
we  need  to  show  that 


4 


t 


I 


lim 

|£|-0 


Z(0) 


z 


A 


First  let  us  define  some  normalized  quantities.  Assuming  that  z  ^  0 

on 

for  any  n,  then  Z  is  invertible.  In  addition,  since  Z  is  diagonal,  there 
o  o 

exists  a  matrix  Z  *3  such  that 
o 


and  where  Z 


-»3 


(Z  ) 


h-l 


(C19) 


(C20) 


I 


I 


Now  we  will  define  the  following  normalized  values  of  ZT  and  Zft. 


z'=z  i5z„z-'3  =  z-1zm=z„z  1 


T  o 


T  o 


(C21) 


since  Z^,  is  diagonal,  and 


/  -*3  -h 

Z  '  =  Z  Z„  Z 
A  o  A  o 


(C22 ) 


With  these  normalized  quantities,  let  us  rewrite  some  of  the  expressions 
previously  derived.  First  consider  S{Z)  from  Equation  (CIO). 

S  (£- )  =  (ZK  -  Z  )  (ZR  +  Z  5'1 
—  A  o  A  o 

*3#  *3  *3/*3  *3^3  —  1 

=  (Z  z'  Z  -  Z  Z  )  (z  z/z  +  Z  Z  ) 
oAo  OO  OAO  OO 


=  z  13  (z/  -  I)Z  15  Z  (Z  '  +  I)-1  Z  _l5 

o  A  o  o  A  o 

=  z  15  (Z /  -  I)  (z/  +  i)_1  z  _>5 

o  A  A  o 


(C23 ) 
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APPENDIX  G 


GRADIENTS  OF  PSUEDO-ANALYTIC  FUNCTIONS  AND 
APPLICATIONS  TO  COMPLEX  GRADIENT  CONTROL 


Appendix  G . 1 :  Theory  of  Psuedo-Analytic  Functions 

In  the  theory  of  complex  variables  analyticity  is  a  very  use¬ 
ful  property  of  a  function.  It  allows  one  to  treat  the  function  as  a 
real  function  of  a  real  variable  when  differentiating. 

There  are,  however,  many  important  functions,  e.g.  squared 
magnitude  and  conjugation  of  a  variable,  which  are  not  analytic.  In 
this  appendix  a  property  called  psuedo-analyticity  is  defined.  Many 
tools  for  the  study  of  this  property  are  derived.  Finally,  with  these 
tools  the  gradient  control  law  for  complex  variables  is  derived. 

From  single  variable  complex  analysis  we  know  that  if  a  func¬ 
tion  of  a  complex  variable  f(z)  =  u(x,y)  +  jv(x,y)  where  (z  =  x  +  jy) 
is  analytic,  then  there  is  a  well  defined  derivative  such  that  f  and  z 
follow  all  the  basic  rules  of  differentiation  of  a  real  function  with 
respect  to  a  real  variable. 


f(z) 


df 

dz 


8u  .  8v  8v 

8x  3  8x  ~  8y 


(Gl) 


A  necessary  and  sufficient  condition  for  analyticity  is  that  the  func¬ 
tion  satisfy  the  Cauchy-Riemann  equations. 

8u  8v  ,  8u  _  8v 

8x  *  8y  "  8y  8x  (G2) 

These  equations  can  be  written  more  compactly  as 

8f  ,8f 

8x  _:)8y  (G3) 

Now  let  f  be  a  complex  scalar  function  of  N  complex  variable, 
f(zx,  z2,...,zN)  =  u(z1,...,zN)  +  jv{z1,...,zN) ,  where  z±  =  x±  +  jy^ 
Also  let  f  be  analytic  with  respect  to  z^  over  the  range  of  all  possible 
values  of  the  other  variables.  We  may  then  define  the  complex  partial 
derivative . 


The  Inequalities  (FI  a  and  b)  can  be  rewritten  and  combined 


more  conveniently  as 


/  2 l 

'22  P2 


-  4m12  p 


1P2 


<  2  -  Re[m22p2] 


Note  if  Re[m22P2]  >  2,  then  the  right  side  of  Inequality  (F2) 
is  negative  and  no  value  of p  ^  can  satisfy  the  constraint.  If,  however, 
Re[m22p2^  <  2  then  exists  a  range  of  values  that  will  satisfy  the 

constraint.  One  such  value  is 

2 

m  D 

P,  =  22  *2 


which  makes  the  left  side  of  (F2)  equal  to  zero. 

We  have  shown  that  an  auxiliary  with  a  reflectivity  that  would 
normally  make  it  unstable,  can  sometimes  be  stabilized  by  the  pressence 
of  a  second  auxiliary  antenna. 


APPENDIX  F 


A  STABILITY  ANALYSIS  EXAMPLE 


Stabilizing  An  Unstable  Element 

Consider  the  2-element  auxiliary  array  illustrated  in  Figure  FI. 

It  is  assumed  the  self  coupling  factor,  return  loss,  of  the  first  antenna 
is  zero,  the  main  antenna  is  absent  or  has  =  0,  and  the  reflectivities 
have  dominant  single  stage  narrowband  filters  in  the  reflectivities  so 
that  Re  [all  eigenvalues]  <  1  is  both  necessary  and  sufficient  for  stability 

Consider  the  case  >  The  second  antenna  then  by 

itself,  or  with  p^  =  0,  would  be  unstable.  It  shall  be  shown  that  under 
certain  conditions,  Re[m22p2^  <  there  are  choices  for  p^  which  will 
stabilize  the  array. 

First  we  need  to  find  the  eigenvalues  of  MP  for  this  array. 


MP  = 

'°  "12" 

"Pi  0' 

"12P2  ‘ 

_m12  "22_ 

1 

CN 

a 

0 
_ 1 

_”12P1  "22P2_ 

Characteristic  Equation  =  det[sl  -  J*P] . 


det 


'"l2P2 
s  -»22P2 


'  3  '  “”22P2  '  "12  P1P2 

The  eigenvalues  are  the  roots  of  the  characteristic  equation. 


Eigenvalues 


m„„p„  ± 
22  k2 


V® 


2  ,  2 

22  P2  ~  4m!2  P1P2 


The  stability  constraint  applied  to  these  eigenvalues  yields 
two  inequalities. 


Re 


m22P2  +V^222p22  "  4m!22plp2 


<  1 


(Fla) 


and 


Re 


m22P2  -V“12'P2'  ‘  4m^_P 
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1P2  <  1 


(Fib) 


Note  the  quantity  within  the  brackets  is  necessarily  positive  because 
the  term  underneath  the  radical  is  the  sum  of  two  squares,  r  cos  26-1 
and  r2  sin  26. 

Substitite  cos  6  *  V(1  +  cos  28)/2  into  (E8)  and  square  both 

sides . 


2 

(1  +  cos2  ) 


+  1  -  2r2cos26  +  r2cos26  -  1 


or  equivalently 


1  +  r2  ?Vr4  +  1  - 


2r  cos26 


(E9) 


Squaring  both  sides  of  (E9)  now  yields 

2  4  4  2 

1  +  2r  +  r  >r  +l-2r  cos26 

2r2  >  -2r2cos28 
1  >  -cos26 


This  will  be  true  if  8  4  i\/2  +  mir  for  m  =  0,  ±1,  ±2,....  But  we  have 

the  constraint  1 6 1  <  TT/2j  therefore,  it  is  always  true. 

Thus,  we  have  shown  for  a  single  stage  narrowband  filter  in 

the  reflectivities  the  stability  constraint  RefA^]  <  1  is  sufficient. 

Necessity  is  shown  by  considering  the  frequency  where  ftjuM  =  1, 

and  noting  that  Re [A  ]  <  1  is  the  necessary  constraint  derived  in  Sec- 

n 

tion  4.2. 


E-5 


The  roots  of  (E4)  are  found  from  the  quadratic  equation. 


roots  -  -  |U  -  An)  ±V[ |(1  -  An)]2  -  Wq2 


He  want  to  show  the  two  values  of  (E5)  are  in  the  left  half  plane  if 

Re  [A  ]  <  1.  Let  x  *=  b(l  -  A  )/2ui  ,  then  equivalently  we  can  show 
n  no* 


Re  [  -  X  +Vx2  -  1  <  0  for  Re  [x]  >  0  . 


(E6)  is  really  two  inequalities,  one  for  each  choice  of  the  ambiguous 
sign.  These  can  be  combined  as  in  (E7)  . 


Re  [x]  >  (Re'oc  -  l| 


Let  x  =  reJ  ,  where  ]8|  <  ir/2.  Then  (E7a)  becomes 
r  cos  8  >  |ReVrV2^  -  l| 


>  |ReVr  cos  28  -  1  +  jr2sin  201 


(E7a) 


(E7b) 


At  this  stage  an  identity  for  complex  numbers  becomes  useful.  If 

0-0  +  jo  where  a  ,  a„  are  real,  then 

R  I  R  I 


/a  -  ± 


✓2  L 


+  °R  +  W  ^ 


I  2  2 

where  y  >=yOR  +  and  all  square  roots  on  the  right  side  aure  positive 
This  identity  can  be  easily  verified  by  squaring.  Using  the  real  part 
of  this  identity,  we  can  rewrite  (E7b) 


r  cos  6 


>  *  —  ~\jr^c os220  +  1  -  2r2cos26  +  r4sin220  +  r2cos20  -  l| 

/2  L 

“  f Vr4  +  1  -  2r2cos26  +  r2cos20  -ll*3 


N+l 


order  term  of  A(*)  will  vanish.  The  highest  order  term  will  become  D 
Since  it  was  assumed  that  D(s)  and  N(s)  had  no  common  factors,  D  (s) 
will  not  equal  zero  when  N(s)  does. 

Thus  for  any  zero  of  (El),  N(s)  ^  0.  The  zeros  of  (El)  will 
then  be  the  zeros  of  A.[f  ^ (s)  ].  The  zeros  of  this  factor  will  be  the 
values  of  s  for  which  f  (s)  equals  an  eigenvalue  of  M+P+.  Therefore, 
the  array  is  stable  if  and  only  if  all  zeros  of  the  following  N+l 
equations  are  in  the  left  half  plane. 


(s) 


f(s)  -  -0,  n  **  0,  ...»  N 


(E2) 


n 


and  X  is  the  nth  eigenvalue  of  M 

n  +  + 


Note  the  similarity  between  these  equations  and  that  for  posi¬ 
tive  feedback  root  locus.  They  are  the  same  if  f(s)  is  considered  the 
open  loop  transfer  function  and  X^  the  loop  gain.  The  difference  is 

that  in  (E2)  X  is  complex  instead  of  positive  real, 
n 

Finally,  it  was  mentioned  earlier  that  if  a  single  stage  nar¬ 
rowband  filter  is  used  in  the  reflectivity,  the  necessary  stability  con¬ 
dition,  ~Re[Xn]  <  1,  becomes  both  necessary  and  sufficient.  This  shall 
now  be  proved. 

What  needs  to  be  shown  is  for  f(s)  given  in  (E3)  ,  the  equations 

of  (E2)  will  have  no  right  half  plane  zeros  if  Re[X  ]  <  1. 

n 


f(s) 


sb 


2  ,  2 
s  +  sb  +  a) 

o 


[s  - 1  -  - 1 +  ’  V“„"  -  (D'J 


(E3) 


Here  <*>o  is  the  center  frequency  of  the  filter  and  b  is  the  3  dB 


It  is  assumed  U)q  >  b/2  and  indeed  b  is  small  enough  that 


bandwidth. 

M+(s)  does  not  change  significantly  across  the  filter  passband.  The  equations 
of  (E2)  become 


sb 


2  2 
s  +  sb  +  to 

o 


1 

X 


n 


or 


s2  +  8b(l  -  Xn)  + 


(E4) 


E-3 


nth  Auxiliary 
Antenna 


Figure  El  Reflectivity  with  Narrowband  Filter 
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APPENDIX  E 


SUFFICIENCY  OF  STABILITY  CONDITION  WITH 
NARROWBAND  FILTER  IN  REFLECTIVITIES 

Let  us  apply  the  general  stability  criterion  of  Section  4.2 

to  the  case  of  narrowband  filters  within  the  reflectivities  (see  Figure 

El)  .  It  is  assumed  that  the  narrowband  filter  is  narrow  enough  so  that 

M+(s)  is  constant  over  the  passband  and  is  approximated  by  M+.  Also 

assume  that  the  filters  within  the  reflectivities  are  all  the  same,  so 

that  the  nth  entry  on  the  diagonal  of  P  (s)  is  p  f (s) .  Here  p  is  a 

+  n  n 

constant  and  f(s)  is  the  same  for  all  reflectivities.  Let  f(s)  ■  N(s)/D(s) 
where  N(s)  and  D(s)  are  polynomials  with  no  common  factors.  In  this 
case  the  det[M+(s) P+(s)  ]  «  p^  ...  PNfN+1(s)det[M+]  is  a  minor  of  order  N  +  1. 
Its  denominator,  D»+1,  is  the  lease  common  denominator  of  all  the  minors. 

Let  P+(s)  =  P+f  (s)  ,  where  P+  with  no  argument  is  the  diagonal  matrix  of 
the  Pq,...,p^.  Then  the  system  is  stable  if  and  only  if  the  following 
polynomial  has  no  right  half  plane  zeros  {16,  p.376]  , 

M+l 

D  (s).det[l  -  M  P  f  ( s)  ] 

+  + 

=  DN+1(s) -fN+1(s) -detll  •  f_1(s)  -M+P+] 

-  NN+1(s)-det[l  •  f_1(s)  -  M+P+] 

«=  NN+1(s)  A[f_1 (s) ] 

«  NH+1{s)  A[D(s)/N(s)J  (El) 

where  A  (•)  is  the  characteristic  or  eigenvalue  equation  of  M+P+. 

Note  this  is  a  polynomial  because  A(*)  is  a  polynomial  of  order 

N  +  1  and  all  the  N(s) 's  in  the  denominator  of  f  1(s)  will  be  cancelled 
N+l 

by  the  factor  N  (s)  out  front. 

Let  us  first  see  what  happens  to  the  polynomial  (El)  at  a  zero 
of  N(s) .  Because  of  the  preceding  factor  of  NN+^(s)  all  but  the  highest 


E-l 


From  the  expression  for  P+  and  Identities  A. 6  and  A. 7  of  Appendix  A,  it  is 
seen  that 


(I  +  P+)  -  2Zt(Zt  +  zoI) 


(D.4) 


(I  -  P+>  «  2zo(Zt  +  zoI) 


(D.5) 


Substituting  Equations  (D.4)  and  (D.5)  into  Equation  (p.3)  obtains 


ST  +  P+  “  <ZT  +  ZoI}  Vt  +  Zoirl2Zo(ZT  *  2oI,_1m+ 

•  (I  -  P+M+)_1(Zt  +  ZQI) _1 


(D.6) 


Since  2zq  is  a  scalar  its  placement  is  arbitrary.  Also  since  (ZT  +  zQI) 
and  2Zt  are  both  diagonal,  they  commute.  Therefore,  let  us  rewrite 
Equation  (D.6)  as  below. 


Sm  +  P, 
T  + 


-  2Zt(Zt  +  zoI)  (ZT  +  zoI)"1(ZT  +  zoD_1M+ 

*  «  “  W'l2zo(ZT  +  V5’1 

-  2Zt(Zt  +  zoI)_1M+(I  -  P+M+)_12zo(Zt  +  zoI)_1 


(I  +  P+)M+(I  -  P+M+)  (I  -  P) 


(D.7) 


It  only  remains  to  be  shown  that 

m+(i  -  p^r1  -  (i  -  M+P+)-\ 

and  the  equivalence  of  Equations  (3.16)  and  (3.21)  will  have  been  proved. 

M+(I  -  P+M+)_1  «=  [(I  -  M+P+)_1(I  -  M+P+)]m+(I  -  P+M+)_;L[(I  -  P+M+)(I  -  P^)'1] 

-  (I  -  M+P+) _1  [(I  -  M+P+)M+]  (I  -  P+M+)_1 

-  (I  -  M+P+)"1(M+  -  M+P+M+) (I  -  P+M+)_1 

-  (I  -  M+P+)_1  [m+(I  -  P+M+)]  (I  ~  P+M+)~1 

-  (I  -  M+P+) 


’»  m  *  -  *  •  *  •  * 


\  *■  *.  /.  *-  *• 


APPENDIX  D 


EQUIVALENCE  OF  EQUATIONS  (3.16)  AND  (3.21) 

This  appendix  will  show  that  the  two  forms  for  derived  in 
Section  3.3  are  equivalent.  These  expressions  are  repeated  below. 

ST  *  (ZT  +  zoI)  (gA  *  P+}  U  "  Va)1(ZT  +  V)_1  (3*16) 

and 

*  (I  +  P+)(I  -  M+P+)_1M+(I  -  P+)  -  P+  (3.21) 

Note  that  it  was  shown  in  the  text  that  M  •=  S,  and  by  an  arbitrary  choice 

+  A 

we  shall  use  the  former  expression.  It  will  be  convenient  to  add  P+  to  both 
equations  and  show  the  equivalence  of  the  resulting  two  forms  for  ST  +  P+ . 

ST  +  P+  ~  (ZT  +  ZQI) (M+  -  P+) (I  -  P+M+)_1(Zt  +  zoI)_1  f  P+  (D.l) 

However,  we  know  from  the  text  that  P+  «  (ZT  -  zqI) (ZT  +  zqI)  But  since 

all  the  matrices  in  this  expression  sure  diagonal,  we  can  also  write 

P+  “  (ZT  +  zoI)P+(ZT  +  zoI)"1  •  (D,2) 


Substituting  Equation  (D.2)  for  the  last  term  of  (D.l)  yields 

Sm  +  PL  -  (Z_  +  z  I) (M.  -  P  ) (I  -  P  M  )”1{Z_  +  z  I)_1 

T  +  T  o  +  +  ++  T  o 

+  (ZT  +  zoDP+(ZT  +  zoI)_1 

"  (ZT  +  zoI)  [ (M+  -  P+)  (I  -  P^)"1  +  P+]  (ZT  +  z^)"1 

“  <ZT  +  ZoI}  [(M+  '  V  +  V1  -  W]  U  •  P+M+)_1{ZT  +  ZoI)_ 

-  (ZT  +  zoI)(I  -  P+2)M+(I  -  P+M+)-1(Zt  +  ZoI)_1 

-  (ZT  +  zqI) (I  +  P+)(I  -  P+)M+(I  -  P+M+)_1(Zt  +  zoI)_1  (D.3) 


D-l 


Now  let  us  determine  v^+j  (0+)  in  terms  of  the  source  voltages  From 

Equation  (C18)  we  see  that  v  (0-)  =  »j  v  .  v, , .  (0-)  is  now  the  traveling 
voltage  wave  from  the  source  incident  on  the  termination  transmission 
line  interface.  1°+)  is  composed  of  the  component  transmitted 

through  the  interface  by  v^+j  (0-) ,  which  is  (I  -  P+)  v^  (0-)  ,  and  by 
itself  first  reflected  off  the  antennas  and  then  off  the  terminations, 
P|M|V( ^  (0+) .  As  an  equation  this  becomes 

v(+) (°+)  =  (I  -  P+)v(+)(0-)  +  P+M+v(+) (0+)  (C30) 

Solving  Equation  (C30)  for  v  (0+)  we  obtain 

v(+)(°+)  -  (I  -  P+M+)-1  (l  -  P+)v{+) (0-) 

“  Ml  -  P+M+) _1  (I  -  P+)v 


Finally  we  may  substitute  this  into  Equation  (C29) . 


!(£)  =  *jZ  1  (D  -  D_1  M  J  (I  -  PM)  _1  (1  -  P  )v 

-  O  +  *f  *r  +  — 


if  all  the  transmission  lines  have  the  same  impedance,  then  Z 


9 

-1 


(C31) 


U/zo)I 


and  (I  -  P+M+)  =■  Q+  .  If  the  lines  are  short  compared  to  a  wavelength, 
then  D  «*.  I.  Under  these  assumptions  the  port  currents  become 


Therefore,  M+  is  symmetric  if  S(H  is  symmetric.  From  Equation  (C8) 


SU-)  -  IZA  -  zo]  [zA  ♦  z^'1  . 


From  Identity  A. 5  of  Appendix  A,  we  see  that  S(4-)  is  symmetric  if 
T 

Z„  z  is.  But  Z„  is  symmetric,  and  if  z_  “  *  for  all  n,  then 
A  o  A  on  o 


Z„  Z  =  z  z. 
A  o  o  A 


which  is  symmetric.  Therefore,  if  all  the  ZQn  are  equal,  then  M+  is 
symmetric. 

Finally,  an  expression  for  the  antenna  port  currents,  i^fc) ,  will  be 
derived  in  terms  of  transmission  line  parameters.  By  definition  and  from 
Rule  1,  we  obtain  an  initial  expression  for  i_(i)  • 


i(4)  =  i +  i  U-) 

=  Z  -V,  U-)  -  Z  y  (£-) 
o  — (+)  —  o  — (-)  — 


(C27) 


For  simplicity  of  notation  let  the  delay  matrix  D(_£)  be  denoted  as  D. 
Noticing  that 


and 


(£-> 

a-) 


■»<*)  <°+> 

D~Xv  (0+) 


Equation  (C27)  can  be  rewritten  as  below. 

i  (4)  =  z;1  Dv(  +  )  (0+)  -  D_1v  }  (0+) 


(C28) 


But  we  also  know  that  (0+)  =  M(v ^ | ^  (0+)  from  the  definition  of  S(0+) 

and  Equation  (C16) .  This  can  be  substituted  into  Equation  (C28)  to  obtain 


iU) 


(D  -  D'1M+)v(  +  ) (0+)  . 


(C29) 


C-ll 


Using  Identity  A. 3  of  Appendix  A,  we  can  put  Equation  (C23)  into  a  form 
% 

that  will  turn  out  to  be  more  convenient.  (It  was  the  use  of  this 
identity  that  required  the  normalization  procedure. 

S (£-)  =  Z  **  (z/  +  I)”1  (Z  '  -  I)  z  15  (C24) 

—  o  A  A  o 

However,  as  |_£|  approaches,  D(£.)  approaches  I.  Therefore,  from  Equation  (Cll) 

S(£-) 

Z  15  (Z  '  +  I)"1  (Z  '  -  I)Z  -1*  (C25) 

O  A  A  o 

We  can  now  show  the  limit  as  |£j  approaches  0  of  Z(0)  is  Zft.  From 
Equation  (C12)  and  noting  that  Z(0)  =  Z(0+),  we  obtain 

lim  Z (0) 

| 


ZA  (C26) 

Now  it  will  be  shown  that  if  all  the  transmission  lines  are  of 

same  characteristic  impedance,  then  M+  =  =  S(0+)  is  symmetric,  i.e. , 

=  M  . 

+ 

M+  =  S(0+)  =  D(£)S(£-)D(£) 

=  dT(£)s(£-)d(£) 
because  D(x)  is  diagonal  and  hence  symmetric. 


the 


=  lim  [I  -  S (0+) ] _1 [ I  +  S(0+)]Z 


£  ho 


'  11  '  ZB  'V  +  I|"1(Za'  -  +  Z*  (ZR»  +  I1_1(ZA'-  «Z0',]Z0 


*  Zo’  [I  '  IZa'  +  -  HJ‘VZo  [I  +  <Z»  *  I,'1(ZA  -  I))Zo"’Zo 


=  z  **  [z/  +  I  -  z'  +  I] _1  (Z  '  +  I)  (Z*  +  I)_1[Z  ‘  +  I  +  z‘  -  I]z  13 


-h 


O  1  A 


-  0  (2I)'1<2Zs'>Z0‘! 


lim  S(0+)  = 

UN 


c-10 


(G4) 


WI  ■  „■ 1  L"  J 


■m-r-r-T 


3f  3u  ,3v  jW  _  .  3u 

3zi  =  9xi  3Sxi  =  9yi  ~  D9Yi 


If  f  is  analytic  for  all  the  variables  z^, 
a  complex  gradient 


z  ,  then  one  can  define 
N 


(G5) 


Consider  a  function  g(z)  =  r(x,y)  +  jq(x,y) ,  which  is  not 
analytic,  but  whose  real  partials  3r/3x,  3q/3x,  3r/3y  and  3q/3y  exist 
throughout  the  region  of  interest.  Such  a  function  is  defined  to  be 
psuedo-analy tic .  For  such  a  function  it  is  useful  to  define  the 
psuedo-derivative  g'(z)  as  in  (G6)  . 


jiSL  « 


3r 

3x 


.3r  _  33. 

Jdy  3y 


(G6) 


Note  from  (G3)  we  see  that  the  psuedo-derivative  of  an  analytic  function 
is  zero.  The  computation  of  the  psuedo-derivative  can  be  a  test  for 
analy ticity . 

Psuedo-analyticity  can  be  extended  to  multivariate  functions 
as  analyticity  was.  Consider  a  complex  scalar  function  of  N  complex 
variables  glz^, . . . , z^)  =  x(.Zy  . . .  ,z^)  +  jq(z1» . . .,zN> .  If  the  partial 
derivatives  3r/3x^,  3r/3y^,  3q/3x^  and  3q/3y^  exist  for  all  possible 
values  of  the  other  variables,  then  g  is  said  to  be  psuedo-analy tic  in 
z^  and  one  can  define  the  psuedo-partial  derivative  as  in  (G7) . 


A 


la. 


i£  +  -jla  _  la.  +  l£ 

3xi  33xi  dyi  dy± 


(G7) 


Similarly,  if  g  is  psuedo-analytic  for  all  the  variables  z^,...,^,  then 
one  can  define  the  psuedo-gradient  as  in  (G8) . 


G-2 


(G8) 


A 

Again,  note  for  an  analytic  function  f,  Vf  “  0. 


To  ease  the  cumbersome  notation,  let  us  define  the  operators 


V  (•)  and  V  (•) • 

x  y 


vx(o 


(G9) 


With  these  operators  Equation  (G8)  can  be  rewritten. 

A 

Vg  ■-  Vxg  +  jVyg 

*  Vr+jVq-Vq+jVr  (G10) 

x  x  y  y 

Note  Vx(g*)  =  Vxr  -  jVxg  =  (Vx9>*  where  (  )*  denotes  conjugation.  Simi¬ 
larly,  Vy  (g*)  =  (Vyg)*.  From  Equation  (G3) ,  we  see  the  gradient  for 
an  analytic  function  f  is  easily  written  as  in  (Gil) . 


Vf  =  V  f  =  - jV  f 

x  y 


(Gil) 


Finally,  we  define  the  operator  V  ,  the  conjugate  psuedo-gradient. 


in  (G12) 


V  g  -  Vxg  -  jVyg 


(G12) 


Note  for  an  analytic  function  f,  the  conjugate  psuedo-gradient  is  not 
zero  but  twice  the  analytic  gradient. 


V*f  -  V  f  -  jV  f  “  2Vf 
x  y 


(G13) 


G-3 


He  will  prove  three  theorems  useful  when  working  with  psuedo- 


gradients . 

Theorem  1:  The  conjugate  of  the  psuedo-gradient 
gradient  of  the  conjugate,  i.e., 


(Vg) 


* 

V  g 


proof:  (Vg) 


(Vxg  +  jVyg) 


(Vxr  +  jVxq  -  vyq  +  jVyr) 


Vr-jVq-Vq-jVr 

X  x  y  y 


*  * 

V  g  -  jV  g 
x  J  y 


*  * 

V  g 


Theorem  2:  For  two  psuedo-analytic  functions  of 
the  product  rule  of  derivatives  holds,  i.e.. 


V(g1g2)  =  gjVg2  +  g2Vg  . 


is  the  conjugate  psuedo- 

(G14) 


N  variables,  g.^  and  g2> 


(G15) 


proof:  V(g1g2) 


Vx{gig2)  +  jVy(gig2) 


Theorem  3:  For  an  analytic  function  f, 

V|  f  | 2  -  2f  (Vf)  *  (G16) 

proof:  V | f | 2  «  V(ff  ) 

=  fVf  +  f  Vf  ,  from  (G15) 

*  * 

«  f(Vf  )  +  0  »  from  note  after  (G8) 

# 

=  f(V  f)  ,  from  (G14) 

=  f (2Vf ) *  ,  from  (G13) 

=  2f  (Vf)  * 


Appendix  G . 2  s  Complex  Gradient  Control 

Consider  a  system  (e.g.  a  parasitic  array)  that  is  controlled 
by  N  complex  state  variables  z-,...,z„  (e.g.  the  reflectivities).  It 
is  desired  to  minimize  some  positive  definite  scalar  quantity,  V (z) 

(e.g.  power  received  at  the  main  antenna)  .  A  function  is  positive  defi¬ 
nite  if  it  is  nonnegative  real  for  all  values  of  its  arguement  and  is 
zero  only  for  an  arguement  of  zero .  Note  that  V (z)  is  not  analytic  in 
any  of  its  state  variables.  If,  however,  v (z)  is  psuedo-analytic  in 
all  of  its  state  variables,  one  can  employ  the  complex  gradient  control 
algorithm  defined  by  (G17) . 

A 

z  =  -k  V  V(z)  ,  where  k  >  0  ,  z  = 

A 

and  V  is  the  psuedo-gradient  defined  in  the  previous  section. 

In  applying  this  control  law  one  equates  real  and  imaginary 
components  of  both  sides.  Let  z_  *  ic  +  j£  where  x  and  £  are  real  vectors. 
Then  (G17)  implies  the  following  two  real  control  laws. 


x  -  — k  Re  [V  V  (z)  ] 
i  «  -k  IM  [V  V(z)] 


(G16) 


(G19) 


From  Equation  (G10) ,  we  can  rewrite  (G18)  and  (G19) . 


x  «  -k  Re  ^  VxV  (x,yr)  +  jVyV  (x,^) 

-  -k  V^V  (x,^r) 

i  “  -k  Ini  (x>iL)  +  jV^V  (x,Z) 

-  -k  VyV  (x,^) 


(G20) 


(G21) 


The  quantity  W{x,y)  separates  into  real  and  imaginary  components  as 

above  because  VXV  and  VyV  are  purely  real  since  V  is  positive  definite. 

But  Equations  (G2Q)  and  (G21)  simply  represent  the  familiar 

real  gradient  control  law  for  the  2N  real  variables  x  x  ,  y,  ,...fy 

1  N  1  N 

The  complex  control  law  of  (G17)  is  just  a  compact  way  of  writing  this 
for  complex  variables. 

Let  us  now  compute  dV/dt  from  the  real  multivariate  chain 

rule. 


dV 

dt 


(V  V)Tx+  (VV)Ti 

x  y 


-  Re 


=  Re 

**  Re 


(V  V  +  jVV)*T(£  +  j£) 

*  y 

(W)  *T2  1 

A  ^  ^  “I 

(W)  (-k  W)J 


-k|Vv| 


2 


(G22) 


We  see  that  dV/dt  is  always  nonpositive.  Thus  the  system  is  stable  in 
a  control  sense.  This  property  of  dV/dt  is  negative  semi -definiteness. 
If  dV/dt  or  equivalently  |W|  is  zero  only  when  z  «*  0,  then  dV/dt  would 
be  negative  definite  and  the  system  would  be  assymptotically  stable  in  a 
control  sense.  This  means  there  are  no  nonzero  local  minimum  of  V  for 
the  control  law  to  get  trapped  in  and  it  will  always  drive  V  to  zero. 

A  useful  special  case  is  when  V(z)  =  |fU)  |2  where  f(z)  is  an 
analytic  function  over  all  the  state  variables .  We  can  then  apply 
Theorem  3  of  the  previous  section  to  Equation  (G17) . 

z  -  -2k  f(Vf)* 


G-6 


(G23) 


This  is  the  form  of  the  complex  gradient  law  we  shall  find  most  useful 
in  our  analysis. 

Lastly,  consider  the  example  of  a  conventional  adaptive  array, 
as  in  Figure  1  of  the  text.  The  array  output  is  given  by  (G24) . 

e  *  yQ  +  (G24) 

where  y^  is  the  signal  at  the  main  antenna,  £  is  the  vector  of  signals 

at  the  auxiliary  antennas  and  w  is  the  vector  of  the  adaptive  weights. 

2 

It  is  desired  to  minimize  |e|  .  We  can  apply  (G23) . 
w  =  -2k  e(Ve) * 


=  -2k 


from  Identity  A. 17 


(G25) 


This  is  the  well  known  least  mean  square  (LMS)  algorithm. 

Appendix  G.3:  Gradient  Control  of  Active  Complex  Parasitic  Terminations 

Let  us  apply  this  theory  to  the  adaptive  control  of  parasitic 

elements.  Consider  Equation  (G17) .  Our  positive  definite  function  is 
2 

|y  |  ,  the  power  incident  on  the  termination  of  the  main  antenna.  To 

apply  (G17)  we  need  state  control  variables  which  are  zero  when  |yQ|  is. 

Let  these  variables  be  z  *  p  -  p  .,  where  p  is  the  nth  reflectivity 

n  n  no  n 

and  Png  is  the  value  of  this  reflectivity  necessary  for  a  null.  From 

Chapter  2,  we  know  I y _ 1 2  as  a  function  of  the  reflectivities  p  .  Let 

0  2  2  ^ 

us  signify  this  by  | yQ (£)  |  *  |Yq<£  +  Pq) |  where  £,  z,  and  are 


respectively  vectors  of  the  p  ,  z  ,  and  p  _. 

n  n  no 


Lastly,  let  V  (  )  repre- 
~  P 


sent  the  gradient  with  respect  to  the p  n  and  )  and  (  )  be  the 

gradient  and  psuedo-gradient  with  respect  to  the  z  .  Equation  (G17)  now 

n 

becomes  (G26) . 


z  -  £  -  -kVz  |yo  (z  +  £o) 


-2k  y  (z  +  p  )  [v  y  (z  +  p  )]* 
0  0  l  Z  0‘  0  J 


(G26) 


Note  with  the  notation  of  Appendix  A, 
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We  can  apply  the  chain  rule  of  A. 20  and  obtain  (G27). 


V0‘*  +  fi0)  ■  ( Vsi)T 


Oy/  £  >  l3£/3i 


(G27) 


But  3g/3z  is  the  identity  matrix.  Equation  <G27)  then  becomes 


V<M  +  "  {3yc/  £T)T 

■  *y0/*£. 


Vo<£) 


(G28) 


This  has  singly  been  an  application  of  the  chain  rule  to  the  set  of  trans¬ 
lated  variables  z  “  p  -  p  _ .  The  result  allows  us  to  state  the  control 
n  Kn  nO 

law  in  terms  of  the  reflectivities  without  the  need  to  know  the  values 


necessary  for  nulls. 


£-  -2k  y0<£>[Vo(£)]  * 


(G29) 


The  vector  £  is  not  used  in  computation,  but  the  reflectivities  are 
arranged  into  the  diagonal  matrix  P.  To  use  Equation  (G29)  we  must 
separate  it  into  N  equations  of  the  form. 


,  n  «  1,  . .  • ,  N 


(G30) 


We  must  now  compute  3yrt/3p  .  The  relationship  of  y  and  p  is  expressed 

u  n  ■*0  n 

by  Equation  (2.3). 

yo  "  yso  +  pomooyo  +  2oTp* 


Z  -  Zs  +  P^o  +  MP Z 


(2.3) 
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Let  us  differentiate  both  sides  of  these  equations  using  differentiation 
formulae  A. 16  and  A. 17. 


.  p  „  .  j.t/e  *  p|x- 

3pn  00  3pn  -0  I  nn%  3pf 


(G31a) 


where  E 


3yo  (  \ 

3pn  +  M  (w  +  p  3pn  ) 

i-n,  j-mj 


(G31b) 


Solve  (G31b)  for  -r*—  and  substitute  into  (G31a) . 
”Pn 


°  +mT  I 
=<) 
n 


0  m  x — 

Ho  oo  3p 


/  T  - 

p  I  m _ +  m  PQ 

H0\  00  —0  * 


V  +  ^  (pcA)  3p^  + 


1  \^0 

^3pn 


+  ^  (l  +  P2"1m)  Enn* 


(G32) 


where  Q  -  I  -  MP 


Note  the  existance  of  Q  is  necessary  for  RF  stability.  Using  matrix 

T 

identity  A. 8  with  B  =  I,  V  =  M,  U  =  P,  and  realizing  MP  =  (PM)  ,  ,  one  can 
-1  -1  T 

show  (I  +  PC  M)  »  (Q  )  .  Lastly,  define  e  *»  vector  of  all  zeros 

^  T 

except  the  nth  element  which  is  1.  Then  E  =  e  e  _  ,  and  E  _y  =  e  y  . 

nn  — xr-n  nn*-  — n  n 

We  can  now  solve  (G32)  for  3y./3p  . 

o  n 


T  -l.T 

2q  (2  ) 


1  "  po  (m00  +  2oTp2_12o) 


T  -1 
e  0 

— n  — 0 


1  -  p0  (”oo  +  SoTpe'  s0) 


a  y  ,  where  this  defines  a 
n  n  n 


(G33) 
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Equation  <G33)  can  be  substituted  into  the  control  law  of  (G30) . 
Pn  =  '2k  *0  (Vn}*  <634) 


Compare  this  expression  to  the  LMS  control  law  for  a  conventional  adaptive 

* 

array  in  (G25) .  They  are  the  same  except  for  the  factor  ,  where 

corresponds  to  w  ,  y.  to  e  and  y  to  the  signal  at  the  nth  antenna.  The 
n  u  n 

physical  significance  of  a  is  discussed  in  the  test. 
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APPENDIX  H 


DERIVATION  OF  GRADIENT  CONTROL  WITH  PILOT 


From  Appendix  G,  Equations  (G30)  and  (G31)  state  the  gradient 
control  law  for  an  adaptive  parasitic  array. 

* 


-2kyo‘Vn’ 


for  k  >  0 


where  a 


n 


„  T  -1 
e  Q  m 
— n  *  z~q 


and  e 


1  ■  p0  (n00  +  2oTp2~1S0) 


(G31) 

nth  element 
(G30) 


In  implementing  gradient  control  it  will  be  necessary  to  estimate  the 

factor  a  .  It  will  be  shown  in  this  appendix  that  a  pilot  signal  sent 
n 

up  the  main  antenna  will  be  received  at  the  nth  auxiliary  antenna  multi- 


T  -1 


Re¬ 


plied  by  the  factor  a^,  under  the  condition  pQ  =  0. 

From  the  expression  for  a  above,  we  see  that  a 

n  n 

Now  consider  a  pilot  signal,  y^Q,  injected  upward  at  the  main  antenna 
(see  Figure  HI)  .  This  upward  voltage  traveling  wave  couples  to  the  N 
auxiliaries  with  the  coupling  factors  of  n^.  This  creates  an  initial 
downward  vector  of  traveling  waves  of  m^y^.  From  Equation  (2.5b)  the 
final  vector  of  downward  traveling  voltages,  is  found  to  be  Q  1®0^po' 


when  p0  *  0.  The  voltage  traveling  downward  at  the  nth  antenna,  y 


is  eTQ-1m  v  . 
n  -0  p0 


pn 


rp0 


T  -1 
e  Q  m 
— n  —o 


(HI) 


given  pr 


It  can  be  shown  for  the  case  pQ  f  0,  through  Equations  (2.3 

a  and  b) ,  that  the  ratio  y  /y  _  is  in  general  given  by  (H2) . 

pn  pu 


“ 

- 

1 

T 

e 

Q  +  (Po 

\  T 

\m_m  P 
1=0=0 

_ 1  ■  moopo _ 

— n 

V1  ~  m00p0 , 

'po 


This  does  not  equal  a  .  Thus  the  constraint  p. 

n  0 


-1 


m 

=0 


0  is  clearly  needed. 


(H2) 


H-l 


END 
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